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In  this  dissertation,  we  present  an  analysis  of  continuous  review 
models  of  a two-echelon  inventory  system  for  recoverable  items.  The 
system  consists  of  a depot  and  a set  of  bases.  Primary  demands  occur 
at  the  bases  for  one  or  several  units  at  a time.  It  is  assumed 
that  demands  arrive  in  a Poisson  manner.  Upon  arrival  of  a demand  for 
certain  units,  a like  number  of  failed  units  are  turned  in  at  the 
base.  An  inspection  of  the  failed  units  is  carried  out  to  decide 
whether  the  units  will  be  repaired  at  the  base  or  at  the  depot  or  will 
be  removed  from  the  system  in  case  repair  is  not  economical.  The 
bases  use  an  (s-1,  s)  policy  for  procurement  of  serviceable  units 
from  the  depot,  and  the  depot  uses  an  (s,S)  policy  to  procure  from 
the  external  supplier.  Demands  in  an  out-of-stock  situation  are 
backlogged.  It  is  assumed  that  all  the  locations  have  infinite  repair 
capacities  and  repair  and  procurement  lead  times  are  constant. 

A common  problem  in  inventory  management  is  to  specify  the  policy 
parameters  that  will  minimize  expected  cost  per  unit  time  for  operating 
the  system  subject  to  constraints  of  certain  performance  measures. 

To  formulate  such  a problem  we  must  find  the  stationary  distributions 
for  inventory  position,  on-hand  inventory,  backorders  and  in-repair 
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inventory.  Our  main  objective  is  to  find  exact  expressions  for  these 
distributions. 
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The  investigation  begins  with  an  extensive  analysis  of  a single 
location  system.  The  procurement  policy  is  a continuous  review  (s,S) 
policy.  The  inter-arrival  times  between  successive  requisitions  are 
independent  and  identically  distributed  random  variables.  The  system 
experiences  two  types  of  demands  - recoverable  and  non-recoverable. 

The  two  demand  processes  may  be  independent  or  dependent.  For  the 
inspection  of  failed  units,  two  models  - batch  and  unit  - are  considered. 
In  the  batch  model,  the  entire  batch  is  either  recoverable  or  non- 
recoverable,  whereas,  in  the  unit  model  each  unit  in  a batch  is 
inspected  independently.  The  special  cases  of  compound  Poisson  demands, 
(s,  nQ)  procurement  policy,  complete  recoverability  and  cmmplete  non- 
reccverability  are  also  considered. 

For  the  two-echelon  system  we  first  consider  the  case  where  demands 
at  the  bases  occur  for  a single  unit  at  a time.  The  approach  is  then 
applied  to  a general  situation  where  demands  at  the  bases  are 
random.  Both  the  batch  and  unit  inspection  models  are  considered. 

For  the  case  when  there  are  no  condemnations  of  the  item,  results  are 
compared  with  the  METRIC  model.  The  METRIC  model  provides  a simple 
but  approximate  expression  for  the  probability  distribution  of  system 
backorders.  The  comparison  indicates  that  there  is  a considerable 
discrepancy  between  the  METRIC  results  and  our  results  when  the  depot 
spare  stock  is  low  or  when  a major  proportion  of  the  repair  is  done 

t 

at  the  depot. 
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CHAPTER  I 


IivTU' •3UCTIQU 


This  study  is  devoted  to  the  analysis  of  recoverable  (repairable) 
item  inventory  systems.  Upon  failure,  a recoverable  iter  is  returned 
to  the  source  of  supply  (inventory  point)  where  a decision  is  made 
either  to  remove  (condemn)  the  item  from  the  inventory  system  cr  to 
perform  repair  on  it  in  order  to  restore  it  to  a serviceable  condition  . 

r 

The  decision  to  repair  or  condemn  the  failed  item  is  based  on  the 
degree  and  the  nature  of  failure,  the  repair  facilities  available,  and 
the  economics  involved.  Once  an  item  is  designated  as  recoverable,  it 
is  presumably  more  economical  to  repair  the  item  than  it  is  to 
dispose  of  it  and  replace  it  with  a new  item. 

Most  inventory  systems  consist  of  consumable  (non-recoverab.le) 
items  that  are  predominantly  low  in  cost.  In  many  large-scale 
industrial  activities,  military  organizations,  for  example,  a large 
proportion  of  the  inventory  investment  is  in  recoverable  items, 
although  precentage-wise  most  items  are  consumable.  Hence  management 
of  recoverable  item  inventory  systems  from  both  the  design  and  control 
viewpoints  is  important. 


1 . 1 Pecoverab.l o Item  Inventory  Systems 

A typical  recoverable  item  inventory  system  consists  of  customers, 
inspection  and  repair  facilities,  supply  (inventory)  points,  and  an 
external  supplier  (manufacturer’) . 

Customers  generate  the  primary  demands  on  the  system.  While 
placing  requisitions  for  the  replacement  of  one  or  several  units, 
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customers  turn  in  a like  nurfc cr  of  failed  units.  In  view  of  tho 

V. 

description  presented  earlier,  the  system  enporienc.es  two  types  of  demand 
raoov arable  and  non-recovcrab.le.  Depending  on  the  processes  generating 
the  failures , the  two  demand  processes  can  be  either  (1)  independent 
cc  (2)  dependent.  Demand  processes  are  independent  when  there  are  two 
independent  processes  generating  the  two  types  of  failures.  Although 
dependence  of  the  demand  processes  may  arise  in  different  ways,  we 
shall  limit  the  consideration  to  the  case  where  a single  failure  process 
results  in  both  types  of  failures.  Thus  for  independent  demand  processes 
there  are  two  types  of  customers  from  two  independent  sources  whereas 
for  dependent  demand  processes  customers  arrive  from  a single  source  . 

We  assume  that  customers  arrive  from  an  infinite  population. 

Also,  the  inter-arrival  times  of  the  customers  and  the  number  of  units 
demanded  upon  an  arrival  are  assumed  to  be  random  variables  with  finite 
means  (known). 

Upon  arrival,  a batch  of  failed  units  is  inspected  to  classify 
the  units  as  repairable  or  n on-repairable . It  is  assumed  that  inspection 
takes  a negligible  amount  of  time  and  the  probability  of  failed 
units  being  repairable  is  known  and  is  the  same  for  all  arrivals. 

We  will  not  consider  the  decision  rules  for  classifying  the  items. 

After  inspection,  the  repairable  units  are  sent  to  repair  facilities 
where  repairs  are  performed  on  a f irst-come,  first-served  basis.  We 
s’,  a LI  consider  only  continuous  repair  process-,  that  is,  no  batching 
is  allowed  at  the  repair  facilities.  It  is  assumed  that  repaired  units 
bcha.o  exactly  like  new  ones  in  tneir  performance  characteristics. 

Upon  completion  of  repairs,  the  unit  immediately  joins  the  stock  of 
serviceable  units  at  the  supplv  point. 
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Supply  points  stock  the  ready -for-isce  units  to  resupply  the 
customers.  They  receive  their  supplies  from  two  sources:  the 

air  facilities  ind  th<  external  supplier.  It  is  i umed  that 
points  can  stock  an  unlimited  number  of  units. 

An  external  supplier  can  supply  an  unlimited  number  of  units  to 
the  system  within  a known  duration  of  time.  This  duration,  known  as 
procurement  lead  time,  is  assumed  to  be  the  same  for  all  orders 
independent  of  the  size. 

Thus  the  main  functions  of  a recoverable  item  Inventory  system, 
in  general,  are  to  fill  the  customer  demands,  to  diagnose  (inspect) 
the  failed  units,  to  repair  the  failed  but  recoverable  units,  and  to 
procure  units  from  an  external  supplier.  We  have  already  described  the 
inspection  and  repair  functions.  We  shall  consider  the  folxowing 
policies  for1  supply  and  procurement. 

The  Supply  Policy: 

Demands  are  satisfied  from  the  ready- for-issue  stock  (also 
referred  to  as  on-hand  inventory)  at  a supply  point.  Upon  arrival  of 
a requisition,  the  customer  is  immediately  shipped  the  quantity 
requested  if  there  is  sufficient  on-’nand  inventory.  If  there  are 
not  enough  units  on  hand,  then  all  the  units  in  stock  are  dispatched 
while  the  balance  is  backlogged.  In  either  case,  the  batch  of  failed 
units  is  sent  for  inspection.  The  backlogged  demands  are  satisfied  on 
a first-come,  first-served  basis  when  the  next  supply  arrives  from  an 
external  supplier  or  from  repair  facilities.  The  assumption  of 
backlogging  may  not  be  appropriate  in  some  cases,  but  is  representative 
of  situations  such  as  military  organizations  where  a captive  market 


situation  exists. 


The  Procurement  Policy: 


In  order  to  ake  up  for  system  losses  due  to  condemnations, 
new  units  are  procured  from  the  external  supplier.  We  shall  consider 
a continuous  review  reorder-point  (s),  reorder-level  (S)  pro- 
curement policy.  The  policy  is  based  on  the  inventory  position 
which  is  defined  as  the  sum  of  the  units  on  hand,  on  order  and  in 
repair  minus  backorders.  When  the  inventory  position  drops  below  the 
level  s,  a procurement  order  is  placed  so  as  to  bring  the  inventory 
position  to  the  level  S.  Thus  the  procurement  decisions  are  cognizant 
not  only  of  those  serviceable  units  on  hand  or  on  order,  but  also 
of  those  in  repair  or  awaiting  repair. 

The  main  reason  for  considering  the  continuous  review  (s,S) 
policy  is  to  deal  with  situations  in  which  the  presence  of  a computer- 
ized control  system  makes  it  possible  to  update  the  inventory  levels 
after  each  transaction.  Also,  this  policy  is  known  [7]  to  be  superior 
to  the  popular  continuous  review  reorder-point  Or),  fixed  order 
quantity  (Q)  procurement  policy  in  terms  of  total  reduced  costs  of 
procurement  and  carrying  inventory,  especially  when  the  order  size  is 
random . 

In  addition,  we  shall  consider  a continuous  review  (s,nQ) 
policy  because  of  its  mathematical  simplicity  [19].  In  this  policy, 

n.Q  units  are  ordered  each  time  the  inventory  position  drops  below  the 

level  s,  where  n is  the  largest  integer  such  that  the  subsequent 
inventory  position  is  between  s + 1 and  s + Q. 

To  provide  a better  understanding  of  two-echelon  systems,  we  shall 
first  study  a single  location  system.  An  outline  of  the  two  systems 
is  presented  in  the  following  subsections. 


I • 1 ■ 1 E ingle  Location  System 


Figxire  1.1:  Single  Location  System. 

Figure  1.1  shows  the  schematic  diagram  of  a single  location 
system.  System  demands  are  generated  at  one  single  point; 
procurements  at  the  supply  point  are  made  directly  from  the  external 
supplier. 

1.1.2  Two-Echelon  System. 

We  shall  investigate  a two-echelon  system  as  depicted  in  Figure 
1.2.  System  demands  are  generated  at  several  locations  called  bases, 
which  in  turn  receive  their  supplies  from  a central  location  called 
a depot.  The  depot  and  the  bases  are  also  called  the  upper  and  lower 
echelon  of  the  system,  respectively.  Each  location  in  addition  to 
being  an  item  stocking  point,  has  facilities  to  perform  repairs.  The 
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Figure  1.2  Two-Echelon  System 
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depot  stock  is  used  only  to  resupply  the  bases,  and  bases  are 

resupplied  only  by  the  depot;  that  is,  lateral  resupply  among  the 

bases  is  not  allowed.  Thus  procurement  of  units  from  the  external  supplier 

is  done  through  the  depot  only.  From  the  viewpoint  of  network  theory, 

the  system  looks  like  a parallel  activity  arborescence  structure. 

Upon  arrival  of  requisitions  at  a base,  the  failed  units  are 
sent  to  the  inspection  facility  whose  function  is  two  fold:  first, 

the  units  must  be  classified  as  either  recoverable  or  non-recoverable, 
and  then  if  repair  is  warranted  where  it  will  take  place  at  the 
base  or  at  the  depot.  The  latter  decision  depends  only  on  the  severity 
of  the  damage  caused  to  the  units  and  the  base  repair  capability. 

1.2.  The  Problem 

The  management  problem  in  both  the  single  and  two-echelon  systems 
is  to  establish  the  operating  rules  that  will  minimize  expected  cost 
per  unit  time  for  operating  the  system.  The  solution  to  the  problem 
is  usually  sought  in  the  environment  of  a limited  budget  and  a 
set  of  goals  to  reach  certain  levels  of  some  measures  of  system 
performance. 

The  major  cost  components  are:  cost  of  acquisition  per  unit  of 

the  item,  a fixed  procurement  setup  cost  independent  of  the  quantity 
ordered,  a fixed  backorder  cost  each  time  a stockout  occurs,  a time- 
weighted  cost  for  each  backorder,  a holding  cost  for  the  units  held  in 
stock,  and  a charge  for  the  units  held  at  the  repair  facilities.  Given 
the  repair  policy  and  the  number  of  repair  facilities,  specifying  the 
operating  rules  includes  determining  procurement  policy  parameters 
s and  S for  an  (s,S)  policy,  and  s and  Q for  an  (s,nQ)  policy. 
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Since  all  recoverable  failures  must  be  repaired  and  all  non- 
recoverable  failures  must  be  replaced  at  costs  independent  of  the 
procurement  policy  parameters,  we  can  ignore  both  repair  costs  and 
acquisition  costs  (other  than  fixed  procurement  or  set  up  costs)  as 
far  as  determination  of  these  parameters  is  concerned. 

Several  measures  of  effectiveness  have  been  used  for  inventory 
systems.  Following  Hadley  and  Whit in  [7],  some  of  these  are:  the 

probability  of  no  stock  on  hand,  the  expected  number  of  backorders  and 
expected  on- hand  inventory  at  any  time.  Feeny  and  Sherbrooke  [U] 
considered  fill  rate,  service  rate,  ready  rate  and  operational  rate  as 
measures  of  effectiveness  for  a base  stockage  system  with  no  condem- 
nations. For  a two-echelon  system  with  no  condemnation,  Sherbrooke 
[17 J found  the  expected  backorders  to  be  the  most  suitable  measure. 

The  Approach- 

In  order  to  solve  the  problems  as  described  above , several 
approaches  have  been  suggested  [3],  viz. expected  cost  analysis, 
stationary  process  analysis,  dynamic  programming  and  dynamic  process 
analysis.  We  shall  use  an  approach  based  on  stationary  process  analysis 
This  approach  is  more  appealing  to  us  since  it  is  applicable  to  general 
situations  and  it  is  also  computationally  less  complex.  The  approach 
uses  techniques  based  upon  Markov  processes  and  elements  of  renewal 
and  queueing  theory.  The  principal  problem  is  to  find  the  stationary 
probability  distributions  for  several  stochastic  processes.  These 
distributions,  if  they  exist,  are  functions  of  the  procurement  policy 

used  and  of  the  demand  distribution,  but  not  of  any  costs.  The  cost 
structure  tc  represent  the  objective  function  - expected  cost  per  time 
unit  - can  be  constructed  using  these  stationary  distributions.  Also, 
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the  various  measures  of  system  effectiveness  described  earlier  can  be 
obtained  from  these  distributions. 

In  order  to  obtain  expressions  for  the  objective  function  and 
system  performance  criteria,  we  must  find  the  stationary  distributions 
o^: 

1.  Inventory  position 

2 . On-hand  inventory 

3.  Number  of  backorders 

4.  In-repair  inventory. 

1.3 Scope  of  the  Study 

The  objective  of  this  study  is  to  obtain  an  exact  expression  for 
the  stationary  distributions  of  the  stochastic  processes  mentioned 
in  Section  1.2  for  the  single  location  and  two-echelon  inventory  systems 
described  in  Sections  1.1.1  and  1.1.2  , respectively. 

With  reference  to  the  dependent  demand  processes  mentioned  in 
Section  1.1,  we  shall  consider  the  following  two  inspection  policies. 

Batch  Model: 

Upon  arrival,  the  entire  batch  is  either  recoverable  or  non- 

recoverable.  Inspections  are  considered  to  be  repeated  independent  Bernoulli 

trials  with  probability  rD(say)  of  a batch  being  recoverable  and 

B 

probability  (1  - r ) of  it  being  non-recoverable . 

B 

From  a practical  viewpoint,  this  model  represents  situations  where 
the  units  of  a batch  fail  simultaneously  for  the  same  reason  and 
the  extent  of  damage  is  the  same  for  all  units  in  the  batch.  For 
instance,  the  maintenance  system  of  aircraft  engines  considered  by 
Muckstadt  [1 1]  in  his  MQD-METRIC  model  can  be  described  by  the  batch 
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model.  The  engine  consists  of  several  modules  which  contain  a large 
number  of  recoverable  units.  In  a failed  engine,  all  the  units  in  a 
module  (batch)  are  considered  to  have  sustained  same  extent  of 
damage  as  far  as  maintenance  (repairs)  is  concerned.  Mere,  the 
inspection  is  carried  out  to  decide  whether  the  module  will  be  repaired 
at  the  supply  point  (base)  or  will  be  sent  to  the  depot  for  repair. 

Unit  Model: 

Each  failed  unit  in  a batch  is  inspected  independently  to  deter- 
mine whether  it  will  be  repaired  or  condemned.  Inspections  are 
considered  to  be  repeated  independent  Bernoulli  trials  with  probability 
r,t  (say)  of  sending  a unit  to  the  repair  cycle  and  probability 
(1  - r ) of  condemning  it. 

From  a practical  viewpoint,  this  model  can  be  applied  to  the 
situations  where  units  failed  under  different  conditions  but  are 
submitted  in  a batch  for  replacement. 

For  the  single  location  system,  we  will  consider  both  the  cases 
of  independent  and  dependent  demand  processes  under  (s,S)  and 
(s,nQ)  procurement  policies.  A general  structure  will  be  provided 
to  obtain  the  stationary  distributions  for  the  cases  of  finite  and 
infinite  number  of  repair  facilities  and  general  repair  time 
distributions. 

The  two-echelon  sytem  will  be  analyzed  for  the  case  where  a 
Poisson  process  generates  the  demands  at  the  bases.  We  will  consider 
a one-for-one  (s-i,  s)  procurement  policy  at  the  bases  and  (s,S) 

(s,nQ)  policies  at  the  depot.  The  repair  times  at  the  depot  and  the 
bases  as  well  as  procurement  lead  times  between  a base  and  the  depot, 
and  the  external  supplier  will  be  assumed  to  be  deterministic  and  known. 
Also,  we  shall  assume  that  the  repair  capacity  at  all  locations  is  infinite. 
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For  both  systems  we  will  first  obtain  the  results  for  the  general 
case  where  an  item  will  either  be  repaired  or  condemned ^ results  will 
then  be  derived  for  the  special  cases  where  (i)  no  condemnations 
occur  and  (ii)  all  failed  units  are  condemned. 

1.4  Organization  of  the  Study 

We  begin  with  a brief  review  of  the  work  of  previous  authors 
in  Chapter  II. 

Chapter  III  is  devoted  to  the  study  of  the  single  location 
system.  We  identify  the  inventory  position  as  a semi-Markov  process 
and  obtain  its  stationary  distribution.  Following  this,  the  stationary 
distribution  of  the  stochastic  process  representing  both  the  backorders 
and  on-hand  inventory  is  obtained. 

Both  independent  and  dependent  demand  processes  under  (s,S) 
and  (s,nQ)  policies  are  examined.  The  case  of  a Poisson  process 
generating  the  failures  is  studied  in  depth.  At  the  end  of  the 
chapter,  some  long-run  averages  are  derived  from  the  stationary 
distributions. 

The  two-echelon  system  is  studied  in  chapters  IV  and  V.  In 
Chapter  IV,  the  case  of  unit  demand  at  the  bases  is  examined.  The 
case  of  an  arbitrary  order  size  distribution  at  the  bases 
is  considered  in  Chapter  V.  In  both  the  chapters,  the  special  cases 
of  no  condemnation  and  complete  non-recoverability  are  considered. 

In  Chapter  VI,  we  assess  the  degree  to  which  Sherbrooke's 
results  [17]  can  serve  as  approximations  to  our  exact  results.  Finally, 
Chapter  VII  contains  concluding  comments.  In  this  chapter,  we  also 
indicate  areas  for  future  research. 


CHAPTER  II 


REVIEW  OF  SOME  PREVIOUS  WORK 

Most  of  the  vast  literature  on  inventory  theory  is  addressed  to 
consumable  items  and  it  is  not  directly  applicable  to  inventory  systems 
of  recoverable  items.  We  shall  present  a brief  review  of  the  work  of 
previous  authors  for  single  location  and  multi-echelon  systems  for 
recoverable  items.  We  will  set  aside  the  numerous  simulation  models 
that  have  been  developed  and  applied  to  specific  situations.  Also,  we 
will  concentrate  on  vrork  dealing  uith  stochastic  models  for  these 
systems . 

2.1  Single  Location  System 

Considerable  attention  has  been  devoted  to  the  analysis  and  design 
of  this  problem  as  a whole  and  to  its  subproblems.  Feeney  and  Sherbrooke 
[5]  investigated  stochastic  recoverable  item  models  that  assume  compound 
Poisson  demand  distributions  and  complete  recoverability  of  failed  items. 
Schrady  [IS]  examined  a deterministic  model  that  permits  condemnations. 

In  a subsequent  survey  report  [16],  Schrady  described  approximate  solu- 
tions to  both  continuous  and  periodic  review  models,  and  continuous  as 
well  as  batch  repairs . 

Allen  and  D'Esopo  [1]  allowing  condemnations  obtained  approximate 
stationary  results  for  expected  number  of  backorders  and  other  measures 
of  effectiveness.  They  assumed  a Poisson  demand  distribution  and 
deterministic  (positive)  repair  and  procurement  lead  times.  In  a sub- 
sequent paper,  Simon  and  D'Esopo  [21]  obtained  exact  results  for  the 
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same  model  with  a relaxation  of  the  assumption  made  in  [1]  that  Poisson 
process  generate  the  failures.  They,  however,  assumed  that  recoverable 
and  non-recoverable  demand  processes  are  independent.  In  all  the  above 
references  repair  facility  was  treated  as  an  infinite  server  queue. 

Francois  Lureau[10]  viewed  the  problem  as  a queueing  process  and 
obtained  stationary  results  for  the  expected  number  of  backorders.  In 
addition,  he  obtained  the  stationary  distribution  for  waiting  time  of 
a customer  before  being  resupplied.  Assuming  complete  recoverability 
throughout,  he  developed  models  allowing  general  repair  time  distri- 
butions and  finite  as  well  as  infinite  number  of  the  repair  facilities. 
Recently5 Richards  [13]  examined  the  problem  with  condemnations.  His 
results  allow  for  random  lead  time  and  alternate  repair  disciplines. 

He,  however,  maintained  the  assumption  that  recoverable  and  non-recoverable 
demand  processes  are  independent. 

Our  analysis  is  more  general  in  that  it  permits  the  recoverable 
and  non-recoverable  demand  processes  to  be  dependent  as  well  as  inde- 
pendent. In  addition,  for  dependent  demand  processes  we  consider  batch 
and  unit  inspection  policies. 

2 .2  Two-Echelon  System 

A fundamental  work  on  the  two-echelon  system  was  the  development 
of  METRIC  (Multi-Echelon  Technique  for  Recoverable  Inventory  Control) 
by  Sherbrooke  [17]  for  a completely  conservative  system  that  does  not 
allow  item  condemnation.  He  considered  the  problem  of  allocating 
several  units  among  a depot  and  several  bases  in  order  to  minimize  the 
total  expected  number  of  backorders  at  bases  within  the  limitation  of 
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a given  budget.  He  assumed  a compound  Poisson  demand  distribution  at 
each  base  and  a cne-for-one  procurement  policy.  Depending  upon  the 
nature  of  failures , repairs  could  be  performed  at  the  base  where  the 
demand  originated  or  at  the  depot.  The  resulting  expressions  are 
approximate , and  a special  case  of  our  model  can  be  used  to  check  the 
accuracy  of  his  results.  Sherbrooke  assumed  an  arbitrary  repair  time 
distribution  though  his  results  depend  only  upon  the  means  of  these 
distributions.  We  shall  assume  deterministic  repair  times.  Sherbrooke 

also  presented  an  approximate  method  for  including  item  condemnation 

• 

m his  model,  assuming  that  procurements  are  made  on  a periodic  review 
basis.  For  Sherbrooke's  model  of  the  conservative  system,  Muckstadt 
Cl2]  developed  a computationally  more  efficient  approach  than  the 
previous  work  by  Sherbrooke,  for  determining  the  optimal  system  stock 
levels . 

A variation  of  the  METRIC  model  was  introduced  by  Simon  [20]  to 
obtain  the  exact  expressions  for  the  stationary  distributions  of  on-hand 
inventory  and  of  the  backorders  at  the  bases.  This  model  is  more  general 
than  METRIC  in  that  it  permits  non-recoverability  as  well  as  recoverability 
with  positive  condemnation  rates.  It  is  less  general  in  that  all  repair 
and  lead  times  are  deterministic  and  demand  distributions  at  the  bases 
are  simple  Poisson.  In  a subsequent  comment  on  Simon's  paper,  Kruse 
and  Kaplan  [9]  pointed  out  that  Simon's  derivations  were  valid  for  the 
two  special  cases  in  which  non  base-repairable  failures  are  either  all 
depot-repairable  or  all  non-depot  repairable.  For  Simon's  model,  they 
suggested  a simpler  method  of  deriving  the  probability  expressions  for 
the  number  of  backorders  at  the  bases.  We  use  the  same  approach  for  the 
case  of  arbitrary  demand  distributions  at  the  bases. 
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CHAPTER  III 


SINGLE  LOCATION  SYSTEM 


3 . 1 Introduction 

In  this  chapter  v:o  study  the  single  location  system  as  described 
in  Section  1.1.1.  The  analysis,  in  addition  to  offering  the  solution 
to  this  system,  is  also  applicable  for  the  upper  echelon  (depot)  of 
the  two-echelon  system  discussed  in  Section  1.1.2.  The  results  can  also 
be  used  to  find  an  approximate  solution  for  an  individual  location  in 
lower  echelons  of  a multi-echelon  system. 

The  stationary  distributions  of  inventory  position,  on-hand 
inventory , number  of  backorders  and  in-repair  inventory  will  be 
obtained  for  the  problem  described  in  Section  1.2.  For  t > 0,  let 
X(t)  = the  inventory  level  at  time  t which  consists  of  the 
units  ready  for  issue  minus  any  backorders, 

Q(t)  = in-repair  inventory  at  time  t , 

0(t)  = the  number  of  units  on  order  at  time  t from  the  external 
supplier, 

Dc(t)  = the  number  of  units  condemned  during  the  interval 
<0,t], 

D (t)  = the  number  of  recoverable  units  turned  in  during  the 
r 

interval  (0,t], 

C(t)  = the  number  of  units  repaired  during  the  interval  (0,t], 

0Q(t)  = the  number  of  units  ordered  during  the  interval  (0,t] 

and 

RQ(t)  = the  number  of  units  received  via  procurement  during 
the  interval  (0,t]. 
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Also,  for  any  stochastic  process  {P(t),t  0},  PCt^,!^)  = 

P(t2)  - P(tp. 

Obviously, 

(3.1)  Q(t)  = 0(0)  + Dr(t)  - C(t), 

X(t)  = X(0)  + C(t)  + RQ( t)  - D (t)  - D (t), 

r c 

and  0 (t ) = 0(0)  + 0Q(t)  - RQ(t). 

The  inventory  position  Z(t)  at  time  t is  defined  as 

(3.2)  Z(t)  = X(t)  + Q(t)  + 0(t). 

From  Eqs . (3.1)  and  (3.2) 

(3.3)  Z(t)  = Z(0)  + 0Q(t)  - D (t). 

We  shall  first  consider  the  continuous  review  (s,S)  procurement 
policy.  When  the  inventory  position  Z(t)  falls  to  the  level  s or 
below,  a procurement  order  is  placed  to  bring  Z(t)  to  the  level 
S (_>  s).  Upon  completion  of  a repair,  Q(t)  decreases  and  X(t) 
increases  by  the  same  amount;  thus  no  change  in  Z(t).  Similarly, 
upon  an  arrival  of  supply  from  the  external  supplier,  0(t)  decreases 
and  X(t)  increases  by  the  same  amount  and  therefore  Z(t)  remains 
unchanged.  Thus  Z(t)  may  change  only  at  the  epochs  of  customer  arrivals. 
Consequently,  procurement  orders  may  only  be  placed  at  demand  epochs. 

By  convention,  whenever  an  order  is  triggered,  Z(t)  is  meant  to 
include  the  demand  just  arrived  plus  the  order  that  the  demand  triggered 
so  that  Z(t)  = S at  such  epochs. 
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The  state  space  for  the  process  {Z(t),  t ^ 0}  is  the  finite 
set  E = {i|s+l  < i < S,  i integers}.  Z(t)  bounces  between 
s + 1 and  S during  each  ordering  cycle.  The  ordering  epochs 
R1,R2’"'’  are  regeneration  points  for  (Z(t),  t >_  0}  since  at  each 
such  epoch  the  inventory  position  restarts  at  S and  the  continuation 
of  the  process  thereon  is  a probabilistic  replica  of  the  previous 
cycle  (Figure  3.1). 


Z(t)  A 


Figure  3.1:  A sample  realization  of  the  inventory  position  with  the 

(s,S)  policy  and  random  order  sizes. 


*"  - • tc  two  possible  approaches  to  obtain  the  stationary 

distribution  of  Z ( ) . One  approach  is  based  on  renewal  theory 
and  * ho  other  uses  the  theory  of  semi-Markov  process.  Though  the 
underlying  at on  .s  remain  the  same,  we  sliall  use  the  second  approach 
in  our  analysis. 

V.'e  restrict  the  use  of  the  term  'limiting  distribution'  in  the 
context  of  the  stochastic  processes  with  discrete  index  parameter 
such  as  Markov  chains:,  whereas,  the  term  'stationary  distribution* 
will  be  used  in  the  context  of  continuous  index  parameter  stochastic 
processes  such  as  {Z(t),t  >_  0}. 

The  stochastic  process  {X(t),t  >_  0}  changes  its  state  at  demand 
arrival  epochs  (decreases),  at  completions  of  repairs  (increases), 
and  at  arrivals  of  procurement  orders  (increases).  The  state  space  of 
this  process  is  the  set  Z = {S,S-i,... ,1,0, -1,-2,...}.  The  positive 
values  of  X(t)  indicate  on-hand  inventory  while  the  negative  values 
indicate  the  existence  of  backorders  B(t)  at  time  t ; that  is, 

B(t)  - max(0,-X(t) ). 

In  order  to  obtain  the  stationary  distribution  of  X(t),  we  assume 
that  the  procurement  lead  time  is  a constant  t . 

For  t t let 

(3.4)  x.  .(t)  = Pr{X(t ) = jls(0)  = i}  i & E,  j e Z 

Anything  on  order  from  external  supplier  at  time  t - x will  have 
arrived  by  time  t and  anything  ordered  after  time  t - r will 
arrive  after  time  t.  Thus  the  number  of  units  received  via 
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procurement  during  (t-t,  tj  is  0(t  - t).  The  number  of  units 
arriving  from  the  repair  shop  during  (t-t,  t]  is  C(t-x,  t). 
Therefore , 


X(t)  = X(t  - t)  + 0(t  - t)  + C(t-T,  t)  - D (t-t,  t)  - I>c(t-i,  t). 
But  C ( t— t , t)  = Q(t  - t)  + D (t-t,  t)  - Q(t)'j  thus,  we  get 
X(t)  = Z(t  - x)  - Q(t)  - D(t-x,  t). 


Then  Eq.  (3.4)  can  be  written  as 


(3.5)  x._.(t)  = l Pr{X(t)  = j|z(t-x)=k;  Z(0)=i>  • Pr{Z(t-x)  = k|Z(0)  = i] 
J keE 


l Pr{Z(t  - t)  - Q(t)  - D (t-t,  t)  = j | Z(t  -x)  = k;  Z(0)  = i} 


keE 


• Pr{Z(t  - t)  = k|Z(0)  = i} 


kel 


= I I Pr{Q(t)  + D (t-x,  t)  = k-j jQ(t-x)  = m;  Z(t-x)  = k; 
(_  m=0 


• Pr{Q(t  - x)  ~ m|z(t  - x) 

• Pr{Z(t  - x)  = klz(0)  = i} 


Z(0)  = i} 
= k,  Z(0)  = i}J 


The  stationary  distribution  x(i)  = lira  x..(t)  is  obtained  by  taking  the 
limit  (as  t-*08)  of  the  right  hand  side  of  Eq.  (3.5).  The  stationary 
distributions  for  on-hand  inventory  and  the  backorders  can  easily  be  obtained 
knowing  x(j).  The  stationary  distribution  of  Q(t)  will  be  derived 
during  the  proe' ~s  of  obtaining  x(j). 
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In  the  next  section,  wo  assume  the  demand  processes  {Dc(t),t  > 0} 
and  {D  (t),t  >_  0}  are  independent.  In  Section  3.3,  these  demand 
processes  are  assumed  to  be  dependent.  In  Section  3.4  results  are 
obtained  for  the  case  where  the  demand  arrival  process  is  Poisson. 

In  Section  3.5,  an  (s,nQ)  procurement  policy  is  studied.  In  Section  3.6 
results  for  the  two  special  cases  of  complete  recoverability  and  complete 
non-recoverability  are  outlined.  In  Section  3.7  some  long-run  averages 
are  derived. 

3 . 2 Independent  Demand  Processes 
3.2.1  The  Model 

The  case  of  independent  demand  processes  arises  when  there  are  two 
independent  streams  of  failure  processes  responsible  for  the  recoverable 
and  the  non-recoverable  demands . 

The  recoverable  demands  arrive  at  the  supply  point  at  the  epochs 
0 12 

of  time  T = 0,  T , T ,...,  where  the  inter-arrival  times 
r r r 

Tr  ” ^r  ^ = ^->2,...)  are  independent  and  positive  random  variables 

with  common  distribution  function  A (t)  = Pr{Tn  - Tn  ^ < t}, 

r rr  — 

(t  >_  0;  n = 1,2,...).  Let  be  the  number  of  the  units  (repairable) 

"th  1 2 

turned  in  for  replacement  at  n epoch.  The  order  sizes  £ ,5^,...  are 
independent,  positive  and  integer- valued  random  variables  with  common 
probability  distribution  ^(j)  = PrU^  = j),  = 1,2,...;  n=l,2,...). 

These  random  variables  are  also  independent  of  the  arrival  process 

• It  is  assumed  that  the  inter-arrival  times  and  order  sizes  have 


finite  means  . that  is. 


21 


rQ  go 

a = / tA(dt)  < - and  dm  = I jtjj)  < ”• 
r o 1 r j=i  r 

Similarly,  the  non-recoverable  demands  arrive  at  the  epochs 

T = 0,  and  the  inter-arrival  times  Tn  - Tn  1 (n  = 1,2,...) 

are  independent  and  positive  random  variables  with  common  distribution 

function  A (t)  = Pr{Tn  - Tn_1  < t> , (t  > 0;  n = 1,2, . . . ).  The 
c c c — — 

• .12 

respective  number  of  the  units  £ »£»•••  demanded  are  independent, 

positive  and  integer-valued  random  variables  with  the  common  probability 

distribution  ^(j)  = Prtsjj  = j>,  (j  = 1,2,...;  n=l,2,...),  and 

these  are  also  independent  of  the  arrival  process  {T^}.  It  is 

assumed  that  <p  (1)  > 0,  and 
c 

CO  00 

a = / tA  (dt)  < » and  dm  = £ ” 

0 c j=l 

The  demand  processes  {D  (t),t  > 0}  and  {D  (t),t  > 0}  are 

r — c — 

independent  since  {Tn},  {Tn},  {£  n}  and  {£  n}  are  independent.  We 

r c r c 

also  assume  that  A and  A are  non-arithmetic  C231. 
r c 

For  any  n = 1,2,...,  let 


A^°(t)  = 0 for  t < 0, 

= / A (t  - y)A^n  l\dy)  for  t > 0, 

c 0 c ~ 


N (t)  = the  number  of  non-recoverabie  requisitions  arriving  during 
the  interval  (0,t], 

U(t)  = the  length  of  time  interval  between  time  t and  the  epoch 
of  the  first  non-recoverable  requisition  arriving  after  t 


v^(t)  = Pr{D  (t)  = k),  k = 0,1,..., 

K C 

v^(t,  t+t)  Pr{Dc(t,  t+x)  = k]  k = 0,1,...,  x ]>  0, 


and 


Then 


lim  w 


c 

k 


(t 


t+x) . 


k 

vf(t)  = 1 Fr{D  (t) 

K -n  c 

n-0 


k|N  (t)  = n}  • Pr{N  (t)  = n} 
' c c 


l 

n=0 


(n) 


(k) 


(A(n)(t) 

c 


t<n+1)ct». 


As  shown  in  the  reference  [23],  we  note  that  Vj_.(t)  satisfies  the 
following  integral  equations. 


(3.6) 


v^(t)  = 1 - Ac(t),  t >_  0 


and 

(3.7)  v°(t) 


k t 

I I <t>  (j)v^  -(t  ' x)A  (dx)’  k = • 

j=l  0 C C 


Since  A is  non-arithmetic,  we  have 
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, u 

(3.8)  lim  Pr{U(t)  < uj  = — J (1  - A (x))dx  (see  p.  97,  [23]). 

t-*»  ac  0 

For  each  t _>  0,  we  have 

(3>9)  w“(t,t  + t)  = Pr(U(t)  > t},  and 

k x 

(3.10)  w^(t,t  + t)  = £ J p (j)v£  *(t  - u)*dPr{U(t)  <_  u} 

k j=l  0 c 3 

k = 1,2,...  . 


From  Eqs.  (3.8  - 3.10)  we  get  the  stationary  distribution 


— / {1  - A ( i ) }du  k = 0 

|u  1 c 

, j c T 

(3.11)  w^  (t)  =( 

~~  l >vk-i <_r  " u),{1  ” Ac(u)}du 

c j=l  o J 

k = 1,2,... 


Similar  results  can  be  obtained  for  the  process  {D^(t),t  >_  0}. 

3.2.2  The  Stationary  Distribution  of  the  Process  {2(t),t  ^ 0} 

From  the  definition  of  the  inventory  position  (F.q.  (3.1)),  it 

is  clear  that  the  stochastic  process  (Z(t),t  >_  0}  changes  its 

1 2 

state  only  at  the  arrival  epochs  Tc,T*,...  . It  remains  unchanged 

1 2 

at  the  epochs  T ,T  ,...  because  in  the  case  of  recoverable  demands, 

r r 

X(t)  decreases  and  Q(t)  increases  by  the  same  amount  with  no 
change  in  Z(t).  Completions  of  repairs  or  arrivals  of  supplies  from  an 
external  supplier  do  not  result  in  a change  in  Z(t). 


It  is  clear  that 
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Z(t)  = Z(Tn)  for  Tn  < t < Tn+1;  n = 0,1,..., t > 0. 

c c — c ~~ 


'or cover  , we  have 


Pr{Z(Tn+-“-)  = j ; Tn+1  - Tn  < t|S(TU),  Z(TX),  . . . ,Z(Tn)  ; 
c c c — 1 c c c 


T°,  T1  , ...  , Tn} 
c c c 


Pr{Z(Tn+1)  = j;  Tn+1  - Tn  < tlzCln);  Tn}  almost  surely, 
c c c — 1 c c 


for  n = 0,1,..., t > 0;  and  i e £. 


Thus  Cz(t),  t > o}  is  a semi-Markov  process  whose  kernel  is  given  by 


Q( i ,3  t)  = Pr{Z(T^+1)  = j;  T^+1  - < t|z(T^)  = i}  i , j e E. 


Define 


P(i,j  ) = lim  Q(i,j  • t) 

-£->oo 

Then  P(i,j)  > 0 and  J P(i,j)  = 1,  so  that  the  P(i,j)  are  the 

jeE  n 

transition  probabilities  of  the  imbedded  Markov  chain  {Z(T^)}.  To 

derive  these  results  we  have 


Q(i,j  t)  = 0 

0(i , j t)  = Ac(t)4>c(i  - j) 

Q ( i > S ■ t)  = A (t)  Y 4>„(k) 
c , k c 
k>i-s 


s + l<_i<_j<_S-l', 
j <i,  s + 2 < i < s ; 
s + 1 < i < S. 


r 


Using  the  fact  that  lim  A (t)  = 1,  we  obtain 

t-x»  c 


P(i,j)  = 0 

P(i , j ) = <t>c(i~j  ) 

P(i,S)  = l 6 00 
k>i-s 


s + 1 < i < j <_  S - 1 , 
j <i-,  s + 2<^i^_S, 
s+1  < i < s- 


Nov;  we  have  the  following: 


Pn  ( i , j ) = Pr[Z(Tn)  = j|Z(TO  = i},  n = 0,1,. 

C W- 


1,  3 


e E. 


Let  Z = Z ( Tn ) , n = 0,1,...:  then  {Z  } is  a Markov  chain  imbedded  in  the 
n c n 

semi-Markov  process  {Z(t),t  > 0}.  The  transition  probability  matrix 

of  this  chain  is  P.  We  first  obtain  the  limiting  distribution  of 

the  chain  {Z  } and  then  the  stationary  distribution  of  the  process 
n 

{Z(t),t  > 0} 

Theorem  3.1.  The  limiting  distribution  v(j)  = lim  P (i,j)  f°r 

n-+°° 

i,j  e E,  of  the  imbedded  Markov  chain  {Z^}  exists  and  is  given 
by 


(3.12) 


...  m(S  -j ) 
v(j)  ~ l+M(S-s-l) 


s + 1 . ,S  - 1 , 


v(S) 


1 

i+mTs-s-1) 


where 


k 

M(k)  = l mOO, 

«.  = ! 


(3.13) 
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and  m(k)  satisfies 


(3.14)  m (1)  = (j>  (1) ; nuk)  = $ (k)  + / A (k  - q)m(q),  k = 2,3,. 

c c , 0 

Q 1 


Proof.  From  our  assumption  that  <(>^(1)  y it  follows  that  the 

chain  { Z ^ } is  irreducible.  From  the  theory  of  finite  Markov 

chains  [8],  we  know  that  in  a finite  irreducible  Markov  chain  all  states 

are  positive  recurrent.  We  consider  the  following  two  cases. 

(i)  0 < $ (1)  < 1. 

c 

In  this  case  the  chain  is  aperiodic.  V/e  know  that  for  an  irreducible, 
positive  recurrent  and  aperiodic  chain,  the  limiting  distribution 
Mj)  > 0}  exists  and  is  given  by  the  equations 


(3.15) 


v ( j ) = [ v(i)P(i,j), 

ieE 


j e E 


(3.16) 


l v(j)  = 1. 
j eE 


Using  Eq.  (3.14)  we  can  reduce  Eq.  (3.15)  successively  for 
j = S-l,  S-2,...,  s+1  to 


v(j ) = m ( S - j )v(S) . 


The  normalizing  condition  given  by  Eq.  (3.16)  and  the  Eq.  (3.13)  lead 

to  the  desired  results  (Eq.  (3.12)). 

(ii)  <}>  (1)  = 1 
c 


In  this  case  the  chain  is  periodic  with  period  S - s.  Also,  the 
matrix  P(i,j)  is  doubly  stochastic.  We  know  that  the  limiting 
distribution  for  this  chain  exists  and  is  given  by  the  Eqs.  (3.15)  and 
(3.16)  which  on  solving  yield 

(3-17)  v(j)  = — — for  all  ] £ E. 

b"S 

For  <t>cC  1)  = 1,  from  Eqs.  (3.14)  and  (3.13)  we  have  m(k)  = 1 for 
all  k e E,  and  M(S-s-l)  = S-s-1.  Substituting  the:  e into 
Eq.  (3.12)  we  get  the  relations  given  by  Eq.  (3.17). 

Q.E.D. 

We  now  find  the  stationary  distribution  of  the  process 

{ Z ( t ) , t > 0}. 

Theorem  3.2.  For  the  process  {Z(t),  t>0} , the  stationary  distribution 
n(j)  = lim  Pr{Z(t)  = j | Z(0)  = i}-,  i,  j eE  exists  and  is  given  by 

t-voc 

(3.18)  II(  j ) = v(  j ) j e E 

where  v(j)  are  given  by  Eq.  (3.12). 

Proof:  According  to  our  assumption,  the  probability  that  a transition 

will  take  place  within  an  amount  of  time  t,  given  that  process  has 
just  entered  state  i and  will  next  enter  j is  independent  of  i and 
j and  is  A^ ( t ) . Let  h(k)  be  the  expected  amount  of  time  spent  in 


a state  k during  each  visit , then 


mm 


2 b: 


= I P(k»j )a  . 

*■>  C 

jeS 


r a < 00  for  all  k e E. 
c 


For  the  case  0 < $ (1)  < lj  when 
c 

(H(j)  > 0}  exist  and  is  given  by 


the  chain 


is  aperiodic  we  know  that 


v(j  )~n(  j ) 

'Hi)  = V v(k)h(k) 
kr.E 


= v(j)  . 


For  the  case  tj> 


(1)  = 1,  when  the  chain  is  periodic  with  period 


S - s,  the  mean  recurrence  time  of  state  k is  given  by  (p.  90,  [14]), 


h(k,k)  = (S-s)  / tA  Cdt) 


0 


= (S-s)a 


for  all  k e E; 


ana 


Q.E.D. 


We  note  that  n(j)  does  not  depend  on  the  recoverable  demand  arrival 

process  { A^}  or  order  size  distribution  $r( • ) . 

3 .2.3.  The  Stationary  Distribution  of  t he_Procgss^X ( t ) , t _>_0j 

To  obtain  the  transient  distribution  of  the  process  X(t)  given 


by  Eq.  (3.5),  we  first  prove 


the  following  theorem. 


Theorem  3.3.  For  independant  demand  processes,  Q(t)  and  Z(t) 
are  independent  for  any  t > 0. 

F^oof : From  Eq.  (3.1)  we  have 


Q(t ) = Q( 0 ) + D (t)  - C(t)  and  Z(t)  = Z(0)  + 0Q(t)  - D (t) 
r*  ^ 


As  mentioned  in  Section  3.2.2,  recoverable  demands  do  not  affect  the 

inventory  position  and  thus  do  not  influence  the  procurement  orders. 

Therefore  0Q(t)  is  independent  of  D^(t).  Similarly,  C(t)  is 

independent  of  D (t).  Since  D (t)  and  D (t)  are  independent,  we 
c r c 

con  write 

Pr{Q(t ) = q(t) | Z(t)  = z(t)} 

= Pr{Q(0)  + Dr(t)  - C(t ) = q(t ) j Z(0 ) + 0Q(t)  - DC(t)  = z(t)} 

= Pr{Q(0)  + d (t)  - C(t)  = q(t )} 

r 

Thus  Q(t ) and  Z(t)  are  independent  . 

Q.L.D. 

Applying  Theorem  3.3  and  using  the  fact  that  D^(t  “ T,t), 

0(t)  and  0(t  - t)  + X(t  - t)  are  independent,  we  can  rewrite  the 

condit-ional  transient  distribution  of  X(t)  as  follows.  For 
t > x , we  have 


(3.19)  rr{ X( t)  = j|Z(t  - t)  = k,  Z(0)  = i) 

00 

= £ Pr{D  (t-t,  t)  + Q(t)  = k—  j | Q ( t - t)  = m;  Q(t  - t) 


m=0 


+ X(t  - t)  = k-m;  Z(0)  = i}  • Pr{0(t  - t)  = m|Z(0)  = i} 


CO 

= l Pr{D  (t-T,  t)  + Q(t)  = k- j | Q ( t - t)  = m;  Z(0)  = i} 
m=0  c 


‘ Pr{Q(t-t,  t)  = m|Z(0)  = i} 


= Pr{Dc(t-r,  t)  + Q(t)  = k-j|z(0)  = i} 


k-i 

= l PrtQ(t)  = k-j-d  |Z(0)  = i}  • Pr{D  (t-T,  t)  = d }, 
d =0  c c c 

c 


since  D (t-x,  t)  does  not  depend  on  Z(0). 

Now  following  the  theorem  of  convergence  of  independent  processes 
[2],  from  Eqs.  (3.5)  and  (3.19)  we  have 


(3.20)  x(j)  = lim  Pr{X( t)  = j|z(0)  = i}  for  j e Z 

t~*«o 

"k-j 

= l I lim  Pr{Q(t)  = k-j-d  | 2(0)  = i} 

keE  _d  =0  t 
c 


• dim  Pr{D  (t-t, 


t) 


• n(k) 


The  above  expression  describes  the  stationary  distribution  of  X(t). 
Eq.  (3.20)  can  be  evaluated  after  obtaining  the  stationary  distribution 
of  Q(t)  and  from  Eqs.  (3.11)  and  (3.18). 

The  stationary  distribution  of  Q(t)  is  obtainable  using 
queueing  theoretic  methods.  As  indicated  in  [6],  explicit  analytical 
results  are  not  available  for  a GI/G/C  (general  inter-arrival  and  repair 
time  distributions  and  a finite  number  of  repair  facilities)  system 
with  bulk  (batch)  input.  We  shall  consider  the  following  special  case: 
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GI/K/cq  (general  arrival  process-,  constant  repair  time  = R,  infinite 
repair  facilities)  system  with  batch  input. 

For  any  t >_  R,  the  number  of  units  in  repair  at  time  t is 
equal  to  the  number  of  the  recoverable  units  that  arrived  in  interval 
(t  - R,t]  ; that  is  , 

Pr{Q(t ) = k - j - d [Z(0)  = i>  = Pr{ D (t  - R,t)  = k - j - d }. 

C T C 

* - . . . C,  . 

Following  the  steps  similar  to  those  used  m deriving  lim  w^(t,  t+r), 
we  have 

(3.21)  lim  Pr{Q(t)  = k'|z(0)  = i} 
t;-x» 

00 

/ (1  - A (u))du  k'  = 0 

r R 

= ( k'  R 

— I J <>  (j')  vf  .(R  - u).  (1  - A (u) }du  k'  = 1,2,... 
a . L , r J k-i  r 

r 1=1  0 J 

V 

k' 

, r ...  r An),,  , .,A(n),  . (n+l).,., 

where  v (t)  = ) <p  (k ’ ) { A (t)  - A (c)}. 

k'  ur  r r r 

n=0 

From  Eqs . (3.11),  (3.18)  and  (3.21),  we  can  obtain  x(j). 


3.3  Dependent  Demand  Processes 
3.3.1  The  Model 

We  view  the  case  of  dependent  demand  processes  arising  from  a 
single  failure  process. 

The  demands  arrive  at  the  supply  point  at  the  epochs  of  time 
T°  = 0,t\t^,...  . The  inter-arrival  times  Tn  - Tn  ^ (n  = 1,2,...) 
are  independent  and  positive  random  variables  with  common  distribution 
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function  A(t)  = Pr{TU  - T^"1  < t},  (t  > 0-  i.  = 1,2,...).  The 
respective  number  of  units  d manded  are  independent,  positive  and 
integer-valued  random  variables  with  the  common  probability  distribution 
$(*)•  These  are  also  independent  of  the  arrival  process.  It  is 
assumed  that  4>(1)  > 0;  and 

CO  00 

a=/tA(t)<°°  and  dm=  [ j <j>(j)  < '*>  . 

0 j=l 


Here  the  demand  process  {D  (t),t  > 0}  and  (D  (t),t  > 0} 

j?  — c — 

are,  in  general,  dependent.  Similar  to  that  in  previous  section,  we  shall 
denote  the  n-fold  convolution  of  >{>(•)  and  A(  • ) by  <J>^n;(')  and 
respectively.  For  t >_  0,  let 
N(t)  = the  number  of  total  requisitions  arriving  in  (0,t], 

U(t)  = the  length  of  the  time  interval  between  t and  the 
epoch  of  the  next  arrival. 


\,k2(t)  - r^c(t>  ■ V Brlt)  * V Vk2  = 


wv  v (t.t  + t)  = Pr{D  (t,t  + t)  = k • D (t,t  + r)  = k„) 

y ^2  c X r*  z 

kl5k2  = 0,1,...  ,t  > 0 , 

wk  k (t)  = lirT1  wt-  v <t5t  + t)  . 


'l’"2  t-x»  kl’k2 


We  shall  derive  these  for  the  two  inspection  models  discussed  earlier. 
Batch  Model: 


In  batch  model,  where  the  entire  batch  of  the  failed  units  is 
either  repaired  or  condemned,  inspections  are  repeated  independent 
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Bernoulli  trials  vd.th  probability  r.,(0  < r < 1)  of  sending  a batch 
to  the  rupiiix’  cyclu  and  probability  (1  - r^)  of  condemning  the'  batch. 
This  divides  the  requisitions  into  recoverable  and  non-recoverable 
types.  For  t > 0,  let 

N^(t)  = the  number  of  non-recoverable  requisitions  arriving 
in  (0,t],  and 

N^(t)  = the  number  of  recoverable  requisitions  in  (0,t]. 

We  have  the  following. 


FrfhMt)  = k|N(t)  = n}'  = (^)  rR  (1  - rR)n  k 


'k ' ‘B  ‘B 

k = 0,1 , . . . ,n. 


kl+k2|  f 

(3.22)  v . (t)  = l 
kl’  2 n=0 


y0  (k)rB(1  - rB)n_k  ^(k)(k2)r(n-k)(k1) 


,[A(n)(t)  - A(n+1)(t)] 


k , k2  = 0 -.1  : • • • • 


Following  an  approach  similar  to  that  used  in  deriving  Eqs.  (3.9) 
and  (3.10),  we  get 


(3.23)  wk  k (t,t  + T)  = Pr{U(t)  > t},  = 0,  and 


2 T 


(3*24)  w k(t,t+T)  = r • l J 4>(j)v.  k . 

kljk2  B j = l 0 W3 


.(t  - u)dFr{U(t)  < u} 


1 T 


+ (1  - r )•  I / <f(j)v  . , (t  - u)dPr(U(t ) < u} 
B j = l 0 kl"3,K2 


kx,k2  = 1,2,...  . 


lf—V.  . *-  • - *K  A -W 


J 
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The  stationary  distribution  is  given  by 


(3-25) v (t  ) 

krk  2 


- f { 1 - A( u ) }du  k ,k  =0 

cl  J. 

T 

_ / , k2  T 

• “ l I k • (t  - u){l  - A ( u ) } du 

B d j=l  0 kl,K2  ] 

1 ^ * 

+ (1  - r )•-  / / <f> ( j > v . , (t  - u){l  - A(u)}du 

C E a j = l 0 V3,k2 

k,  ,k  -1 , 2 , . . . 

j.  /. 


Unit  Model : In  the  unit  model,  each  unit  in  a batch  is  inspected  to 

determine  if  it  will  be  repaired  or  be  scrapped.  The  inspections  are 
repeated  Bernoulli  trials  with  probability  r (0  r^  < 1)  of 
repairing  a unit  and  probability  (1  - r^)  of  scrapping  the  unit. 
This  divides  the  units  into  recoverable  and  nor.-recoverable  classes. 


We  can  write 


(3.2G)  v,  (t) 
kl’k2 


Vk2  kl  k2 

( )(1  ~ Fy)  ^ 


kl+k2 

y 

L 

n=0 


i(n)(kl+k2)(A!(t)-A(n+1)(t)) 


kx,k2  = 0,1,2,... 


Also, 


(3.27)  w (t,t  + t)  = Pr(u  (t ) > t} 
kl  ’ k2 


krk2  = 0 


kl  k2  f j,+j„  11  j2  ] 

= l . I \(  j - ru)  (ru)  ‘ l *<h  + 12h 

j^l  j =1  K.  J1 


(3.28)  w , (t,t  + t)  = 
kllk2 


0 


v,  . . (t  - u)dPr{U(t) 
kf-1  ik^-] 2 


,k^  “ 1 *>2  9 • • • 
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The  stationary  distribution  is  given  by 


r 

00 

— J {1  - A(u) }du 

a t 


(3.29)  w*  (t)  = 

Vk2 


C?!(h*h 


{ ‘ V1  h 


k^^  = 0 


n ]2 

(1  - r ) r 
u u 


/ <Hj1+j0)v  . . (t-u){1  - A(u) }du ' 

0 1Ki  31’  X2  J2 

= 1,2,...  . 


3.3.2  The  St at ion ary  D i s t ribufjon  of  the  Process  (Z(t),t  >_  0} 

When  a single  process  generates  failures,  {Z(t),t  0}  changes 

1 2 

its  state  possibly  only  at  arrival  epochs  T ,T  ,...  . Completions 
of  repairs  or  arrivals  of  procurement  orders  do  not  change  Z(t). 
Following  the  arguments  similar  to  those  presented  in  Section  3.2.2, 
we  see  that  {Z(t),t  ^ 0}  is  a semi-Harkov  process  whose  kernel  is 
given  by 


Q ( i , j ' t)  = Pr{Z(Tn+1)  = j Tn+1  - Tn  1 t|Z(Tn)  = i} 

i,j  e E. 

In  the  following  two  subsections  we  obtain  the  kernel  and  the 
stationary  distribution  of  (Z(t).t  > 0}  for  the  two  models. 

3. 3. 2.1  Batch  Model 

Upon  a demand  arrival,  the  following  two  outcomes  are  possible: 
(i)  the  ba^ch  is  condemned  (prob.  = 1 - r0),  X(t)  decreases 
and  Z(t)  cl.  nges,  and 
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(ii) 


the  batch  xj  found  repairable  (prob.  = r ),X(t)  decreases 
and  Q(t)  increases  by  the  same  amount,  and,  therefore, 
Z(t)  remain  unchanged. 


Then , 


Q(  i , j t ) = 0 s + l£i<j£S-l 

Q(i,j  t)  = A(t)rB  s -f  - 1 £ j = i<S. 

Q(i,j-  t)  = A(t ) ( 1 - r )(|){ i - j)  j < i,s  + 2 < i < S 

D — — 

Q(i,S  t)  = A(t)(l  - r ) £ 4>  ( k ) s + l<i<^S-l 

k^i-s 

Q(S,S  t)  = A(t){rB  + (1  - rB)  £ <p( k)). 

k>S-s 


Using  the  fact  that  lim  A(t)  =1,  we  obtain  P(i,j)  = lim  Q(i,j;  t). 
The  transition  probability  matrix  P(i,j)  of  the  imbedded  Harkov 
chain  is  shown  on  next  page.  We  have  the  following  theorem. 

Theorem  3.4.  The  stationary  distribution  Il(j)  = lim  (Z(t)  = j|Z(0)  = i}, 

t-x» 

i,j  e E exists  and  is  given  by  Eq.  (3.13). 

Proof.  From  the  matrix  P(i,j)  it  is  clear  that  the  chain  is 

irreducible  since  <j>(l)  > 0 and  0 < r < 1.  Therefore  all  the  states 

are  positive  recurrent.  For  0 < <f>(l)  < 1,  the  chain  is  aperiodic 

and  limiting  distribution  is  given  by  Eos.  (3.15)  and  (3.16)  which  on 

solving  yield  Eq.  (3.12).  For  <j>(l)  = 1,  the  chain  is  aperiodic  for 

0 < r < 1 and  is  periodic  with  period  S-s  for  r =0. 

3 B 

Therefore  the  limiting  distribution  is  given  by  Eq.  (3.17). 

Proceeding  as  in  the  proof  of  Theorem  3.2.  we  obtain  the  desired 
result . 

Q.E.D. 


J 
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3_.  3_12_._2__Unit_  Model 
Loti 


p(k)  = Pr{the  number  of  non-recoverable  units  in  a batch  = k} 


s £ Pr{the  number  of  non-recoverable  units  in  a batch 
d=k  = k|hatch  size  = ajPrioatch  size  = d>. 


(3.30) 


= l 


d=k 


(d)  (d-k)(  _ 

Kk  U 'J 


4>(d). 


Then  upon  the  arrival  of  a demand  , Z(t)  changes  v;ith  probability 

oo  00 

V p(k)  and  remains  unchanged  with  probability  p(0)  = 1 - £ p(k), 

k=l  k=l 

Moreover,  for  0 < r^  <1  0 < p(0)  < 1.  We  have 


Q(i,j  t)  = 0 

Q(i,  j ; t)  = A(t )p( 0 ) 

Q(i , j t)  = A(t )p(i  - j ) 

Q(i,S  t)  = A(t ) l p(k) 
k>i-s 


s + l<_i<j^S- 
s + l<^j  = i<S- 
j<i,s+2<i<S' 
s + l<i<S  — 1- 


and 


Q(S,S;  t)  = A(t)  • (p(0)  + l P ( k ) } . 

k>S-s 


We  can  obtain  the  transition  probability  matrix  of  the 

imbedded  Markov  chain  Z = Z(T*  ) ,n  = 0,1,...  . (noting  that  lim  A(t)  = l) 

n t-x» 


Theorem  3.5.  The  stationary  distribution  TI ( j ) 
i,j  e F.  exists  and  is  given  by 


lim  Pr{Z(t)  = j |Z(0) 
■£-►00 


i} 


5.  2iat 


3 j 


(3.31) 


n(j)  = T 


S(S-j) 

+G(S-s-l) 


, j = s + 1,. . . ,S  - 1; 


n (s)  = 


l+G(S-r.-l) 


where 

(3.32) 


G(k)  = l gU) 
1=1 


and  g(k)  being  determined  recursively  from  the  equations 


(3.33) 


g(l)  = p(l)/(l  - p( 0) ) 
k-1 

gOO  = Cp(k)  + J p(q)g(k  - q)]/[l  - p(0)] 
q-1 

k = 2,3,...  . 


Proof:  Suppose  0 < ry  < 1 then  0 < p(l)  < 1.  Following  arguments 
similar  to  those  in  the  proof  of  Theorem  3.1,  we  note  that  the  chain 
is  irreducible  and  all  the  states  are  positive  recurrent.  Therefore 
the  limiting  distribution  (v(j)  > 0}  of  the  chain  exists  and  is 
given  by  Eas.  (3.15)  and  (3.16).  Using  Eq.  (3.33)  we  reduce 
Eq.  (3.15)  successively  for  j = S - 1,S  - 2,...,s  +1  to 

v( j ) = g(S  - j )v(S) . 

The  normalizing  condition  (Eq.  (3.16))  and  the  Eq.  (3.32)  yield 


g(S  - j)  4 
1+G( S-s-i ) 


/ • \ r\\  ^ 


= s + 1,...,S  - 1 . 


v(S)  = 


I+G(S-s-l ) * 


■ 
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Now  applying  Theorem  3.2  we  obtain  the  desired  results. 

Q.E.D. 

The  special  cases  of  complete  recoverability  and  complete  non- 
recoverability  are  discussed  in  Section  3.6. 

3.3.3  The  Stationary  Distribution  of  the  Process  {X(t),t>.  0} 

We  consider  a constant  repair  time  R and  an  infinite  number 
of  repair  facilities.  It  is  assumed  that  R < r.  In  order  to 
obtain  the  transient  distribution  x.^(t)  t >_  t,  given  by  Eq.  (3.5), 

t j 

we  first  obtain  the  joint  probability  distribution  of  {Q(t),t  > 0}  and 

{Z(t),t  > 0}  . If  the  demand  processes  are  independent  then  the  two 

processes  are  independent  (Theorem  3.3).  This  offered  a considerable 

simplification  in  evaluating  x..(t)  in  Section  3.2.3. 

i] 

Since  only  the  repairable  units  received  during  (t-R,t]  are 
in  the  repair  cycle  at  time  t (t>_R),  we  can  write, 

(3.34)  Pr{Q(t)  = m;  Z(t)  = k)Z(O)  = i} 

= l Pr{Q(t)  = m Z(t)  = k|Z(t  - R)  = 2 • Z(0)  = i} 
zeE 

• Pr{Z(t  - R)  = z|z(0)  = i) 

= l Pr{D  (t  - R,t ) = m j Z(t)  = k|Z(t  - R)  = z;  Z(0)  = i} 
zeE  r 

• Pr{Z(t  - R)  = z | Z( 0 ) = i},  m = 0,1,...,  and  k c E. 

Let 

N^(t  - R,t)  = the  number  of  procurement  orders  placed  during 
(t  - R,t], 


Figure  3.2:  A sample  realisation  of  the  process  {Z(t),  t>0} 

during  the  interval  (t-R,  t]. 


The  random  variables  Y^,  Y^-Y^  (i=2,3,...)  are  mutually 

independent.  Furthermore,  we  have 

(3.35)  PrtYi  < y|z(t-R)  = z}  = £ Pr{Dc(t-R,  r-R+y)  = J^}; 

k. >z-s 

and 


Pr{Y.  - Y.  < y}  = £ Pr(D  (y)  = k },  i=2,3,...  . 

1 k^S-s 
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ihe  process  { YJ  is  a rerewal  process  when  z = S,  and  it  is  a so- 
cs+led  delayed  renewal  process  when  z ^ S (see  Ross  [14]).  For 


p 1,2, ; let 


&!  (t,y)  = Pr{  N°(t-R , t)  = p,  Y < y;  D (t-R,  t-R+Y  ) = 
x.,z  r ’ p — J r p 

We  have. 


(t,y) 


and 


Pr{Dc(t-R, 


t-R+y)  = k.^; 


D (t-R,  t- 
r 


(3.35) 


F(p) 

a,Z 


v Z 

<t,y)  -I  l 
0 k =0 

4- 


F(P-1) 

^-k2,S 


(t, 


y-u)F^  (t,du) 


For  0 y _<  R and  p=l,2 , . . . ; we  have, 

(3.37)  Pr{Z(t)  = k;  Q(t)  = m|l(°(t-R,  t)  = p;  Yp  = y; 

Z(t-R)  = z} 

= Pr {D^( t-R+y , t)  = S-k;  Dr( t-R+y,  t)  = m-JllN 

YP  = y;  Dr(t_R’  t-R+y>  = 

and  for  p = 0 we  have 

(3.38)  Fr{ Z( t)  = k;  Q(t)  = mj Z(t-R)  = z} 


Pr{D  (t-R,  t)  = z-k ; D (t-R,  t)  = m} 
c r 


fw  (t-R,t ) , 
I z-k,m 


■ A J Z(t-R)  = z}  . 
R+y)  = &} 
p 2,3,...  . 


(t-R,  t-R+y)  = i- ; 

°(t-R,  t)  = p; 
Z(t-R)  = z}. 


for  z >_  k, 
otherwise  . 


From  Eqs.  (3.36  - 3.38)  we  get 


(3.39)  Pr{D  (t-R,  t)  = m;  Z(t)  = k|z(t-R)  = z;  Z(0)  = i} 
r 


m R 


z-k,m 


(t-R,  t)  + I l ( 
p=l  &=0  y=0 


.Substituting  Eq.  (3.39)  into  Eq.  (3.34)  we  obtain  the  joint  probability 
distribution  of  Q(t)  and  Z(t),  given  by 

(3.40)  Pr(Q(t)  = nr,  Z(t)  = k|Z(0)  = i} 


m R 


V 

= L 
zeE 


. (t-R,  t)  + y y / wq  v m-o 

VZ"k>m  p=l  a=0  y=0  S-k’m  * 


'•  Pr{Z(t-R)  = z|Z(0)  = i}. 


We  can  now  obtain  the  expressions  for  the  two  probability  terms  in  Eq.  (3.5). 

Since  D (t-x,  t)  is  independent  of  Q(t-x)  and  Z(t-x), 
c 

and  Q(t)  depends  on  Dc(t-x,  t)  only  through  D (t-R,  t),  we  have 


(3.41)  Pr{Dc(t-x,  t)  + Q(t ) = k-j|Q(t-r)  = m;  Z(t-x)  = k-  2(0)  = i} 

, . d 
k-3  c 

= l l Pr{D  (t-x,  t-R)  = q;  D (t-R,  t)  = d -q;  D (t-R,  t) 
d =0  q=0  c c r 

c 

= k-^-dc} 


, . d 

k-3  c / 00  i / 

III  wa  k (t"T’  t-R))  Td 

d =0  q=0  'k  =0  q,k2  ' ' a 

c 2 


• « * - 


(3.42)  Pr{Q(t-T)  = m|Z(t-x)  = k;  Z(0)  = 1) 

Pr{Q(t-T)-m;  Z(t-i )=k |Z(0)  = i} 
Pr{Z(t-x)=k|Z(0)=i} 


The  numerator  of  Eq.  (3.42)  can  be  obtained  from  Eq.  (3.40). 
Substituting  Eqs.  (3.41)  and  (3.42)  into  Eq.  (3.5)  we  get 


w , . J (t-R,  t)  , I ; w (t-x-R,  t-x) 

| KV, 


°°  m R 

♦ l l I 

p=l  £=0  y=0 


dy)/ 

Pr{Z(t-x-R)  = z|z(0)  = i} 


The  stationary  distribution  x(j)  can  be  obtained  from  Eq.  (3.42) 
using  the  Laplace  transform  approach  as  suggested  by  Sivazlian  [22] 


3.4  The  Casa  of  Comoo aid  Poisson  Demands 


'.’e  now  consider  the  case  where  the  failure  processes  are 
Poisson  processes.  Both  the  cases  of  independent  and  dependent 
demand  processes  outlined  in  Sections  3.2  and  3.3  are  investigated. 

3.4.1  Independent  Demand  Processes 

Suppose  that  non -recoverable  and  recoverable  failures  are 

generated  by  two  independent  Poisson  processes  with  rates  X ( > 0 ) 

and  A (>  0),  respectively.  Then,  in  terms  of  the  model  described 

r -At  -X  t 

in  section  3.2.1,  A (t)  = 1 - e c and  A(t)=l-e  ,t>0; 

c r — 

a = 1/A  (<  “■)  and  a = 1/A  (<  “).  Let  <J>  (•)  ($  (1)  > 0)  and 

c c r r rc  c 

• ) be  the  order  size  distributions  of  non-recoverable  and  re- 
coverable demands,  respectively.  Then  {Dc(t),t  >_  0}  and  {Dr(t),t  >_  0} 
are  independent  compound  Poisson  processes  with  parameters  A and 

A with  compounding  distributions  $ (')  and  4>  (•),  respectively, 
r c r 

Consequently,  from  Eqs.  (3.7)  and  (3.8)  we  have 

-A  t 

k , v e C (A  t)n 

(3.44)  V°(t)  = l ^n)(k)  

n=0 

and 

u -A  x -A  u 

lim  Pr{U(t)  £u}=A  f e C dx  = 1 - e ■ c . 

c 0 

Thus  from  Eq.  (3.10),  w£(t,t  + t)  a v£(t).  Therefore 


(3.45)  w£  (t)  = v£(t)  for  x >_  0. 
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Since  A (t)  = 1 - e c , t>0,  and  $ (*)  satisfy  the  ccjri.it ion c 
c — c 

of  Theorems  3.1  and  3.2,  the  stationary  distribution  II Cj  ) , j £ E 

is  given  by  Eq.  (3.18).  To  obtain  the  stationary  distribution 

x(j),  ) e Z,  given  by  Eq.  (3.20),  we  first  obtain  the  stationary 

distribution  of  the  process  (Q(t),t  > 0}. 

Let  G denote  the  repair  time  distribution  of  a unit  and  k 

be  the  number  of  repair  facilities.  Then  using  the  nomenclature 

from  queueing  theory  as  given  in  Gross  and  Harris  [6],  the  repair 

l<>  ] 

system  is  equivalent  to  a M r /G/k  queueing  system,  and  the 
stationary  distribution  of  the  process  {Q(t),t  >_  0}  is  equivalent 

tc  {p  ) , the  stationary  distribution  of  the  number  of  customers 

in  the  queueing  system.  As  mentioned  by  Gross  and  Karris,  no 

analytical  results  for  {p  } are  available  for  the  systems  like 

tV  " ro  

M .'G/l,  H /M/k  and  M 1 /K/l.  The  results  are  availctle  m 

the  form  of  a generating  function  or  a Laplace  transform.  In  our 

[*r3  . ' 

following  discussions,  we  shall  consider  a M /GO’)/00  repair 
system,  where  T is  the  mean  (finite)  of  the  repair  time  distribution. 
It  is  assumed  that  repair  time  is  the  same  for  all  units  in  a batch 
and  Q(0)  = 0.  Feeney  and  Sherbrooke  [5]  have  shown  that  for  such 
a system  {p^}  is  a compound  Poisson  with  parameter1  A^T  and 
compounding  distribution  <j>  ( • ) ; that  is, 

-A  T 

a , , e r (A  T)n 

(3.46)  lira  Pr{Q(t)  = q}  = £ <j>r  (q) 

n=0 

q = 0,1,. . . . 


We  consider  certain  special  cases. 


r 
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(I) Geometric  Order  Size  Distributions : 

The  geometric  order  size  distribution  has  been  extensively 
considered  in  inventory  theory  in  the  context  of  random  order 
size,  (for  example,  see  Hadl.y  and  Whitin  [7]).  With  geometric 
order  sizes,  the  resulting  demand  processes  {D  (t),t  >_  0}  end 
vD  (t),t  >_  0}  have  stuttering  Poisson  distributions.  As  noted 
by  Sherbrooke  [18],  a stuttering  Poisson  process  offers  considerable 
analytical  simplification  and  is  a 'natural'  candidate  for 
describing  the  demand  processes  in  a number  of  practical  situations. 
Let 

<(>  (0)  = 0,  o (k)  = (1  - a )uk_1  , k = 1,2,...  . 
c c c c 

and  i (0)  = 0,  $ (k)  = (1  - ct  )aK-1,  k = 1,2,...  . 

-n  V TT 

where  0 < a <1  and  0 < a <1. 

— c — r 

From  Lqs.  (2.14)  and  (3.13)  we  have 
m(l)  = 1 - a 

c 

m(2)  = (1  - a )a  + (1  - a )(1  - a ) = 1 - a 
c c c c c 

m(3)  = (1  - a )a“  + (1-a  )a  (1-a  ) + (1-a  )(1  - a ) 
c c c c c c c 

= 1 - a . 

c 


It  cen  be  easily  shown  that 


4 


( 3.47) 


n(j) 


1-a 

c 

l+(S-s-l)(l-ac  ) 


j = s + l,s  + 2,. . . ,S  - 1 


and 


ri(s) 


1 

l+(S-s-l) (1-a  ) ' 

c 


The  distribution  given  by  Eq. 

distribution  and  is  given  by 

\ -X  T 

C 

e 


(3.48)  w. 


(t)  = ' 


-X  T 

, c 


l 

3 = 1 


(3.45)  is  a stuttering  Poisson 


= 0 


> 1. 


that 


L J 

For  a M r /G(T)/“  repair  system,  it  follows  from  Eq.  (3.46) 
lim  Pr(Q(t)  = q}  also  has  a stuttering  Poisson  distribution 


■t-KO 

given  by 


(3.49) 


lim  Pr{Q(t)  = k-j-d  |Q(0)  = 0} 


\ k-i-d 

c 

a 

r 

L 


-X  T k-j-dc 

* " l 

n=l 


n! 


n 


dc  = k - 3 
dc  < k - j. 


Substituting  Eqs.  (3.47  - 3.49)  into  Eq.  (3.20)  we  can  obtain  the 
stationary  distribution  x(j),  j e Z. 

(2)  Unit  Order  Size 

When  units  are  demanded  one  at  a time,  that  is,  when  4^(1)  “ 1 
and  <!>_(!)  = 1,  it  follows  from  Eqs.  (3.18)  and  (3.45)  that 


"<i>  - sr? 


j e E 


i‘S 


ar.d 


Wd  (T) 


-A  T c 

e C( A x) 
c 

d : ' 
c 


d > 0. 
c — 


For  a H/G(T)/oo  repair  system,  from  Eq.  (3.46)  it  follows  that 

-A  T , . , 

e r (A  T)k_J*'ic 

lim  Pr{Q(t)  = k-  j-  d}  = — ~ --r—  , k > j + d 

t-xo  c (k  - -j  - d ).  - c 


Substituting  the  above  results  into  Eq.  (3.20)  we  get 


-At  d -AT  ... 

. S k-j  e C (A  t)  c e r (A  T)k  1 c 

(3.50)  x(j)=^-  I l c r 

S-s  L 


k-s+1  d =0 


d : 


(k  - j - d ).* 
c 


= J-  l 

Q_o  L 


-(A  t +A  T)  , . 

S e C r (A  t+A  T)k_: 
c r 


k=s+l 


(k  - j): 


As  noted,  by  Richards  [13],  it  is  interesting  to  compare  the  results 
given  by  Eq.  (3.50)  with  the  results  given  by  Hadlev  and  Whitin 
([V],  pp.  183-187)  for  a continuous  review  consumable  item  inventory 
system  with  Poisson  demands  with  parameter  A and  a constant  procurement 
lead  time  x , where 


(3.51) 


x(  j ) = ~ 


s 

l 

k=s-hl 


-Ax.,  ,k-i 
e (Ax) 

(k-]Ti 


Thus  a simple  change  of  parameter  from  Ax  to  (A^x  + A^T)  makes 
Eqs.  (3.51)  and  (3.50)  identical. 

He  conclude  by  emphasizing  that  results  obtained  here  hold  for  the 
case  where  an  item  has  n possible  independent  failure  modes.  For  mode  i, 
the  failure  process  is  a Poisson  process  with  rate  A^  and  the  order  -size 
distribution  is  <j>^(').  Furthermore,  if  the  failure  mode  belongs  to  a 


set  P.,  the  units  are  recoverable;  but  if  the  node  belongs  to 
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the  complement  set  C,  the  units  are  non-recoverable.  for  this  cane, 

the  results  obtained  in  this  section  hold  with  A = £ A., 

r . „ 1 
ieR 

AC  = V V}  = ir  X Ai<fi(’)  and  = xi  X Vi(,)* 


r ieR 


3.4.2  Dependent  Demand  f r ocesses 

Suppose  the  failures  are  generated  by  a Poisson  process  with  rate 

-At 

X ( >0 ) ; that  is,  A(t)  = 1 - e , t >_  0.  Also,  let  $(•)  denote  the 
order  size  distribution.  We  examine  this  case  for  batch  and  unit 
inspection  models. 

Batch  Model : 

Let  r be  the  probability  that  a batch  of  failed  unit  is  recoverable, 
B 

then  from  Eq.  (3.22)  we  have 


(3.52)  Pr{B  (t)  = k D (t)  = k.} 

1 J.  C c. 


kl+k2 


(!Vk(l  - r.)n-yk)(k,H(n-k)(k.) 


n=0  k=max(0,n-k^) 


'k  B 


k2  e A*at(Xrpt)k 

Y B 

L 1. 1 

k-C  k‘ 


-X(l-r_)t 


1 e 

4<k)(k.>  l - 

^ n-k=0 


(X(l-rB)t)‘ 


(n-k)I 


(n-k) 

<t>  (k  ) 
1 


From  (3.52)  it  follows  that  the  demand  processes  {0^(0,  t 0}  and 

{D  (t),  t > 0}  are  independent  compound  Poisson  processes  with  par'ameterc 
c — 

r X and  (1  - r )x,  respectively,  and  have  a common  compounding 
B B 

distribution,  $(•)•  Therefore,  the  results  of  Section  3.4.1  hold  for 
the  batch  model  with  X^ 

*(•). 


r_A , X = (1  - r )X  and  d>  C " ) = $ (*)  = 
B e B r c 
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Unit  Model: 

Let  rn  be  the  probability  that  a unit  in  a batch  is  recoverable. 
Then  from  Eqs.  (3.26)  and  (3.26)  the  demand  distribution  is  given  by 


(3.53)  v , (t)  = Pr{P  (t)  = k, . D (t)  = k0} 

x r * 


kl+k2 


\ K 


. , I k +k 

KJ  k2  ^ ^ 


(n). 


-At 

e (At)r 


(1  - V ru  j l * ‘VV-Tfr 

l_n=0 


and 


(3.54)  k (t)  = vk  k (t)  for  kn  ,k2  = 0,1, . . . 


The  stationary  distributions  H(j),j  e E and  x(j),j  £ Z can  be 

obtained  from  Eqs.  (3.31)  and  (3.43),  respectively.  The  derivation  of 

x(j),  however,  is  computationally  complex  for  a general  $■(•).  We 

consider  the  special  case  of  unit  order  size,  that  is,  $(1)  = 1. 

Using  the  fact  that  $^n\k^  + k^)  = 1 for  n = k^  + and 
^(n)^  +]<.)  = 0 for  n i k^  + 1<2»  it  follows  from  Eq.  (3.53)  that 


(3.55)  Pr(D  (t)  = k ; D (t)  = k.} 

c , . 1 r , , _ k, +k„ 

1*M„  ,kl  k2  e-»(At)  1 2 
) (1  ' V rU  Oc^p: 


h 

\ 


-V1  k 

e (TyXt)  2 e 

k < 

2' 


■(1-r  )A  k 

U ( (1-Pjj ) At ) X 

k.’ 


Thus  Dc(t)  and  P (t)  are  independent  Poisson  processes  with 
parameters  (l  - r^,)A  and  i^jA,  respectively.  Therefore,  the  results 
for  the  case  of  unit  order  size  obtained  in  Section  3.4.1  hold  in  this 
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situation  with  X = r,,A,  X = (1  - r,,)X,  and  a>  (1)  = <j>  (1)  = 1. 

r U c U r c 

In  addition,  Eqs.  (.3.52)  and  (3.55)  lead  to  the  obvious  conclusion 
that  for  the  unit  order  size  distribution,  the  batch  and  unit  models 
are  identical. 

3.5  The  Uniform  Distribution  of  Inventory  Position  and  the  (s,nO) 

Frocuromont  Policy . 

From  the  previous  sections,  we  make  the  following  observations 
on  the  stationary  distribution  of  inventory  position  under  an  (s,S) 
procurement  policy.  Tor  the  case  where  the  demand  processes  at'e 
independent,  it  is  clear  from  Fh  .or  :i  3.2  that  n(j),  j e E,  is 
uniform  over  {s  + 1,...,S}  for  <{>c(l)  = 1»  and  is  independent  of 

* ) . Similarly,  when  the  demand  processes  are  dependent  and  <^(1)  = 1, 
Jheorems  3.4  and  3.5  imply  that  H(j)  is  uniform  over  (s  + i,...,S}. 
Furthermore,  we  can  derive  the  conditions  for  which  the  uniform  distribu- 
tion is  obtained.  This  is  done  in  the  following  lemma. 

Lemma  3.1:  Under  an  (s,S)  policy,  Jl(j)  is  uniform  over  (s  + 1,...,S) 

if  and  cnly  if 

(a)  <f>  (1)  = 1 under  the  conditions  of  Section  3.2, 

c 

and  (b)  £(1)  = 1 under  the  conditions  of  Section  3.3. 

Proof:  We  have  already  discussed  the  ’if'  part.  For  the  'only  if 

part,  assume  H(j),  is  uniform,  that  is, 


(3.56) 


"«>  * STJ 


for  j e E. 


(a)  From  Lqs.  (3.13)/ and  (3.56)  it  follows  that  m(s-j)  = 1, 

s+l<j<S-l;  and  M(S-s-l)  = 3-s-i.  Substituting  into  Eq.  (3.14) 

we  get  <j>  (1)  = 1. 
c 

(b)  From  Theorem  3.4,  it  follows  that  above  proof  also  holds 
for  the  batch  model;  that  is,  (J>(1)  = 1-  for  the  unit  model, 

Eq.  (3.56)  implies  that  g(S-j)  = 1 for  s + l<^j<_S-l  and 
G(S-s-l)  = l . Lq.  (3.33)  then  implies  that  p(l)  = 1 - p(0)  and 
pOO  = 0 for  k = 2,3,...  . Substituting  this  into  Eq.  (3.30). 
we  get  <Kl)  = 1. 

Q.E.D. 

Because  of  its  mathematical  simplicity,  the  uniform  distribution 
of  inventory  position  has  been  extensively  considered  in  inventory 
theory.  For  example  see  Hadley  and  Whitin  ([7],  pp.  181-183),  Simon 
[19],  and  Sivazlian  [22],  But  Lemma  3.1  suggests  that  under  an  (s,S) 
policy,  II(j)  cannot  be  uniform  unless  (a)  <(>c(l)  = 1 or  (b)  $(1) 

However,  the  uniform  stationary  distribution  is  recaptured  under 
an  (s,nQ)  policy  for  arbitrary  order  size  distributions.  We  first 
describe  this  policy. 

The  (s,nQ)  Policy : 

Under  a continuous  review  policy,  when  the  inventory  position 
falls  to  the  level  s or  below,  n Q units  are  ordered  where  n is 
the  largest  integer  such  that  the  subsequent  inventory  position  is 
between  s + 1 and  s + Q.  Figure  3.3  shows  a sample  realization  of 
inventory  position  under  this  policy.  The  state  space  of  the  process 
{Z(t ) ,t  >_  0}  is  the  set  E'  = (i|s+1  < i <_  s+Q,  i integers}. 
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Figure  3.3:  A sample  realization  of  the  inventory  position  with  the 

(s,nQ)  policy  and  random  order  sizes. 


Tor  a consumable  item  inventory  system,  it  was  shown  by  Simon  [9] 
that  the  stationary  distribution  of  inventory  position  (units  on  hand  + 
units  on  order  - number  of  backorders)  exists  and  is  uniform  over 
(s+1, . . . ,s+Q  }.  In  the  following  theorem,  we  extend  his  results  for 
our  version  of  a recoverable  item  inventory  system. 


Theorem  3.6.  Under  an  (s,nQ)  policy,  n ( j ) = — , j e £’  both  for 

^ 

independent  and  dependent  demand  processes  under  the  conditions  of 

Sections  3.2  and  3.3,  respectively. 

Froof . We  obtain  the  transition  probabilities  P(i,j),  i,  j e E' 
of  the  imbedded  Markov  chain  for  the  process  (Z(t),  t>0}  under  an 
(s,nQ)  policy. 
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(a)  Independent  demand  processes : Proceeding  as  ln  Section 

3.2.1  for  an  (s,S)  policy,  we  obtain 


p(i,j)  = I <t>_(nQ  - <:->>,  i 1 U 

H.--1 

CO 

and  P(i, j ) = [ <j>  (nQ  + (i-j)),  j < i- 

n=0  c 

From  the  assumption  that  4>  (1)  > 0,  the  Markov  chain  {Z  } is 

c ** 

irreducible  and  its  limiting  distribution  v(j),  j e E'  is  given 
by  Eq.  (3.12).  Solving  these  equations  we  get  v(j)  = By 
theorem  3.2,  II ( j ) = ^-  . 

(b)  Dependent  demand  processes:  In  this  case,  we  shall  consider 
both  batch  and  unit  models  of  inspection.  For  the  batch  model, 

proceeding  as  in  Section  3. 3. 2,1,  we  have 

00 

P(i»j ) = (l-rfi)  l <j>(nQ  + i - j ),  j > i; 
n=l 

00 

P(i,j)  = r + l > j = 

n=l 

00 

and  P(i,j ) = (1-r  ) £ 4>(nQ  + i - j)  , j < i. 

u n=0 

The  matrix  P here  is  easily  seen  to  be  doubly  stochastic  , and  the 
rest  of  the  proof  is  similar  to  that  in  case  (a)  above. 

Proceeding  as  in  Section  3. 3. 2. 2,  w=  obtain  the  following  for  the 


unit  model. 
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00 

= l p(nQ  + i - j) 

j > i 

n=l 

and 

00 

p(i jj ) = l p(nQ  + i - j) 

j < i. 

n=0 


From  the  assumption  that  <j>(l)  > 0,  it  is  clear  that  Eq.  (3.30) 

implies  0 < p(l)  < 1 for  0 < < 1.  The  rest  of  the  proof  is 

similar  to  case  (a).  For  r^  = 0,  p(k)  = <j)(k)  and  the  resulting 

situation  is  the  same  as  case  (a),  with  <f  (')  = and  <j>  (k)  = 0 

c r 

for  k > 1. 

Q.E.D. 

When  the  demand  processes  are  independent,  the  stationary  distri- 
bution of  the  process  (X(t),  t>0}  under  an  (s,nQ)  policy  is  the 
same  as  obtained  in  Section  3.2.3  for  an  (s,S)  policy.  Also,  from 
the  discussion  in  Section  3.4.2,  it  follows  that  this  stationary 
distribution  is  the  same  under  both  policies  when  a Poisson  process 
generates  the  failures  and  a batch  model  is  used  for  inspection. 

Similar  results  hold  fer  the  stationary  distribution  of  the  process 
(Q(t),  t>0}.  The  derivation  of  these  stationary  distributions  under 
an  (s,nQ)  policy  for  general  dependent  demand  processes  is  not 
included  in  this  study. 

We  emphasize  that  the  (s,nQ)  policy  has  received  an  appreciable 
acceptance  in  practice  [19],  because  in  addition  to  the  mathematical 
advantage  of  a uniform  distribution  of  inventory  position,  it  permits 
the  use  of  an  economic  lot  size  Q.  The  difference  between  an  (s,S) 
and  an  (s,nQ)  policy  is  that  Z(t),  in  an  (s,nQ)  policy  is  in  E’ 
immediately  after  placing  an  order,  whereas  in  an  (s,S)  policy  it  is 
always  S immediately  after  placing  an  order.  The  two  policies  are  the  same 


Jt  i 
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if  <{>^(.1)  = 1 and  $(1)  = 1 for  independent  and  dependent  dcx.tnd 
process,  respectively. 

3.6  Special  Cases. 

3,6.1  Complete  Recover  ability  (no  condennat.i cr.s )_. 

In  this  case,  all  failed  units  are  recoverable  with  probability 
one.  The  system  experiences  only  one  type  of  demand  (recoverable) 
and  there  is  no  distinction  between  independent  and  dependent  demand 
processes.  There  are  no  procurement  and  inspection  functions  in  the 
system  and  consequently,  the  external  supplier  and  inspection  facilities 
are  eliminated  from  the  list  of  entities  in  Figure  1.1.  The 
resulting  system,  also  referred  to  as  a conservative  system,  resupplies 
itself  from  the  repair  facilities. 

System  demands  can  be  considered  to  arise  from  a single  arrival 
process  {A^Ct),  t>0}  with  order  size  distribution  4>^(  * ) as 
described  in  Section  3.2.  The  inventory  position  Z(t)  remains 
constant  for  all  t >_  0.  Let  Z(t)=S,  t > 0.  The  problem  of 
specifying  the  system  operating  rules  reduces  to  finding  a 
value  of  S that  minimizes  the  total  expected  holding  and  backorder 
costs  per  time  unit.  Eq.  (3.2)  reduces  to 

X(t)  + Q(t)  = S , for  t >_  0. 

Clearly  0 <_  Q(t)  < S indicates  inventory  on  hand  while  Q(t)  S 
indicates  the  existence  of  backorders  at  time  t.  Thus,  for  a 
conservative  system,  we  need  to  study  only  the  process  (Q(t),  t>0} 
to  obtain  the  above  expected  cost. 
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The  stationary  distribution  of  the  process  { 0 (t ) . t>0}  for 
both  a general  arrival  process  and  for  a Poisson  arrival  process  can 
be  easily  obtained  following  the  methods  developed  in  Sections 
3.2.3  and  3.4.1,  respectively. 

3.6.2  Complete  Non-recoverability. 

This  is  the  classical  inven toi'-y  problem  of  a consumable  item. 

The  inspection-repair  loop  is  eliminated  from  the  system  shown  in 
Figure  1.1.  All  the  supplies  are  received  from  the  external  supplier. 

The  system  demands  can  be  considered  to  arise  from  a single  arrival 
pi'ocess  {A  (t),  with  ordei1  size  distribution  <j>(')  as 

described  in  Section  3.2.  In  this  case,  Eq.  (3.2)  reduces  to 

Z(t)  = X(t)  t 0(t),  t >_  0. 

It  can  bo  easily  seen  that  H(j),  the  stationary  distribution  of 
{Z(t),  t>0  } , for  the  (s,S)  and  the  (s,nQ)  policies  are  given  by 
Theorems  3.2  end  3.6,  x’vspectively.  For  obtaining  x(j),  j e Z,  the 
stationary  distribution  of  (X(t),  t>0} , Eq.  (3.5)  for  an  (s,S) 
policy  can  be  simplified  to 

x . . ( t ) * l Pr{X(t)  = j | Z(t-r ) = k;  Z(0)  = 1}  • Pr{Z(t-t)  = k|Z(0)  * ij 
11  kr.E 

= l Pr{D  (t-T , t)  = k-j | Z(t-T ) * k;  Z(0)  = i] 
keE 


• Pr(Z(t-t)  - k|z(0)  * i). 


59 


Taking  the  limit  as  t we  get 


t(j)  = l n(k)  w,"  (t), 
keE 


where  wf  . (t)  is  obtained  from  Eq.  (3.11).  The  resulting  expression 
k-u 

for  x(  j ) is  equivalent  to  that  given  by  Tijns  ([23],  pp.  100-101). 
Similarly,  it  can  be  easily  shown  that  for  an  (s,nQ)  policy 


*(:>*£  l wk_i(T)’ 

Q kcE'  k ] 


where  v'  . (t)  is  given  by  Eq.  (3.11). 


3.7  Certain  Long-run  Averages. 

Combining  the  results  of  previous  sections,  we  can  determine 

certain  long-run  averages  which  may  be  used  to  structure  an  o1  ~ ectrve 

function  of  total  expected  cost  and  to  express  constraints  on 

system  performance  measures  such  as  mentioned  in  Section  1.2. 

We  first  obtain  the  stationary  expected  ordering  cost  per  unit 

time.  Let  K.(t)  represent  the  expected  number  of  orders  placed  in 

(0,t].  Also,  r .furring  to  Sections  3.2  and  3.3  let 

3 = a and  n = dn  for  independent  demand  processes, 
c c 

and  3 = a and  r,  = d™  for  dependent  demand  processes. 

From  elementary  renewal  theory  [23],  wc  have 


K.(t) 


Furthermore,  following  the  analysis  given  ->y  Tijr.s  (ir22j,  pp.  113-1111) 
the  expected  size  of  the  rf 'a  (n  > 2)  order  is  equal  to  If 

the  ordering  cost  of  k units  is  Kfi(k)  + c • k where  f(0)  = 0, 
and  <5 ( k ) = 1,  for  k > 0,  then  the  stationary  expected  ordering 
cost  per  unit  time  is  equal  to 


cn  KH(S) 

3 + e 


Under  an  (s,S)  policy  the  stationary  expected  number  of  backorders 
0 

is  given  as  } jx(j).  The  stationary  expected  number  of  the  units 

j--CT  q 

held  at  the  supply  point  is  l jx(j ).  The  expected  fraction  of 

j=l  s 

time  the  system  is  out  of  stock  is  1 - > x(j). 

3 = 1 

Similar  results  can  bo  easily  obtained  under  an 


(s»nQ)  policy. 


CHAPTER  IV 


TWO-ECHELON  SYSTEM  - UNIT  ORDER  SIZE 


4.1  The  Model. 

Wc  consider  a two-echelon  system  as  described  in  section 
1.1.2  with  J bases;  the  bases  are  numbered  from  1 to  J,  and 
the  depot  is  indexed  as  0.  The  failures  which  generate  the  system 
demands  occur  in  a Poisson  manner  with  known  parameter  X..  at 
base  j ( j = l,2 , . . . , J) . Upon  such  failures,  one  unit  of  the  item  is 
demanded  for  replacement.  A failed  unit  turned  in  at  base  j is 
repaired  at  the  base  with  probability  r..  and  is  shipped  to  the 
depot  for  repair  with  probability  (l-r.!p.  Thus  with  probability 
(l-r^Ml-p),  the  unit  is  condemned. 

We  further  specify  the  following  assumptions  used  in  the  medal. 

1.  The  bases  use  an  (s-1,  s)  policy  for  procurement  of  units  from 
the  der  >t.  The  depot  procurement  policy  is  a general  (s,S)  policy. 

2.  Backlogged  demands  at  each  location  are  supplied  on  a first -come, 
first-served  basis. 

'< . There  are  an  infinite  number  of  repair  facilities  at  each  location. 
The  base  repair  times  R.  and  the  depot  repair  time  R,.  are 
deterministic  and  independent  of  the  arrival  process  and  the  number  of 
units  in  repair.  is  the  same  for  all  the  units  received  from  all 

the  bases. 

4.  The  time  to  ship  a depot  repairable  unit  to  the  depot  from  a 
base  is  assumed  to  be  negligible.  In  reality,  it  can  be  absorbed  in 
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Rq.TIig  procurement  lead  times  for  the  depot,  and  for  the 

bases  arc  deterministic. 

Furthermore  va  assume , as  does  Simon  [20],  that  Rq  <_  tq;  that 
is,  the  depot  repair  time  does  not  exceed  the  depot  procurement  lead 
time.  The  analogous  results  for  the  case  > tq  c^n  le  deri'/ed, 
although  this  case  is  less  realistic. 

We  use  the  following  nomenclature. 

Nomenclature 


D 


r . 
1 


P 


(s.-l 

1 


s. ) 
1 


R. 

1 


T . 

1 


T 


0 


D.  (t) 
1 

D®(t) 


= total  demand  rate  at  location  j (j=0,i,. . . ,J;  j=0  denotes 
the  depot). 

= the  probability  that  a unit  that  fails  at  base  j will 
be  repaired  at  base  j . 

= the  probability  that  a failed  unit  that  is  not  base 
repairable  will  be  depot  repairable,  p is  the  same  for 

all  bases. 

a-  the  procurement  policy  used  at  base  j ( j=l  ,2, . . . ,J) . 

■u  the  procurement  policy  used  at  the  depot. 

= the  deterministic  repair  time  at  location  j ( j=0,l, . . . , j) . 

= the  deterministic  delivery  time  from  the  depot  to  base 
j ( j= 1,2 , . . . , J) . 

= the  deterministic  procurement  lead  time  from  the  external 
supplier  to  the  depot. 

-■  the  number  of  units  demanded  at  location  j ( j=0,l, . . . ,J) 
during  the  interval  (0,t]. 

- the  number  of  units  which  were  declared  base  repairable 
at  base  j during  the  interval  (0,t]  (j=l,2,. . . ,J). 

- the  number  of  units  sent  to  the  depot  for  repair  from 
base  j during  the  interval  (0,t]  (j=l,2,...,J). 


rr{u. 

3 


%T  ( 
* • V 

Lowe 


B?(t) 

J 


D?(t) 


“o(t) 


Qj(t) 

0,(t) 

3 


U.(t) 

3 


(*)=k) 


PCnlr..] 

v case 


= the  number  of  units  condemned  at  base  j during  the 
interval  (0,tj  ( j=l,2, . . . ,J). 

- the  number  of  units  demanded  from  the  djpot  by  base  j 
during  the  interval  (0,t]  (j=l»2, . . . »J) • 

~ the  number  of  units  rec  ived  at  the  depot  for  repair 
during  the  interval  (0,t]. 

= the  number  of  units  demanded  from  the  depot  as  a result  of 
condemnations  at  the  bases  during  the  interval  (0,t]. 

= the  inventory  position  at  time  t at  location 
j (j  = 0,l,...,.I). 

= the  number  of  units  in  repair  at  time  t at  location 

j (5*0,1, ..-.J). 

- the  number  of  units  on  order  at  time  t at  location 
j (j*0,l,...,J). 

= the  number  of  backorders  at  time  t at  location 
j (j~0.1, . . . ,J).  Negative  backorders  denote  on-hand 
inventory . 

- total  number  of  units  on  order  plus  in  repair  at  time 
t at  location  j ( j=0,l, . . . ,J).  Thus,  U_.(t)  = 

CK(t)  + JMt). 

= {sQ+l,  Sq+2  , . . . ,S0 } , the  state  space  of  the  process 
{Z  (t),  t>0}. 

- lim  PrO'.(t)  = k),  ( j-0,1, . . . ,J) . 

t>°°  ^ 

^-m,  »n 

r - — n=0,l,...  (Poisson  distribution  with  mean  m). 
ni 

-■  N(t  ) - N(t*)  for  the  process  (N(t),  t>0}. 

letters  are  used  to  denote  a particular  realisation  of  a 


r 
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random  variable. 

We  note  the  following  implications  of  ovr  assumptions. 

(a)  Because  the  bases  follow  an  (s-1,  s)  policy  and  the  time  to 
place  an  order  from  a base  to  the  depot  and  the  time  to  ship  a depot 
repairable  unit  to  the  depot  are  negligible,  we  see  that  for  t 0 


■-C(t)  = I D?(t ) , and  D^(t)  = l D?(t) 

j=l  3 j = l 11 


that  is , 


co!t)  = Do 


(t)  t DX(t)  . I D?(t). 


1=1 


(b)  For  j=l,2,...,J;  the  demand  processes  ^D^Ct),  t^o), 

^D^(t),  t>0^  and  tB?(t).  t>o)  are  mutually  independent  Poisson 
J - 3 ~ 

B D 

processes  with  parameters  A.  = r.A.,  A.  = (l-r.)PA.  and 
1 3 ] ] 3 J 3 

A.-  r (l-r__  )(1-P)A^  , respectively.  The  process  tD?(t),  t>o)  is  a 

Poisson  process  with  oaramoter  = (l-r.)A.. 

3 : 3 

D ° 

(c)  The  depot  demand  processes  (D0(t),  t^o)  and  {Py(t),  t^p)  are 

D J ^ 

independent  Poisson  processes  with  parameters  A _ = l A . and 
J 0 -1=1  'J 

C V ° J -*■ 

= Z respectively.  The  depot  total  demand  process 

j=l  : 

^D.Ct),  t>o}  is  a Poisson  process  with  the  parameter  A^  = 

/,  1?  = f x°. 

(d)  Because  of  the  infinite  number  of  repair  facilities  and  constant 
repair  times,  the  units  in  repair  at  base  j at  time  t(>_  Iw ) will 
be  due  to  the  base  repairable  failures  occurring  only  in  (t-R. , t]; 

g 

that  is,  Q.(t)  = P.(t-P..,  t).  Thus  for  t > R. , Q.(t)  is  a Poisson 

J 3 3 -33 

■ft 

variable  with  mean  A.R. , j=l,2,...,J.  Similarly  for  t R , 

p 

Q.(t)  is  a Poisson  variable  with  mean  A„R„. 

0 u u 


blL 


•»  *.  • - -*  > : 
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(e)  Because  depot  demands  are  for  a single  unit  at  a time,  (s,S)  arc 
(s,nQ)  policies  for  the  depot  are  the  same  with  Q = S-s. 

In  view  of  the  problem  described  in  Section  1.2,  our  goal  is  to 
obtain  the  stationary  distributions  of  the  processes  {Z_.(t),  t>0}, 
(B_.(t),  t>0}  and  {Q^(t),  t>0},  for  j=0,l,...,J.  The  overall 
abjective  of  the  model  is  to  find  the  policy  values  which  minimise 
total  expected  base  backorders.  Depot  backorders  are  of  interest 
only  insofar  as  they  affect  the  base  backorders.  The  basic  approach 
for  determining  the  stationary  distributions  is  described  in  Section  4.2. 
In  Section  4.3,  the  distributions  are  obtained.  The  results  for  the 
cas«s  of  complete  recoverability  and  compare  non-recoverability  are 
derived  in  Section  4.4. 


4.2  The  Basic  Approach  Fm?  Stationary  Distributions 

In  this  section  we  describe  the  basic  approach  for  determining 
the  stationary  distributions  of  the  processes  lZA t),  t>0), 

(Bj(t),  t>0}  and  {Q^(t),  t>0}  (j=0,l,. . . ,J)  when  the  bases  use  an 
( s— 1 , s)  procurement  policy.  The  approach  will  also  be  applied  in 
the  next  chapter  where  the  situation  of  a random  order  sice  at  the 
bases  is  dealt. 


4.2.1  The  Depot 

To  obtain  the  stationary  distributions  of  the  processes 

{z  (t),  t>0},  { Bq ( t ) , t_>0}  and  {Q0(t),  t>0} , we  first  find  the 

D C 

distribution  of  the  processes  (B0(t),  t>0}  and  (D~(t),  t>0}.  The 
depot  can  be  viewed  as  a single  location  system  and  the  results  of 
Section  3.4.1  apply. 
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6.2.?  The  Bases 

Because  the  bases  follow  an  (s-1,  s)  policy,  the  inventory 
position  is  constant;  that  is,  Z_.(t)  = Sj  for  all  t > 0 and 
j=l,2, . . . ,J. 

To  find  the  stationary  distribution  B. (*),  we  first  obtain 

J 

Pr{B_.(t)  = b}  for  b e { -s;  ,-s..+l, . . . ,0,1, . . . } , and  then 
3 j J 

evaluate  lim  Fr{B .( t)=b}.  ■ Tin,  appi'oach  is  the  one  given  by  Kruse  and 

p-,io  J 

Kaplan  [9]. 

Referring  to  Figure  4.1,  the  oniy  units  that  can  arrive  at  base 
i from  the  depot  by  time  t arc  those  on  order  by  tins  1 0 • 


*-VTj'  t_VTj  t“Tj  t " 

Figure  4.1:  The  sequence  of  events  at  base  j. 

This  depends  on  the  total  assets  (ready-for-issue  units)  available 
at  the  depot  by  time  t , the  total  demand  at  the  depot  during  the 

O 

interval  ( ,t„l  and  the  sequence  of  arrivals  of  requisitions  at  the 
depot  from  the  bases  during  the  interval  (t^t^].  This  is  so 
because  the  units  on  order  through  time  will  have  arrived  at 

the  depot  from  the  external  supplier  by  time  , and  any  units  not 
on  order  by  :.ime  t^  will  not  arrive  by  time  t^ . The  total  assets 
available  at  the  depot  by  time  include  the  units  on  hand  minus 

any  backorders  at  time  t^,  the  units  on  order  at  time  t^,  the 
units  in  repair  at  time  t^  and  the  units  received  for  repair 
during  the  interval  (t^t^.  This  equals  ZqO^)  + do^tl,t2^‘  Now, 
the  following  two  mutually  exclusive  situations  are  possible: 
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CASE  A:  The  total  depot  demand  during  the  interval  (v^t^J  does 

not  exceed  the  total  assets  available  by  time  t,(. 

In  this  case. 


d0^1’^2^  + do^l’^2^  — VV  + <^0^1’^'2^ 


d0(W  < W* 


D C 

Thus  all  the  depot  demands  ^q^I’1^  = d^(t1,t2)  + dgtt^t,,)  are 
satisfied  by  time  t^.  Only  the  depot  demands  dr(t2,t„)  against 

c 

the  stock  of  z„(t, ) - d„(t.,t„)  units  available  by  time  t„  determine 
how  many  demands  could  possibly  remain  unsatisfied  by  time  t~. 


CASE  B:  The  total  depot  demand  during  the  interval  exceeds 

the  total  available  stock  by  time  t^. 

In  this  case. 


dS(tl’t2)  + d0(tl’t2)  > W + d0(tl’t2) 


d0(tl’t2)  > z0(tl) ‘ 


Thus  there  is  no  stock  available  at  the  depot  at  time  t0  to  satisfy 

the  demands  d^it^t^).  Also,  there  is  no  guarantee  that  all  the 
D 

dn(t  ,t  ) demands  will  be  satisfied  since  this  depends  on  the 

V 1 4 

sequence  of  arrivals  of  d^(t^,t2)  and  d^(t^,t0).  Hence,  the  total 

Q 

demand  d (t  ,t  ) dravm  against  the  amount  of  z (t  ) + d (t  ,t  ) 

U 1 J 

determines  hov;  many  demands  will  remain  unsatisfied  by  time  t . 

O 


As  a consequence  of  the  (s.-l,  o . ) policy  at  base  j, 

3 3 

z.(t)  = ; that  is  -Bj (t ) + U.(t)  = s.  for  all  t > 0.  Th  n for 


any  b e {-s^  , -s.,+1,. . . ,0,1,. . . } and  for  any  t > rQ  + r. , 

event  B.(t)  = b occurs  if  and  only  if  U.(t)  = s.  + b.  Thus 
i J J 


the 


(4.1) 


Pr(B.(t)  = b)  - Pr{U. (t)  = s.+b}. 
1 ] 3 


In  view  of  the  previous  discussion,  Eq.  (4.1)  can  be  rewritten 


as 


(4.2)  PrtUj (t)  - s.+b} 


l 


Z0(tl)eE0  d0(tl,t2'1  ° 


Pr{Uj(t)=sj+b|DE(t1,t2)=dQ(t1,t2); 


0 


lV°VV)!*tD0<tl’t2)'dO(VVi  WW1} 


Vs  investigate  Pr{U.(t)  - s^+b J D^(t1,t2)=dQ(t^,t2 ) , ZQ(t^)=z0(t1 ) } 
for  the  following  two  mutually  exclusive  events, 


(A)  Dd(trt2)  < z0(tl), 

and  (B)  D^t^t^  > z0(t1). 

Using  the  independence  of  Dg(t  jt^)  and  Zg(t^)  WG  can  eiiVrcSi:  Eq.  (4. 
corresponding  to  events  (A)  and  (B)  as  follows, 


(4.3)  Pr(U. (t)  = s.+b} 


z0(ti)e-0 


I [Pr{U. (t)=s.+b>.  + Pr{U.(t)=s  +b)  3 

3 1 A J 3 b 


For  a given  z.Ct,),  Eqs.  (4.4)  and  (4.5)  represent  Pr{B  (t)  = b}  for 
Q J 

cases  A and  B5  respectively. 

To  evaluate  Eq.  (4.4),  let 

u.  (t)  = uht)  + U?(t)  , 

5 1 ) 

where 

U^"(t)  = the  sum  of  the  units  in  repair  at  base  j at  time  t and 
1 

the  units  for  which  orders  were  placed  on  the  depot  by  base 
j during  the  interval  (t3»t], 

and 

U^(t)  = the  units  ordered  from  the  depot  by  base  j during  the 

interval  (t  .t^]  that  remain  unfilled  by  time  tg. 

1 2 

Because  the  arrival  process  is  Poisson,  IK(t)  and  U_.(t)  are 
independent.  As  mentioned  earlier,  all  the  demands  levied  on  the 


7U 


depot  fi'om  base  j during  the  interval  (t„  ,t]  will  remain  on 
order  at  time  t.  Also,  the  base  repairable  demands  occurring 


only  during  the  interval  (t-Rj,  tj  will  be  in  the  repair  cycle  at 

base  j at  time  t.  Thus,  U1(t)  = P7(t-R. , t)  h D7(t„,t)  + 

D ] J J J ° 

D^.(t3,t)  and  the  probability  distribution  of  U^vt)  can  be  easily 

2 

ootained.  The  probability  distribution  of  U((t)  requires 

considering  the  sequence  of  arrivals  of  requisitions  from  the 

bases  during  the  interval  (t„,t  ].  Suppose  the  probability 

distributions  of  LtI" ( t ) and  U7(t)  ar’e  obtained  then  Eq.  (4.4) 

J 'J 

1 2 . 1 

can  be  evaluated  by  oonvoluting  U,(t)  with  U. (t).  Since  U7(t) 

j 3 j 

. . C 

is  independent  of  Zrt(t  ) and  D„(t,,t_), 

u JL  0 J. 


zo(V 


(4.G)  Pr(U.(t)  = s,tb}A  = 
J j A 


I 

4o!  W=°  L 


.+b 

i 

I Pr{U.(t)  = s.+b-d) 
d=0 


• Pr{Uj (t)  = d|D^(trt2)  = d^(trt2),  ZQ(t1)  = z0(t;L)} 

• Pr{DQ(t15t2)  = dy(t],t2)}. 


To  evaluate  Eq.  (4.5),  wo  r te  that  all  the  bane  demands 
levied  on  the  depot  during  the  interval  (t0,t]  remain  unfilled 
by  time  t.  In  addition,  some  demands  from  base  j placed  during 
the  interval  (t^,t2]  may  remain  unfilled  by  time  t^.  Let 

E'^(t)  = the  sum  of  the  units  in  repair  at  base  j at  tin. a t 
and  the  units  for  which  the  orders  Were  placed  on  the 
depot  by  base  j during  the  interval  (t2,t], 
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and 

2 

Uj(t)  = the  units  ordered  from  the  depot  by  base  j during  the 

interval  that  are  unfilled  by  time  t^. 

• , 1 2 

Obviously  U^v.t)  = U.(t)  + U.(t).  Because  of  Poisson  arrival 
1 1 1 

process,  U..(t)  and  U". (t)  are  independent.  Here  U7(t)  = 

3 3 

B _ _ 

, t)  + D.(t_,t)  + D7(t_,t)  and  we  can  readily  determine 

the  probability  distribution  of  l)}(t).  The  probability  distribution 
2 

of  th(t)  involves  consideration  of  the  sequence  of  ari'ivals  of 

requisitions  at  the  depot  during  the  interval  (t^t^L  Eq.  (4.5) 

1 2 

can  be  obtained  through  the  convolution  of  U.(t)  and  U7(t). 

1 C 

Again,  L\(t)  is  independent  of  ZQ(t^)  and  * *i:iUS 


(4.7)  PHU.(t)  = s . +b)  _ = 
J J B 


d0'tl’t2^=Z0^tl'+1 


s . +b 

7 Pr{u;}(t)  = e.+b-d} 

d=0  J ] 


• Pr{U7(t ) = d|D^(t1,t2)  = dQ(t15t2>;  = Zq^)} 

' ’W'S'VV  51  ❖ vh”' 


Upon  substituting  Eqs.  (4.6),  (4.7)  and  the  results  for 
Pr{ZQ(t^)  = Zp(t,)}  into  Eq.  (4.3),  we  obtain  Pr{lh(t)  = s_.+b), 
The  probability  distribution  for  the  process  (Q.(t),  t>0) 
is  easily  obtained  as 


Pr{Q.  (t)  = q}  = Pr{D^(t-R^,  t)  = q}  , for  t >_  R. . 
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4.3  The  Stationary  Distributions^ 

4,3,1  At  Th  j Ds  pot 

n 

Since  (D0(t),  t>0)  and  t^_0)  are  independent 

D C 

Poisson  processes  with  parameters  AQ  and  AQ  , respectively 
it  follows  from  Section  3.4.1  that 


(4.8)  lim  Pr{ Z_(t)  = k}  = n (k)  = — - — , for  k e E • 

t*»  0 VS0  0 


(4.9)  in  wV.i  - v • l , 


for  ~Sq ,-Sq+1  9 • • . j 


and 


(4.10)  lim  Pr(Qo(t)  = qQ}  = P[q Ja^P.  ],  for  qQ  = 0,1,...  . 
t-*» 


4.3,2  At  the  Bases 


To  obtain  Pr{B.(t)  = b}  we  compute  Pr{U.(t)  - s.+b} 
3 3 3 

and  Pr{U.(t)  = s.+b}  by  evaluating  Eqs.  (4.6)  and  (4.7), 
D 3 b 

respectively. 

4. 3. 2.1  Case  (A):  Pr{U.(t)  = s.+b}. 

3 J 


Since  U*(t)  = D?(t-R. , t)  + D?(t,,t)  + D^(t-,t)  and  the 
2 2 2 J 0 J ° 

processes  (Dj(t),  t>0},  (D?(t),  t>0}  and  {D?(t),  t>0}  are 
independent,  Ut(t)  is  a Poisson  variable  with  mean 

J 

+ (Ap  + A^)x_.  • that  is. 


(4 


.11)  Pr(U*(t)  = s.+fc-d}  = P[s_. +b-d|  A?R.  + (A?  + A°)r.]. 
J 3 J 
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* . ? 

•Je  now  proceed  to  find  the  probability  distribution  of  Ur(t). 

3 

Suppose  we  are  given  DQ(t2,t3)  = dQ(t2,t3)  and  D°(t2,t3)  = 
dj(t2,t„).  Then  U?(t)  = d if  and  only  if  d?(t0,t_)-d  units  of 

base  j demands  are  filled  by  time  t . In  other  words,  d.(t„,t,)-d 

O ]]  o 

C 

of  the  first  z (t  ) - d (t  ,t  ) depot  arrivals  after  time  t„ 
u X U X 2 2 

must  come  from  base  j . 

Let  Rfl  3 0>°(t2,t3)  = dj(t2,t3);  00<t2,t3)  = tytj.tj); 

We  shall  obtain  Pr{U?(t)=d|D3(t.|  ,>t2)=dp(t1,t0)-J  Z^t.^  = z^t.^} 

2 A 

by  enumerating  Pr{U.(t)  = d|RA}  over  D (t0,t  ) and  D.(t  ,t„). 

J U 2 ci  J A o 

Because  of  the  properties  of  the  Poisson  process  [20],  we  have 


(4.12)  Pr{D°(t2,t3)  = d?(t2,t3)|D0(t2,t3)  = dQ(t2,t3)} 


(VVV^U°  a 


.0  -|d0(t2’t3)  dj(t2,t3) 

Aj/AoJ 


for  d°(t2,t3)  = 0,1,. . . ,dQ(t2,t3). 


Substituting  Pr{DQ(t2,t3)  = dQ(t2,t3)}  = P[dQ(t2 ,t3> J X0RQ] , in 
Eq.  (4.12)  we  get 


(■..13)  Pr(D°(t2,t3)  = d°(t2,t3);  Dgltj.tg)  = VVV1 

d0(t2,t3)\,.0.  ,dj<t2"t3)  r.  .0.  ,<Vt2,t3)_dj(t2't3) 

v ° ' ■ 3 0 


for  d^  (t2,t3)  £ dQ(t2,t3)(  >_  0). 
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To  obtain  Pr{U?(t)  = d|RA}  we  consider  the  following  two 
mutually  exclusive  cases. 

(i)  DQ(t2,t3)  < z0(t1)  - dg(t1,t2):,  that  is,  all  the  depot 

demands  during  interval  (t2,t3]  are  satisfied  by  time 
Therefore 


(4.14)  Pr{UT(t)  = d | RA}  = 

] I 0 


d = o; 

d > 0 , 


(ii)  DQ(t2,t3)  > z0(t1)  - d3(trt2)  ; that  is,  the  last 

dQ(t2,t3)  - (z0(t1)  - d3(t1,t2))  depot  demands  that  arrived 
during  the  interval  (t0,t,,l  cannot  be  satisfied  by  time  t^. 
Because  of  a first-come,  first-served,  resupply  policy  at  the  depot, 

using  the  results  of  Corollary  1.1  of  Appendix  A,  we  have 


(4.15)  Pr{Ut(t)  = djRA} 

d?(t„,t„) 
1 *■  J 


d0(t2’t3)  " dj(t2’t3) 

0, 


\d? (t2,t3)  - d / \ z0(~t:1)  ~ WV  ~ (dj(t2’t3)~d), 


d0(t2,t3) 
,Z0(tl)  " d0(tI’t2) 


From  Eqs.  (4.13  - 4.15)  wo  obtain, 

(4.16)  Pr{U?(t)  = 0|Dp(t1,t2)  - d^(tl5t2);  Z^)  = zQ(t1)} 


Z0(tl)'d0(tl’t2) 


d(/t2 


PCd0(t2’t3)lX0R0]  ■ 
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d0(  Vt3)=Z0(tl)“d0(tl’t2)  +1 


d0(t2’t3^ 


VvV'^'W  \('/  VVV  Yv 

?<t‘t,)=o]VzO<tl)-dO(^-t2)-d^t2't3)iVo:tl)‘dS<tir2vJ 


j ' 2’  3 


. .¥,2't3,[ 

\d.(t  „,t3>ycvv 


n ,0..  .VW'VW 
li  - yy 


and 


P^do(t2’t3^lX0K0^: 


(4.17)  Pr{U^(t)=d|DQ(t1,t2)=dQ(t1,t2);  Z0(tl)=Z0(tl>} 


l C 

d0(t2’t3)=Z0(tl)_d0(tl’t2)+1 


d0(t2,'!V  f,  , , x 


/a, 

i X 1 


dj(t2>t3)=0^dj(t2»t3)-d 
dQ<t2  *t3) 


d0'^2 3 C j ^2  ’^3^ 

zo(tl^"do(tl’t2^_dj  (t2’t3 \^0*tl^”<Vtl’t2^ 

/d0(t2jt3)N\  d°(t2,t3)  0 d (t  t )-d?(t2,t3) 

Ut2,t3))  3^-s0/v°  3 


P[d0(t2,t3)i  a 0R()3> 


for  d = 1,2,. . . 


As  mentioned  in  Section  4.1,  {DQ(t),  t>0}  is  a Poicson  process 

c 

with  oarameter  Therefore 
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_ c 

(4,18)  Pr{D^(t1,t2>  = d£(t  *,)>  = pCd0(t1,t2)Uo(VR0):i- 


Now  substituting  Lqs.  (4. 11) , (4 . 16) ,(4.17)  a, id  (4.10)  into 
Eo.  (4.6)  wo  obtain  PriU.(t)  = s.+b>  . After  simplifying  (see 
Appendix  F)  we  get 


(4.19)  lim  Pr{U.(t)  = s_.+b)A 


)(ti}  ! 


. . /, C , ,D, 


l ( p[sj+b|x:Rj+(A:txj)Tj) 

c£(tl5tj=:  ■ 


0 1’  2 


L 


z0(t1)_d0(tl’t2) 

l . rtdo(t2’t3)'XOr'0J 


WV*0 


r 0/1  1d0(t2’t3)-Z0(tl)fd 
Cl-XT/A.] 


+ V Vt3)=z0(tl,'d0<tl>t2,+1  ’ ° 

PCSj+b-a|x®Rj»'.>.9,x?)T. 

S d =1 


P[d0(t2 


l 


0 ,,  ,d  P.  ,0 ,,  ,d0(t2’t3)-d"z0(tl)+d0(tl’t2) 


fo^t2,t3^=Z0^tl^  d0'tl,t2 

• tx?/x03d  Ci-x?/x03 

J 3 

p[d0(t2*t3^lXoR0^ 


,t„)+d  \ 


id0(t2,t3)_Z0(tl)+d0!'tl,t2) 


1 c c 

( • ptdo(t1,t2)|X0(x0-R0)]. 


o o 


4. 3. 2. 2 Case  (B):  Pr{U.(t)  = s.+b) 

3 3 5. 


Pore  uj(t)  = D®(t-R^  , t)  + D?(t2,t)  + D?(t2,t);  hence,  uj(t) 

B p d 

is  a Poisson  variable  with  mean  A.R.  + (AT  + A.)(t.  + R_);  that 

3 3 3 3 3 0 

is. 


(4.20)  Pr{U . (t)  = s.+b-d)  = Pi"  s.+b-dj  A .R.  + (A.  + AT)(x.  + Rn)]. 

3 3 3 3 3 3 3 3 0 


To  obtain  the  probability  distribution  of  U.(t)  in  this  case, 
we  follow  the  approach  similar  to  that  used  in  case  (A).  Suppose 

D0(tl’t2)  = d0(tl’V  and  Dj(t1>t2)  = ^VV*  Then  Uj(t)  = d 

Q 

if  and  only  if  d_.  (t^,t2  )-d  units  of  base  j demands  can  be  shipped 
by  time  t^;  that  is,  d?(t  ,t2)-d  of  the  first  zQ(t^)  + dg(t^,t2) 
depot  arrivals  during  the  interval  (t^,t2)  must  come  from  base  j. 
Let  RB  = {^(t^)  = d°(trt2b  D°(t2,t2)  = djJ^.tj); 

D0^tl’t2^  “ d0^tl,t2r  ^0^1^  = z0^tl^* 

Again  because  of  the  properties  of  the  demand  process. 


(4.21)  Pr{D°(t1,t2)=d°(t1,t2)|Dj(t1,t2)=d®(t1,t2);  D^t^)  = 


d0(tl,t2)+d0(tl,t2)  \ 0.  -|dj(tl,t2) 


d.  (tx,t2) 


[A:r/A0] 


r,  o.  1aS(Vt2)rto(Vt2)-dj(Vt2,> 
n - VV 


for  d?(t^,t2)  = 0,1,... ,dp(t1>t2)+dp(t1,t2). 


Also 


(4.22)  Pr{Dg(ti,t2)=d^(t1,t2);  Do(t1»t2)r;do(tl,t2)} 


P^d0^tl’t2^X0^T0_R0^  ' r*'d0^tl,t2^A0^T0-}':0  )‘*' 


Using  Corollary  i . 1 of  Appendix  A 

(4.23)  Pr{Uj (t)  = d | RB} 


dj(t1’t2) 


d0(trt2)+d0<tl’t2)-dj(tl’t2) 


fj(tl’t2)_d/  \Z0(tl)+d0(tl’t2)“dj(tl,t2)+d  / 


1 dOvtl’t2^H'do''tl’t2^ 


V Z0(tl)+d0Ctl’t2) 


Multiplying  Eqs.  (4.21),  (4.22)  and  (4.23)  and  then  enumerating 

over  D?(t1,t2)  and  D^t^,^),  we  obtain  Pr{U?(t)  = d|Dp(tl5t2)  = 
r 

dcCtfst^);  ZQ( t^)  = zQ(t  )}.  Substituting  this  probability  and 
Eq.  (4.20)  into  Eq.  (4.7)  (see  Appendix  B)  we  get 


(4.24)  lim  Pr{U.(t)  = s.+b}n 

1 1 3 

+->00  J J 


'l  P[s.+b-d|X®R.  + (X?+XI?)(T.+Rn)] 
-n  3 3 3 3 3)0 


d0(tl’t2)::Z0(tl)+d 


d0(tl’t2)=° 


d0(tl’t2)‘Z0(tl) 
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rx>x0]  [i-A3r/A0] 


d_(t.,t0)-z„(t.)-d 


0„  0 1’  2 O'^l'  r .D,  . N , , D,  _ 


PCd0(tl't2>l>0<VV1 


Substituting  Eqs.  (4.19),  (4.24)  and  (4.8)  into  Eq.  (4.3) 
ve  obtain  i'r  { ( * ) = b). 

For  the  case  where  the  depot  follows  an  (Sq-1,Sq)  policy, 
Eq  = {SQ}  and  II(Sq)  = !•  Pr(3.(*)  = b}  can  now  be  obtained  by 

substituting  zQ(t. ) = S into  Ecs.  (4.19)  and  (4.24). 


4 .4 Special  Cases 


4.4.1  Complete  Recoverability 

C 

For  the  case  of  complete  recoverability,  P = 1;  therefore  X^.  = 0 

( j=0 ,1, . . . J)  and  X®  = (l-r.)X.  = X?  ( j=l,2 , . . . ,J) . Consequently 
D J 1 3 9 9 

X . = X.  = J (l-r.)X..  Since  there  are  no  condemnations  and  the 
0 0 jfx  3D 

system  is  conservative,  no  procurement  is  made  by  the  depot  from 
the  external  supplier.  The  inventory  position  at  the  depot  remains 
at  a constant  level  5^  (say),  that  is,  ZQ(t)  = SQ,  for  all 
t ^ 0.  The  stationary  distribution  of  (BQ(t),  t>0}  is  obtained 
in  a manner  similar  to  that  discussed  in  Section  3.6.1  and  is  given 
by 


lim  Pr{BQ(t)  = 
t-+=° 


V ■ p[so  * bolxoV 


for 


8 -V  -A'o 


+ Is ♦ • • 
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Since  D^(t)  = 0 for  all  t 0,  case  (B)  discussed  in 

section  4.2.2  will  not  arise  and  consequently  Eq.  (4.3)  reduces  to 

Pr(u.(t)  = s.+b}  " Pr{U.(t)  = s.+b}  for  z_(t1 ) = S . Upon 
•j  3 3 D A jiu 

substituting  X^  = X0  = X0’  d0^tl’t2^  = ° and  zo'tl^  = ^0 

Eq.  (4.19),  we  obtain 


(4.25)  Pr{B_.  (*)  = b}  =P[sj+b|X^Rj  + X?t.] 


. J t ,.0  «VVVl»oV 

d0(t2*V'° 


. I [1  - l°/J  ]d°<t2’t3>  S°  PWQ(t  ,t  )|A0R0) 

WV'V1 


s_.+b 

+ l P[s.+b-d|X?R.  + X^t-] 
d=l  ] ] : 


d0^t2’t3^  ^C+d  \ d 


VVV^oV  -,d 


0 d')(t2’t3)_d"S0 

• [1  - X°/XQ]  0 2 3 ° FLd 


1 

0(t2’t3) iX0R0] 


We  note  that  Eq.  (4.24)  is  equivalent  to  the  result  obtained 
by  Simon  [20]  for  the  case  p = 1. 


4.4.2  Complete  Non-recoverability 

Wien  an  item  is  consumable,  p = 0 and  r^  = 0 ( j=l,2 , . . . ,J). 

The  repair  loop  is  absent  at  each  location  in  the  system  and  Q^(t)  = 0 
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for  t >_  0 and  j = 0 , 1 , . . . , 3 
{Z  (t),  t>0}  and  (D0(t),  t>G) 

respectively.  Setting  r,  = 0 
and  using  Eq.  (4.3)  we  obtain 


The  stationary  distributions  of 

are  given  by  Eqs.  (4.8)  and  (4.9) 
J 


and 


in  Eq.  (4.24), 


(4.26)  Pr{Bj*)  = b} 


zo(ti)eEo 


P[s^+b|X.T.] 


zo(ti} 


1 PCd0('t2’t3)  |X0R0^ 


^0(t2’t3)  ° 


d0(t2’t3)=Z0(t1)+1  J 


d (t  ,t  ) ~z  (t  ) . 

C1-A./X.3  0 2 3 0 1 rt^VVlVc 


s.+b 

-% 

+ T Pts.+b-djX.T.2 
0 i 1 i i 

d=l  J J 


(wV-o  “Aa.nJ 
lV*!'t3S(tlW  V d / ’ 

(t  ,t,)-d-Z  (t  ) J 

• I“W0<t2,t3)l  oVl 


u - yy 


!J 


(S0  so) 


CHAPTER  V 


TWO-ECHELON  SYSTEM-RANDOM  ORDER  SIZE 


5.1  The  Model 

In  this  chapter  we  study  the  two-echelon  system  as  described 
in  Section  4.1.  Requisitions  arrive  in  a Poisson  manner  with  known 
parameter  at  base  j (j=l,2,... ,J).  Upon  arrival  of  a 

requisition,  a batch  containing  one  or  several  failed  units  is 
turned  in,  and  a like  number  of  new  units  is  demanded  for  replacement. 

We  shall  consider  both  batch  and  unit  models  for  the  inspection  of 
failed  units.  In  the  batch  model,  a batch  as  a whole  is  either  base 
repairable,  depot  repairable  or  condemnable,  whereas  in  the  unit 
model,  each  unit  in  a batch  is  inspected  independently  to  find 
whether  it  is  base  repairable,  depot  repairable  or  condemnable.  We  shall 
use  the  same  assumptions  about  procurement  policies,  repair  facilities, 
repair  and  lead  times,  as  specified  in  Section  4.1.  In  addition, 
partial  backlogging  of  the  demands  is  allowed-  that  is,  upon  arrival 
of  a requisition,  if  the  base  does  not  have  the  number  of  units 
demanded,  then  all  the  units  on  hand  are  supplied  while  the  balance 
is  backlogged.  Partial  backlogging  is  also  allowed  at  the  depot. 

We  shall  use  the  following  nomenclature  in  addition  to  that 
introduced  in  Section  4.1 

Nj(t)  = total  number  of  requisitions  that  arrived  at  location 
j during  the  interval  (0,t]  (j=0,l,. . . ,J). 
j ( k ) = Pr{upon  arrival  of  a requisition  at  location  j,  total 
number  of  units  demanded  = k},  k>l t (j=0,l, . . . , J) . 
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4> j C k ) = Pr{upon  arrival  of  a requisition  at  base  j,  the 

number  of  base  repairable  units  = k}  , k>_l , 

( 3 ~ 1 , 2 , • • • , J ) • 

C 

<)>j(k)  = Pr{upon  arrival  of  a requisition  at  location  j, 
the  number  of  condemnable  units  = k ) , k>l, 

(j=0,l, . . . ,J). 

<j>?  (k)  = Pr{upon  arrival  of  a requisition  at  base  j , the 
number  of  depot  repairable  units  = k },  k>_l , 
(j=l,2,...,J). 

00  , .n  / x 

CP[k|At,f]  = £ - ^ — — r n vk)  (compound  Poisson  distribution 

n=0 

with  parameter  A and  compounding  distribution  f). 

In  Section  5.2,  we  consider  the  inspection  of  the  failed  units 

•under  the  batch  model.  In  Subsection  5.2.1,  we  derive  the  results  for 

the  case  where  order  size  distribution  is  the  same  at  all  bases",  that  is, 

<!>,(•)  = $0(‘)  = ...  = <}> T ( • ) . Subsection  5.2.2  examines  the  case  of 
A.  Z J 

different  order  size  distribution  at  the  bases.  In  Section  5.3,  we  consider 
the  inspection  of  the  failed  units  under  the  unit  model. 


5.2  The  Batch  Model 

In  this  model,  upon  arrival  of  a requisition  at  base  j,  the 
entire  batch  of  failed  units  is  repaired  at  the  base  with  probability 
r j , is  shipped  to  the  depot  for  repair  with  probability  (l-r.)p,  or 
is  condemned  with  probability  (l-r^Hl-p).  This  divides  the  requisitions 
into  three  types:  base  repairable,  depot  repairable  and  condemnable. 


For  t>0  and  j=l,2,...,J,  let 


1 - *3 i. 


and 


&4 


IJ. (t)  = the  number  of  requisitions  at  base  j during  the  interval 
(0,t’J  for  which  the  entire  batch  was  declared  base 
repairable, 

c ... 

fL(t)  = the  number  of  requisitions  at  base  j during  the  interval 
(0,t]  for  which  the  entire  batch  was  condemned. 


N.(t)  = the  number  of  requisitions  at  base  j during  the 

interval  (0,t]  for  which  the  entire  batch  was  sent 


to  the  depot  for  repair. 

Obviously,  N.(t)  = bP(t)  + Np(t)  + hP(t),  for  all  t>0. 

3 'J  j : 

Following  the  arguments  used  in  deriving  Eq.  (3.52),  we  see  that 

the  processes  (N?(t).  t>0},  {N?(tb  t>0}  and  {N?(t),  t>0)  are 
1 - 3 — 3 

mutually  independent  Poisson  processes  with  parameters 


= r.*.,  = (1-r. )(l-p)A.  and  = (1-r. fox.,  respectively 

3 3 3 3 3 3 3 3 J 

B C 

( j = l.  2 , . . . ,J ) . The  demand  processes  {D^Ct),  t>0},  {D_.(t),  t_>0}  and 

{ Dp  ( t ) .,  t>_0}  are  compound  Poisson  processes  with  parameter 

Xj,  Ap  and  A?,  respectively,  and  have  a common  compounding  distribution 

*jC). 

It  is  clear  from  the  above  that  the  depot  receives  two  distinct 

types  of  requisitions  from  each  base.  One  type  requires  depot 

repair  for  the  entire  batch,  while  the  other  corresponds  to  a 

condemned  batch  of  items.  For  t>0,  let 

N?(t)  = total  number  of  requisitions  placed  at  the  depot  by 

base  j during  the  interval  (0,t]  ( j=l ,2 , . . . ,J) , 

C 

bQ(t)  = total  number  of  requisitions  at  the  depot  during  the 
interval  (0,t]  as  a consequence  of  condemnations. 


»D 


and 
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N^(t)  = total  number  of  requisitions  at  the  depot  during  the 
interval  (0,t]  for  which  the  batch  of  failed  units 
was  found  depot  repairable. 

Because  the  bases  use  an  (s-i,s)  policy,  N?(t)  = N (t)+Nt’(t),  for  all 

t>0  and  j=l,2,...,J.  Consequently,  {N?,  t>0}  is  a Poisson  process 

with  parameter  A?  = A?  + A^.  Furthermore,  because  the  bases  operate 

independently,  {N^(t),  t>0}  and  (N°(t),  t>0}  are  Poisson  processes 
C J C D J D 

with  parameters  A„  = J A.  and  \ = J A.,  respectively. 

0 j=l  3 0 j=i  3 

Obviously,  Ng(t)  = N^(t)  + N^(t),  for  all  t>0.  Thus  (N^Ct),  t>0} 

^ 0 

is  a Poisson  process  with  parameters  An  = £ A..  Using  these  results, 

j=l  ] 

we  now  proceed  to  find  the  probability  distributions  for  the 
processes  {DQ(t),  t>0},  {D^(t),  t>0}  and  (D^(t),  t>0}. 

Let  ij).(w)  be  the  characteristic  function  of  that  is, 

CO  3 t 

ijj.(w)  = l 4!.(k)elwk.  Then  <Jj?(w),  the  characteristic  function  of 

o k=1  3 3 

Dj(t),  is  given  by  [18], 


(5.1) 


ifjlw)  = e -1  ] 


Thus  ( D j ( t ) , t >0  } is  a compound  Poisson  process  with  parameter  A? 
and  compounding  distribution  <^(-).  Since  the  bases  operate 
independently,  DQ(t)  is  the  sum  of  J independent  compound  Poisson 
processes.  Therefore,  ^.^Cw),  the  characteristic  function  of 
Dfl(t),  is  given  by 


8 


V<t>(“) 


= P i(j.(w) 
j=l  3 


J -A?t[l-Mw)] 

= n e 3 3 

j=l 

-x0ttl-r'  E x?Vw):i 

= e A0  j.l  ] ] 


Thus-  {DQ(t),  t>0)  is  a compound  Poisson  process  with  parameter 
X p and  compounding  distribution  i+,^(  • ) whose  characteristic  function 
i^0(w)  is  given  by 


D»0(w)  = r~  I X%.(w). 
0 j=l  3 3 


From  the  additive  property  of  characteristic  functions , it  follows 
that 


(5.2) 


,0, 


j J 

k(k)  = i xrt.(k), 

0 xo  j=l  33 


k>l. 


c 

Similarly,  we  can  show  that  the  demand  processes  {DQ(t),  t>0)  and 

{D°(t),  t>0}  are  independent  compound  Poisson  processes  with  parameter 

CD. 

XQ  and  Aq,  respectively.  Their  respective  compounding  distributions 
are  given  by 


(5.3) 


oo  * 4-  I 


»o  3 3 


k>i. 


and 


(5.4) 


♦>>  * ~F  j,  lh(k): 


k>l. 


The  depot  can  now  be  analyzed  as  a single  location  system 
where  recoverable  and  non-recoverable  demand  processes  are  independent 
compound  Poisson  processes.  Therefore,  the  results  derived  in 
Section  3.4.1  apply.  From  Eqs.  (3.13),  (3.46)  and  (3.20)  we  have 

m(S-k) 


(5.5)  lim  Pr{Z0(t)=k|Z3(0)=i,  i>s  }=nQ(k)  = 

■£-*co 


l+M(S-s-l) 


1 

1+M(S- 


s-1) 


s+l<k<S-l, 


k=S; 


wnere 


c r k-i  k 

m(l)=c|>  (1),  m(k)=<J;  (k)+  ][  4>  (k-q)m(q),  k>2,  and  M(k)  = £ m(H),  k>l; 

q=l  Jt=l 


(5.6)  lim  Pr{Q0(t)  = q0|Q0(0)=0}  = CPtqjA^,  qQ>0; 


and 


(5.7)  lim  Pr{BQ(t)  = b0|zc(0)  = i,  i>s0> 


fk+b 


I (l  “ cp[ktb0-dc|x°R0,  ^]crtdc!xc0 

kfcto  ' 


’ $Wk) 


bg  ~Sq ,-Sq+1 , . . . , 0 , 1 , 


b 3 


Let  U3  consider'  a subgroup  of  the  bases,  containing  more  than 
one  base.  Then  following  the  ar'guments  used  in  deriving  Eqs. 

(5. 2-5. 4)  wa  can  easily  determine  the  probability  distributions 
for  the  resultant  demand  processes  of  this  subgroup.  For  the 
purpose  of  obtaining  the  stationary  distribution  of  the  backorder 
process  at  base  j , we  view  the  demand  at  the  depot  arising  from 
two  sources.  One,  the  base  j for  which  the  distribution  in 


being  determined,  and  the  other  being  set  of  the  remaining  bases. 

Let  us  denote  this  set  by  o,  that  is,  o = {l,2,...,j-l,j+l,...,J}. 
Since  the  bases  operate  independently,  the  two  sources  are 
independent.  For  source  o,  u ^ shall  use  notation  similar  to 

that  used  for  an  individual  base.  The  processes  {N‘  (t),  t>0}, 

C D 

{ ( t ) , t>0}  and  {i!o(t),  t>0}  are  mutually  independent  Poisson 


processes  with  parameters  A^  = J A*?,  A^j  = I A^  and 


B 


A^  = £ aY,  respectively.  Following  the  arguments  used  in  obtaining 

ieo  1 

Eqs.  (5. 2-5. 4)  it  can  be  easily  seen  that  the  demand  processes 

(D^(t),  t^o),  {b^(t),  t^O}  and  (D^(t),  t^O)  are  compound  Poisson 

B C D 

processes  with  parameters  Aq,  A^  and  A , respectively.  Their 


ICO 


ICO 


, B. 


i->  J.  r B 

respective  compounding  distributions  are:  ^(k)  = — — > A.$.(k>, 

o -i  D .“'ll 

Aa  ieo 

C(k)  = 4r  l A?ij>.  (k)  and  yD(k)  =—  £ A°<^ . (k) ; k>l.  Also, 

a , C . 1 1 , D . x 1 — 


A ieo 
a 


A*  ieo 
a 


. M-jU 
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the  process  (N^(t),  t>0 } is  a Poisson  process  with  parameter 
A"1  = i A?  and  the  demand  process  (D^(t),  t>0)  is  a compound 
Poisson  process  with  parameter  A^  and  compounding  distribution 

I A°di(k),  k>l. 

A ieo 
o 

In  the  next  two  subsections,  we  obtain  the  stationary 
distribution  for  the  number  of  backorders  at  the  bases  using  the 
approach  as  described  in  S-ection  4.2,  In  Section  5.2.1,  the  results 

are  derived  for  the  case  where  the  order  size  distribution  is  the 

sane  at  all  bases,  while  in  Section  5.2.2,  the  case  of  different 

order  size  distributions  at  the  bases  is  examined. 

5.2.1  The  Same  Ord^.r  Size  Distribution  At  The  Bases 

Let  $(•)  be  the  common  order  size  distribution  at  the  bases. 
From  Eqs.  (5. 2- 5. 4),  it  follows  that  </>(•)  = 4>q(')  = • ) = $(•); 

that  is,  the  deirand  processes  at  the  depot  have  a common  compounding 
distribution,  <)>(•).  The  stationary  distributions  of  the  processes 
(Zp(t),  t>0},  (Q^Ct),  t>0}  and  {BQ(t),  t>0}  can  be  obtained 

from  Eqs.  (5.5),  (5.6)  and  (5.7),  respectively , upon  substituting 

CD  B 

<f>Q(  * ) = ^q^-)  = ('(•).  Also,  we  note  that  in  this  case,  $ ( • ) = 

♦Jit*)  = ♦)!(•)  = ♦<•). 

o o 

To  obtain  the  stationary  distribution  of  the  backorder  process 
(Bj(t),  t>C},  we  first  compute  Pr{lK(t)  = s^+b).  In  order  to  do 
so,  we  consider  the  two  cases  A and  B as  described  in  Section  4.2. 

We  recall  that  case  A corresponds  to  the  situation  where  the  total 
depot  demand  during  the  interval  (t^jt^]  (see  Figure  4.1)  does 
not  exceed  the  total  assets  available  at  the  depot  by  time  t^. 


r 


Case  13,  cn  the  other  hand,  represents  the  situation  wh'.re  the  total 

depot  demand  daring  the  interval  (t^,t^]  exceeds  the  total  stock 

available  by  time  t0.  The  probability  distribution  of  U.(t)  for 

o 3 

the  case  A and  B arc  derived  in  Sections  5. 2. 1.1  and  5. 2. 1.2, 
respectively. 


5. 2.1.1  Case  A:  Pr{U.(t)  = a . tb } 

3 1 A 

As  .mentioned  in  Section  4.2,  in  case  A all  the  depot  demands 

during  the  interval  (t  ,t„]  are  satisfied  by  time  t . Only  the 

depot  demands  during  the  interval  (t  ,t  ] may  remain  unsatisfied 

Z o 

by  time  t . To  obtain  {Pr  U.(t)  = s.+b}  in  case  A,  we  evaluate 
J 3 3 

Eq.  (4.6)  by  computing  the  probability  distributions  of  lh(t)  and 
2 

U.(t),  the  two  independent  components  of  U.(t). 

3 3 

By  definition,  Ud(t)  = D®(t-R.,  t)  + D?(t  ,t).  Thus  the 

3 3 3 3 o 

random  variable  Lh(t)  has  a compound  Poisson  distribution  with 
B CD 

parameter  + + ^j^Tj  and  comPounding  distribution  <{>(•). 

Therefore, 


(5.6)  Prttrht)  = s.+b-d}  = CP[s  . +b-d  | A^R . + ( + A^)t  , <*>]. 
3 3 3 33333 


The  random  variable  LL(t)  represents  the  units  ordered  from 

the  depot  by  base  j during  the  interval  (t^t^l  'that  remain 

2 C 

unfilled  by  time  t^.  We  shall  obtain  Fr{LL(t)  = djD^t^jt^)  = 

c 0 

dC^tl’t2)’  Z0^1^  ~ zo(tl)}  by  conditioning  on  Nj/t^tg),  i=j  , o . Let 


RA={N.(t2,t3)=nj  ; ^0(t2’t3)  =no;  V V V = VVV ; VV  = Z0(V ' 


m »•  +< 
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To  obtain  Pr{uf(t)  |PA},  we  need  to  know  the  number  of  requisitions 
at  the  depot  placed  during  (t^tg]  and  completely  satisfied  by 
time  t . Let  N.'(t0,t  ) denote  the  number  of  requisitions  from 

O 1 L.  O 

source  i that  arrived  during  (t  ,t_]  and  are  completely  satisfied 

2.  o 

by  time  t , i=j,  a. 

O 

Now,  suppose  we  are  given  RA  and  N ! ( t ,t  ) = n'. . Then 

3 / o 3 

2 

(h(t)  = d if  and  only  if  the  sum  of  the  demands  due  to  the  unsatisfied 
-n!  requisitions  and  unsatisifed  units  of  possibly  a partially 
satisfied  requisition  (if  from  source  j)  equals  d.  Let  EX  denote 
the  number  of  the  units  supplied  to  the  requisition  whose  demand  is 
only  partially  met.  The  range  of  the  random  variable  EX  is  from 

<3 

0 to  Zg(t^)  - d^t^,^).  When  EX  = 0,  there  is  no  partially 

C 

satisfied  requisition  and  when  EX  = z3(t^)  “ c*o'tl,t2^’  no 

C 

requisition  is  completely  satisfied  and  zQ(t^)  - d0(t^5t2)  units 
are  supplied  to  the  first  requisition  (if  any)  (see  Figure  5.1) 

Let  us  introduce  an  indicator  variable  I such  that  I " i 
if  the  partially  satisfied  requisition  is  from  source  i,  i=j,  o . 

We  have  the  following. 


(5.9)  Pr(U?(t)  = d |RA;  EX  = 0;  N.'(t2,t3)  = nj;  N^t^tg)  = n^} 
(n.-nj) 

= 4>  3 J (d); 


(5.10)  Pr{U?(t)  = djRA;  EX  = e > 0;  N^V^)  = n^';  N0'(t2»t3)  = n^ ; I = o) 


= <t* 


(nj-nj)(d); 


(5.11)  Pr{UT(t)  = d | RA;  EX  = e >0;  = n . ; N (t  ,t  ) = n • I = ji 

3 l*o  1020  o 

d (n.-n.-l) 

= £ <j> ( k+e ) • 4>  3 -1  (d-k). 

k>0 


The  Pr{U..(t)  = d|RA}  can  be  obtained  from  Eqs.  (5.9  - 5.11)  by 

knowing  Pr{EX  = e(>  0);  I=i;  N.(t  ,t  )=n.;  N (t  ,t  )=n  IRA}, 

for  i = j,o.  Let  Y^  denote  the  number  of  units  demanded  from 

the  depot  by  the  kth  (k_>l)  requisition  during  (t^tg].  Using 

the  results  of  Theorem  1 of  Appendix  A,  we  obtain  the 

Pr{EX  = e,  I=i,  N.(t  ,t  ) = n.;  N (t  ,t  ) = n | PA}  . We  consider 
3^0  3 ^ ^ ^ ^ 

the  following  two  situations.  One  where  not  all  the  requisitions  during 
(tjjtg]  are  satisfied  by  time  t^  ; and  the  other,  where  they  are  all 
satisfied. 

i i 

(i)  0 < n.  + n < n.  + n : n.  + n >1,  not  all  the 

- i q 9 0 3 ° - 

requisitions  in  (t  ,t  ] are  satisfied  by  time  t . 

2 3 J 


(5.12)  PrtEX  = e;  I=j;  N.ft^tg)  = n ; Na(t2,t3)  = nQ|RA} 

= Pr((Y  +Y  +...+Y  , ,)  = z (t  ) -d3(t1,t2)  - e;  Y , , > e; 

n.+n  n.+n  +1 

3 ° 3 0 


out  of  the  first  (n!  + r/  ) requisitions  at  the  depot. 

n'.  are  from  source  i,  and  the  (n!  + n'  + l)th  requisition 
3 jo 

is  from  source  i I R"A } 
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. n.  . n x 


n .+1 
_J —2 


"a 


(n.+n  ) 


n.+n  +1 

3 o 


n.+n  \ 

’D  ' ) 
n . +n  +1 / 

T rv  ' 


i 0 [ 


St>e 


3 o 


for  e = 0,l»....z0(t1>-«io(t1»t2)  “ (nj  +nj);  0 inj  ln; 


Similarly, 


» ii  1 , . 

(5.13)  Pr{EX  = e;  I =o;  (tg.tg)  = n.. ; iMt^tg)  = nJRA} 


' 11  j \ / na  \ 


na+! 


n : 
D 


U{no+1> 


» t 

(n  .+r.  ) 


n . +n  +1  (n.+n  ) 

-j  a 3 a 


° (z0(tl)"d0(tl*t2>-e)  * £ 


ft>e 


for  e = 0,1,.  ..,z0(t1>  - dQ(tl5t2> -(n. +n^;  0 <_n  <_i 


Summing  Eqs.  (5.12)  and  (5.13)  we  have 


(5.14)  Pr{EX  = ej  N.(t  ,t,) =n,;  N^(t0,t,)  = nn | RA} 


j’  a 2 3 ‘ 


( no  . ■ 

' - ' - (n.+n  ) 


/ n . + n \ 

unn) 


— $ \ 


z0(tl)  "d0(tl,t2)  " e)'  ^ ‘J’O  , 

£>e 


for  e = 0,1, . . . ,Zq(1i^)  - » 


i i 

and  0 < n.+n  <n.+n  . 
-303a 


<p(z). 


4.00, 


(5.15)  Pr{EX  = e;  N^t^tg)  = n^;  VW  = nQ  [RA} 


= Pr{(Y1+Y2+. ..+Y  , , ) <_  zQ(t1)  - d3(t1,t2)} 

n.+n 
3 a 

f^0(tl,“d0<tl’t2}  (n!+r/) 


l 


i i 

Z= n .+n 
j a 


la'  ' ' 

J (Z\  for  e = OjTiy-nyn^n^, 

n.+n  >1; 

] c- 


for  e = zQ(t1)  - dQ(t1,t2), 

I » 

n.  = n = n.  = n = 0; 
3 a J a 


, otherwise 


From  Eqs  . (5.9  - 5.1i)  and  Eqs.  (5.12  - 5.15)  we  obtain 

2 

Pr{Ut(t)  = d|RA).  We  consider  the  following  two  cases:  one  where 

all  the  demands  from  base  j during  (t^tg)  are  satisfied  by 
time  t (d=o);  and  the  other, where  some  demands  from  base  j during 

O 

(t2,t3J  remain  unsatisfied  by  time  t3  (d^l). 

2 

(i)  Pr{LK (t)  = 0 |RA}  . From  Eqs.  (5.9)  and  (5.10)  we  conclude  that 

O * f 

given  RA,  Ut(t)  = 0 if  and  only  if  n_.  = n^(>_0).  For  the  case  when 
all  the  requisitions  at  depot  in  (t0,t33  are  satisfied  by  time  t3, 


that  is,  n.  + n = n.  + n } then  Pr{U?(t)  = 0 1 RA } = 1.  From 
3 ° 3 ° 1 

Eq.  (5.15)  we  have 


(5.16)  Pr{Uj(t)  = 0 |RA}  = / z0(t1)-do(ti’t2) 


, for  n . + n =0 

3 0 


0 1 0 1 " ( n . +n  )'  for  n . + n >1 

I ♦ 3 ° 3 °~ 


?.=n.+n 

1 


On  the  other  hand,  when  n_.  + n^  < n^  + ncj  that  is,  not  all 
requisitions  at  tne  depot  during  (t0,t^]  are  satisfied  by  time 

i i 

t_,  then  n.  = n.  implies  that  n < n . From  Eq.  (5.13),  after 
3 3 J on 

' -olifying  the  combinatorial  expressions  we  have 


(5.17)  Pr(lMt)  = Oj  RA} 


• r 


V1  ng!(nj*nq)!  JZC(tl)_d0(tl’t2)_(nj+na) 

\ n !(n.+n  )’.  , \ 

« =n  0 3 0 V ~ = 0 


(n.+n  ) 

4,  ^ 0 (z0(t1)-dQ(t1,t2) 


l ♦<*>> 

H>e 


(ii)  Fr{U?(t)  = d jRA},  d_>l  From  Eqs.  (5.-3)  and  (5.10)  it 

is  clear  that  given  RA,  U?(t)  = d > 0 if  n^  < ik  and  thus 

n + n < n.  + n for  n.  >1.  From  Eqs.  (5.9  - 5.11)  and  (5.12 
j a 3 a 3 “ 

5.14),  it  follows  that 
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(5.18)  Pr£Ut(t) 


n. -1  n 

d|RA)  = ) l 
n.  = 0 n =0 
3 a 


> n . v / n \ 

j )(;) 
' n n ' 


n . 
3 


'n . + n 
3 a 

i i 

a.  + n 
3 


r 


(“0(tl)“d0(tl‘t2)-(nj+na)  ( n ! +n ’ 


V. 


I 

e=0 


3 a 


(20(tl)"d0(tl’t2)_e) 


n ,-n . 
J 3 


d (n_.-n.-l)  'N 

, . t i \ <j>(He)<i>  J J (d-k)  i 

W(n  +n03-Uj+n0)  / 


n -n 
a a 


. ( r. . +n  )-(n.+n  ) 
3 0 3 u 


(n,-n!)  \| 

6 -1  -1  (d)  l (p  (£)))>. 


£>e 


J 


J 


2 C 

From  Eqs.  (5.16),  (5.17)  and  (5. IS)  we  obtain  Pr{L\(t)  = d j 
= d^(t  , t.) j Zn(t.,)  = z„(t  )}  by  enumerating  over  N.(t  ,t  ) 

U _l  2 U u,  U ± jzo 

and  li  (t_.t„).  Substituting  the  resulting  expression  and  Eq.  (5.8) 
a 2 3 


into  Eq.  (4.6)  we  obtain 
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5. 2.1. 2 C-F-'j  (D):  Pr(U.(t)  = s.<-b}n 

1 1 B 

In  this  case,  none  of  the  depot  de.nands  during  (t  ,t  ] are 

Z o 

filled  by  time  t . The  sxock  of  z (t  ) + d^(t  ,t0)  units 
o U J.  U J.  Z 

available  at  the  depot  by  time  t determines  how  many  requisitions 

o 

at  the  depot  during  (t  ,t  J can  be  filled  by  time  t . To  obtain 

-L  Z O 

Pr{U^.(t)  = s_.+b}  in  case  B,  we  evaluate  Eq.  (4.7)  by  computing 

the  probability  distributions  of  U^(t)  and  UZ(t)  as  described 

3 

in  Ssction  4.2.  Here,  U^(t)  = D^t-R.,  t)  + D^t^.t);  therefore, 

3 3 J 3 ^ 


(5.20)  Pr{U^(t)  = s_.+b-d}  = CP[s^+b-d  j A^R..  + A°(t..  + R0), 

Ut(t),  in  this  case,  represents  the  units  ordered  from  the  depot 
by  base  j during  that  31,0  unfilled  by  time  tg.  We  shall 

obtain  Pr{U?(t)  = dl^t^t^  = ; Z^)  = z^)}  by 

conditioning  on  N?(t^,t2),  N^(t,  ,t2)  and  D_.(t^,t2),  i = j,o  . We 

c c 

temporarily  denote  1^(1^, t^)  by  d^.  Let 

RBI  = (N^(t1,t2)  = nj;  N°(tx,t2)  = n°;  ^(t^)  = n°;  = nD; 

zo(ti)  = 2o(ti)}’ 

RB2  = {D^(t1,t2)  + D^(t1,t2)  = d^;  ^(t^)  + D^,^)  = d°} , and 


RB  = RBI  \J  RB2. 

Then 

(nC+nC)  c (n^+n^)  D 

(5.21)  Pr{RB2|RBl}  = <P  3 ° (dQ)  * 4>  1 (dQ). 

To  obtain  Pr{U?(t) |RB},  we  need  to  consider  the  number  and  type 
(depot  repairable  or  condemr.able)  of  requisitions  placed  at 


IOC 


the  depot  daring  (t^,t2]  from  each  source  and  completely  satisfied 


by  time  t . For  i = j,a,  let 
3 

' C 

N^  ’(t7 ,t = the  number  of  condemnation  type  requisitions  that 

arrived  at  the  depot  during  (t^,t2]  and  are  completely 
satisfied  by  time  t^; 

and 

' D ... 

N.  (t, ,t„)  = the  number  of  depot  repairable  type  requisitions 
1 x 2 

that  arrived  during  (t^,t9]  and  are  completely 
satisfied  by  time  t^. 

'C,  . 'C 

Now  suppose  we  are  given  RB  and  know  that  ( t^ , ; = n^ 

and  N* D(t  ,t  ) = n'.D,  i = j,0.  Then  U?(t)  = d if  and  only  if 
1 1 2 1 3 

. ^ . C D 'C  'D. 

the  sum  of  the  demands  due  to  the  unsatisfied  (n  ^ + n..  - n ^ - n_.  ) 

requisitions  and  unsatisfied  units  of  possibly  a partially  satisfied 

requisition  (if  from  source  j ) equals  a.  Let  EX  denote  the  number 

of  units  supplied  to  the  partially  satisfied  requisition.  For 

i = j,a,  let 

iC,  if  the  partially  satisfied  requisition  is  of 


condemnation  type  and  is  from  source  i. 


i- 


iD , if  the  partially  satisfied  requisition  is  of  depot 
repairable  type  and  is  from  source  i. 

Then  similar  to  Eqs.  (5.9),  (5.10)  and  (5.11)  we  have 


(5.22)  Pr{U? (t)  = d|RB;  EX  = 0;  (If.0^,^)  = n^;  NiD(t1»t2)  = nj.D’  i= ^ ’ o)} 


= <t> 


, C D 'C  'D, 
(n.+n.-n.  -n.  ) 

1 3 3 ] (d); 
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(5.23)  Pr{uf (t)  = d Jr3,  EX  = e > 0;  = n^,  N.0^,^)  = n^i 

i = j , a)-,  I = 6 } 


, C D 'C  * D. 
(n.+n.-n.  -n . } 

i 1 3 3 3 (d). 


for  6 = aC,  oD^ 


and 


(5.24)  Pr{uf(t)  = d|RB,  £X  = e > 0;  (N.^t^tg)  = n^0, 

i = j , o) , 1=6} 


, , C D ’C  'D 

d (n.+n.-n.  -n.  -1) 

= l <J>( k+-e ) <t>  3 ^ J ^ (d-k),  for  6 = jC,  jD. 

k>0 


, . , *C  'D  'C  'D  C D C D 

In  this  case,  we  note  that  0<n.  +n.  +n  +n  <n.+n.+n  +n  . we 

-3  3 a o ] 3 o a 

1 C 1 c 1 D 

no;/  need  to  obtain  Pr\EX  = e,  (N^  (t^t^)  = ru  ; (t^tj)  = 

* D 

n.  , i ! j,  o)j  I = 6|RB}.  This  will  be  done  following  the  approach 

used  to  obtain  Eq.  (5.12),  and  using  the  results  of  Theorem  2 and 

Corollary  2.2  of  Appendix  A.  Let  (Yi ,Y„, . . . ,Y  q)  be  the  sequence 

J-  ^ n . 

1 C ... 

of  the  number  of  units  demanded  from  the  depot  by  the  n^  requisitions 

from  source  i (i  = j,  o).  Similarly,  let  {Y  ,Y  , . . . ,Y  _ } be  the 

ni 

sequence  of  the  number  of  units  demanded  from  the  depot  by  the 
n*?  requisitions  from  source  i (i  = j,  o). 


Let 


PS  = 


, c c. 

(n .+n  ) 


where 


'C  C ' D D 

n,  n <n  (_>1);  k = j ,o ; and 


n V ' n i n in 

e = 0,1, . . . ,z  (t  ) + d - (n . +n.  +n  +n  ). 

0 1 0 3 3 o a 


The  ranges  for  and  k 2 are, 


in  'n  *r  «r  c r r • r 1 r 

max(n.  +n  , n.  +n  -z  (t, )+e)}  < k,  <d'T-z.(t,  )+e-(n.+n  -n.  -n  ), 
lo  3 a 0 1 —1—001  303a 


,D  , D D ’D  ’ D. 

l<k  <d„-k  -e-  n.+n  -n,  -n  ). 
— 2—01  3°!  o 


The  probability  distribution  of  U^(t)  can  now  be  obtained 


from  Eqs.  (5.20  - 5.26).  Let 


PS1  = 


. C D C D.  , ’C  1 D 'C  'D, 

(n.+n.+n  +n  ) - (n.  +n.  +n  +n  ) 
] 3 o o 330a 


s .+b 

(5.27)  Pr{U.  (t)  = s . +b } = V "V  PrtU^t)  = s.+b-d} 

3 3 b „ 3 3 

dc>zo(ti}  d=0 


l l I l l l l l l l 

D C D C 'D  'C  'D  'C  iPSl/.D  e 
e n n.n.  n n n.  n.  V / d_ 


where 


l (nR-r/C)US(o,C)  ♦ (nJ-^Vlo.D)  * (n°-»JC>US(J,C)  * (n%j“)US(j,D|| 

L j 

p[n?|»^(T0-R0)]  • PCn”|x“<T0-R0>3  • PC##  VB0)]  ' ■’O 'c1  VE0)]  ^ 


• pip  , 1 D ' D, 

/ ^ ^ \ (n.  +n  ) 

US ( o , C ) = l ^ ♦ Dj  ° (^(t^+^-e)  • * 3 ° <V  | kX  li>(e+k2) 


{ 


d (nC-n’.C) 

l ♦ 3 3 (k3}  ‘ 

k3=° 


(nC-n'C-l)  c 

, » 0 (dR-,0(t1)-k1-k2-k3) 


D 1 D> 


(n"-n.“)  (na"n0D)  D All 

3 ^ (d_k3)  • $ (d0"k1_d+k3)J  } / ’ 


/ 'Cx  'C. 
(n . +n  ) 


, 'EL  'I\ 

(n.  +n  ) 


vn . Tii  > i a \ 

US(a,D)  = m 3 (z0(ti)+kre)  * * (V 


[ $(e+k  ) 


{ 


d (iK-n’.C) 

I <l>  3 3 (k3) 

V° 


(Vn0C)  c , 

> (d0-zo(ti)"krk3+e) 


. D 'D. 
(n.-n.  ) 


, 0 'D  n 
(n  -n_  -1) 


3 3 (d-kj  • 4 ° ° <Vkrk2+e-d+k3) 
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j (»'.C«'C)  (n:D«'D)  \ a 

US(j,C)  = l { Ip  3 (Zo(t1)+k1-e)  ’ <p  3 0^)  < £ lP(e+k2) 


C / C ’C  , . , C ’C. 

(n.-n.  -1)  (n  -n  ) n 

•(  l 4>  33  <k3)4>  ° (do“zo(tl)'krk2'k3) 


, D ' Dv  , D 'D, 

(n.-n.  ) (n  -n„  ) 


• *'  3 3 '(d-k2-k3)$  a ° (d°-k1-d+k2+k3) 


/ 'C  'C.  , 1 D . 'D.  \ , 

(n.  +n  ) (n.  +n  ) d 

3 0 (z„(t, )+k, -e)  * $ 3 (k. ) ( I 


US(  j ,D)  = M $ 3 0 (zQ(t1)+k1 


!>  { l ^e+k2) 

k2=0 

V. 


, C 'Cx  , C 'C, 

(n.-n.  ) (n  -n  ) r 

k 3 3 rv.  w a 0 <aC 


* ( l <f>  3 3 (kgH 


(d0-zo(ti)-krk3+e) 


, D 'D  ,x  , D 'D. 

(n.-n.  -1)  ( n — n ) n 

|>  3 3 (d-k3)4>  (do'kl~k2+e"d+k2+k3) 


The  ranges  for  the  variables  in  the  above  enumeration  are  given  in 
Eqs. (5.25)  and  (5.26). 

Since  the  right  hand  sides  of  Eqs.  (5.19)  and  (5.27)  are  independent 
of  t,  we  obtain  Pr(Bj(*)  = b}  by  substituting  these  equations 
and  Eq.(5.5)  into  Eq.  (4.3). 

The  special  cases  of  complete  recoverability  and  complete 
non -recoverability  can  be  analyzed  in  a manner  similar  to  that 
outlined  in  Section  4.4.1  and  4.4.2,  respectively. 


1 
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5.2. 2 Different  Order  Size  Pic tribut ions  At  the  Ba ser. 

Up  to  now,  all  the  models  studied  for  the  two-echelon  system 
have  assumed  that  order  size  distributions  at  the  bases  are  identical 
In  this  section,  we  allow  for  a more  general  situation  where  the 

order  size  distributions  at  the  bases  are  different. 

At  the  depot,  compounding  distributions  for  the  demand  processes 
(Dg( t) , t>0} , {Dp(t),  t>0}  and  (D^t),  t>0}  are  given  by  Eqs.  (5.2), 
(5.3)  and  (5.4),  respectively.  Also,  Eqs.  (5.5),  (5.6)  and  (5.7) 
provide  the  stationary  distributions  for  the  process  {ZQ(t),  t>0}, 
{QQ(t),  t>0}  and  {By(t),  t>0),  respectively. 

To  obtain  the  probability  distribution  of  the  process 
(B^(t),  t>0} , we  consider  the  two  cases  A and  B in  the  same  context 
as  we  did  in  Section  5.2.1.  In  order  to  evaluate  Eqs.  (4.6)  and 
(4.7)  representing  PrlB^(t)  = b)  for  the  cases  A and  B,  we  obtain 
Pr{U.(t)  = s.+b}„  and  Pr{U.(t)  = s.+b}D  in  Sections  5.2. 2.1  and 

3 : a i 3 B 

5. 2. 2. 2,  respectively.  To  derive  these,  we  use  the  schemes 
developed  in  Sections  5. 2. 1.1  and  5. 2.1. 2.  Also,  we  shall  use  the 
same  notation  as  used  in  these  Sections. 

5.2. 2.1  Case  (A):  Pr{U.(t)  = s.+b}. 

] 3 A 

To  obtain  Pr{U.(t)  = s.+b}  , we  find  the  probability  distributions 

3 3 A 

1 2 

for  its  two  components  U^(t)  and  U^(t)  as  described  in  Section 
5. 2. 2.1.  We  derive  these  by  modifying  the  related  expressions  of 
Section  5. 2.1.1,  for  the  present  case  of  different  order  size 
distributions  at  the  bases. 

Substituting  $ = ^ in  Eq.  (5.8),  we  get 


3.09 


(5.23)  Pr{U^(t ) = s.+b-d}  = CP[s  .+b-d  ! A®R . + ( A^+A1?)!  . , $.]. 

3 3 3 *33  3 3 3 3 


Furthermore,  we  see  that  Eqs.  (5.9  - 5.11)  describing  the 
conditional  probability  distribution  of  lh(t)  hold  for  this  case 
with  4>  = <{> . . Following  the  arguments  used  in  obtaining  the 

i i 

expressions  for  Pr{EX  = 0(  >0);  I=i;  N.(t  ,t  ) = n. ; 

3 ^ ^ 3 

» » 

N (t0,t„)  = n RA}  in  Eqs.  (5.12  - 5.13)  we  have  the  following, 
a z.  3 o' 

» 1 

(i)  0<n.+n  < n.  + n ; (n.  +n  > 1) 

-3  a 3 o 3 q ~ 


Here,  = z^)  - ^(t^)  - e) 

3 a 


z0(tl,-d0(t't2)-e  f.  , , Mn\  c 

i,  'A  ^ 

kF"i 


^tp-e-kpj 


(n.+no)  c 

>0  (z0(tl)"d0(tl*t2)"e)* 


Then 


(5.29)  Pr(EX  = e;  I=i;  ^(tj.tg)  = n^;  = no|RA} 


(n 


..  (Of?) 

,-n . ) \n  . / ' na 

L 1 ° 


1 1 

■j  r 

(n.+n  )-(n.+n  ) /n.+n 


3 a ]o 


3 o 


\ n .+n 
\ 3 O' 


f . i » . 

I (n.+n  ) r 

*^0  (z0(tl)“d0(tl,t2)'e) 


i ♦;(£)> . 


£>e 


for  e = 0,1, . . . ,z0(t1)-dQ(t1,t2)-(nj+no) , 
» 

0 < n . < n. , i = j,o. 

— i—i  J ’ 
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)_d0(tl’t2 V Z0(tl)_d0(tl’t2)  (n.-fn 


112 


5. 2. 2. 2 Case  B:  Pr(U.(t)  = s..+b}B 


Here  again,  the  probability  distributions  of  LL(t)  and 


U.(t)  for  casu  B ar«  derived  using  the  notation  and  the  method 


described  in  Section  5. 2. 1.2. 


Substituting  for  i)>  into  Eq.  (5.20),  we  have 


Pr{U^ (t)  = s^+h-d)  = CF[Sj+b-d|x®R_.  + A°(t_.  + RQ),  fj]- 


C d 

Similarly,  substituting  $0  and  4>  q in  Eq.  (5.21),  we  get 


, C C, 
(n.+n  ) 


(n  ,+ii  ) 


\ II  .r  ii  J r\  n n 

Rr{RB2|RBl}  = Gq  3 ° (d^)  ’ <!>0  3 (dQ). 


Also,  Eqs.  (5.22  - 5.24)  hold  for  this  case  with  $(•)  = 
Following  the  steps  used  in  deriving  Eq.  (5.25)  we  have 


PriEX  = e;  (f/C(t1,t2)  = n/3;  Np^^.tg)  = n^;  i=j,a)i  I-kC|RB},  k- j , a 


, C *C, 
(nk  " nk  ) 


c D C D,  / ' C 'D  'C  'D. 
(nj  + nj  +n0  + no)  " (nj  + "j  + "a  + no  } 


‘ PS 


f r C(n’.Ctn,C) 

■{  i !»S  3 0 <vv  * h - e) 

lki  t 


C C 'C  ’C 


/ U 4/  W 4/  . v 

(n.+n^-n^  -n^  -1)  c 


( l 1 ° ’ ° <d0-ZC<tl>-Vk2,i 


^ ' i'-A 


113 


. , D(n'D+n'D)  D (nD+nD-n'.D-n'D)  D . n 

(*0  ] 0 (k!}  • 'i>0  3 y 3 a (d0-kx)r 


/ C (n^+nC).  C.  D(n?+nD)  D.l 

\*o  1 0 (do)<po  3 0 (V/’ 


k = j,a 


where  the  ranges  for  k^  and  k^  are  the  same  as  given  in  Eq.  (5.25). 

. 'C  'C  'D  'D 

Similarly,  we  can  obtain  Pr{EX  = e;  (N^  (t^t^  = ni  ’ Ni  ^tl,t2^  = ni  ’ 

i = j,o);  I = kDjRB},  k = j,a.  Then  proceeding  in  a manner  similar 

to  that  in  deriving  Eq.  (5.27),  we  can  obtain  Pr{U.(t)-s ,+b} . . This 

1 1 B 

involves  a substantial  amount  of  -numeration . The  detailed  expression 
is  emitted  here. 

The  special  cases  of  complete  recoverability  and  complete 
non-recoverability  can  be  dealt  in  the  same  way  as  mentioned  in 
Section  4.4.  We  emphasize  that  approach  developed  in  this  section 

can  also  be  used  for  a more  general  case  where  the  three  demand  pro- 

B C D 

cesses  {Lu(t),  t>0} , (Dj(t),  t>0}  and  (D^Ct),  t_>0}  at  base 

j ( j = 1, 2, . . . , J)  are  independent  compound  Poisson  processes  with 

different  compounding  distributions. 


5.3  The  Unit  Model 

In  this  model,  upon  arrival  of  a requisition  at  a base,  each 
unit  in  the  failed  batch  is  inspected  independently.  There  are 
three  possible  outcomes  of  each  inspection.  At  base  j,  a unit 
is  base  repairable  with  probability  r_.,  is  depot  repairable  with 
probability  (1-r  ,)p,  or  with  probability  (br.Kl-p)  is 
condemned. 

We  shall  use  the  following  nomenclature  in  addition  to 


those  introduced  in  Sections  4.1  and  5.1. 


Dj  = the  nu.T±>er  of  units  in  a batch  failed  at  base  j 

(j=l,2,...,J). 

B 

0^  = the  number  of  base  repairable  units  in  a batch  failed 

at  base  j (j=l,2,...,J). 

C 

Dj  = the  number  of  condemned  units  in  a batch  failed  at 
base  j ( j=l,2, . . . ,J). 

Dj  = the  number  of  depot  repairable  units  in  a Latch  failed 
at  base  j (j=l,2,...,J). 

0 

D..  = the  number  of  units  demanded  from  depot  upon  a failure 

J 

at  base  j ( j=l,2, . . . ,J). 

0 C D 

Obviously,  D.  = D.  + D.. 

3 3 3 

As  an  implication  of  inspection  under  the  unit  model,  we  have 


(5.32)  r r(D?  = d®;  D°  = d?;  D?  = d?|  D.  = d.} 
3 3 3 3 3 31  3 3 


d t d d d 

-R  - Cr.]  j [(l-r.)(l-p)]  j[(l-r.)p]  j , 

d.MTld.:  3 3 3 

3 3 3 BCD  BCD 

for  dT,d.,d.>0,  and  d.=d7+d.+d.; 
3 3 3-  3 3 3 3 


j otherwise . 


It  can  be  easily  shown  from  Eq.  (5.32)  that  the  conditional 

g g 

probabilities  of  D.,  DT  and  D.  given  D.  = d.  are  B[d.  ■ r.], 
3 3 3 3 3 3 3 


B[dji  (l-r.,)(l-p)]  and  B[d^;  (l-r^)c],  respectively;  where 
BCn,  p3  denotes  the  binomial  probability  distribution  with 


parameters  n and  p (n=number  of  trials,  p = probability  of 


success 


0 < p < 1).  Therefore 
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<»  d . k d .-k 


B(k)  = Pr{DB  = k}  = l (k])tr.]  Cl-r.]  *A&.) 

3 3 d.=k  K J J 

1 


d.-k 


^(k) 

1 


?-ADC.=k}  l ( j)[(l-r.)(l-p)]  [r.+(l-r  )p]  : ^(dJ, 

3 d .=k  k ] 

3 


and 


<»  , d . 


d.-k 


4.B(k)  = Pr{DB  = k>  = l ( -0  C(l-r.)p]  [l-(l-r.)p]  1 ^ dj ) 


d.-k  k 
3 


Furthermore,  let  cji^k^jk^)  = Pr{D^  = k^;  - k^}.  Then 

it  follows  from  Eq.  (5.32)  that 


(5.33)  «fr.(k  ,k  ) = l 

J X * A -V 


d.: 

JL 


"j=kl+k2  kiIk2:  ^dj-kl_k2^: 


k 

• [(l-r.)(l-p)]  [d-fj)p] 


ko  dn"kl"k2 

2r-  1 2 4>.(d.), 


k1>k2  1 °> 


and,  therefore. 


(5.34)  4>?(k)  = Pr{D^  + = k)  = l <f>,(k  , k-k  ) 

3 3 3 k =0  J 


00  d 


- «.x  . d.-k 

y ( ^(l-r.)  r.  ^ 4> . (d . ) , k>_0. 

-u  v l.  ' 3 3 3 3 


d.=k  k 
3 


11G 


Us  note  that  in  this  case  D?  is  a non-negative  variable,  whereas 

in  the  case  of  batch  model  it  is  strictly  positive.  In  our  analysis, 

we  shall  also  include  the  base  requisitions  levied  on  the  depot  for 

which  the  order  size  is  zero.  Consequently,  the  depot  demand 

process  {D  (t),  t>0}  is  a compound  Poisson  process  with  parameter 
J i J o 

Aq  = £ A.  and  compounding  distribution  ^(k)  = ” 1 ^cr 


j=l 

n>l,  let 


0 j=l 


(n), 

1 


(kl»k2)  = J . ^ ^jn"1)(kl"j1»  k2"j2)*j(jl,j2)*  jl’j2-°’ 


j^o  j2=o 


vr’ncre  (f^0^(0,0)  = 1 and  $P(k..  ,k0)  = 0 for  k^,k2  i !•  Then  the 

r p 

process  {D~(t),  D (t),  t_>0 } is  a compound  Poisson  process  with 
‘ 0 0 ' J 

parameter  Aq  and  compounding  distribution  = — £ Aj^Ck^jk^) 

that  is. 


(5.35)  Pr{D0(t)  = dQ;  D (t)  = dQ)  = £ (d  ,d  ) 

n=0 


-At 

0 f , . .n 

(n),X  .D.  ' ° ’ 


dc-  d0  i °- 


C D 

In  addition,  the  processes  {DQ(t),  t>0}  and  {DQ(t),  t>0}  are 
compound  Poisson  processes  with  common  parameter  A and  compounding 

J J 

distributions  k ) = £ A.^(k),  k>0,  and  <Ak)  = t—  I A.  d>  ? ( k ) , 

0 A0  j=l  33  0 A0  j=l  3 3 

k>0,  respectively.  The  depot  can  now  be  analyzed  as  a single 

location  system  where  recoverable  and  non-recoverable  demand  processes 

are  dependent  compound  Poisson  processes.  Therefore,  the  results 

derived  in  Section  3.4.2  apply. 

C 

From  Theorem  3.5,  upon  substituting  C * ) for  p(>),  we 


have 


r 
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(5.36)  lim  Pr{Z  (t)  = kjz(O) 
t;-x» 


g(S-k) 


1+G( 5-3-1 ) 
1 


l+G(S-s-l) 

V. 


) k-s^ 1 j ... *5— 1 j 


k=S, 


k-1 


where  g(l)  = <^>0(l)/  Cl-4>0C  0)3,  g(k)  = [$Q(k)  + V 4>0(q)gO<-q)  ]/ 

[1-^(0)],  k=2, 3, . . . , and  G(k)  l g(l). 

1=1 

Also,  it  easily  follows  from  the  previous  discussions  that 


(5.37)  lim  Pr{Q0(t)  = q0|QQ(0)  = 0}  = CPCq^X^,  ♦£],  qgO. 


t-XX> 


The  transient  distribution  of  the  process  {Bg(t),  t>0}  can  be 
obtained  from  Eq.  (3.43).  As  mentioned  in  Section  3.3.3,  the  derivation 
of  stationary  distribution  is  computationally  complex.  The  results, 
however,  can  be  obtained  using  Laplace  transforms.  The  explicit 
derivation  is  omitted  here. 

To  obtain  the  stationary  distribution  of  the  process  (B^(t),  t>0), 

we  assume  that  the  base  repair  time  does  not  exceed  the  base  procurement 
lead  time-,  that  is,  R.  > t..  The  analogous  results  for  the  case 
j-w  > Tj  can  be  derived  using  the  approach  described  here. 

Here  again,  we  compute  Pr{lK(t)  » s^+b}  for  the  two  cases 
A and  B described  in  Section  4.2. 


5.3.1  Case  A:  Pr(U.(t)  = s.+b} 

3 3 

To  obtain  Pr{U.(t)  = s.+b}  , we  find  the  probability 

3 3 A 

distributions  for  U^(t)  and  U?(t),  the  two  components  of  LR(t) 
in  case  A.  Under  the  assumption  <_  it  is  clear  that 


A i »•  ■ 


it  is  clear  that 
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!T(t)  = DC(tQ,  t-R.)  + D°U 3,  t-R.)  + D (t-R.,  t) 

3 30  j J jj 

= D°(t_,  t-R  ) + D (t-R  , t). 

D J 3 J J 


Because  the  arrival  process  is  Poisson,  Dj(t3>  t-R..) 
D.(t-R.,  t)  are  independent.  Therefore 


To  obtain  Pr(U?(t)  = d },  we  proceed  as  in  section  5.2. 2.1. 
that  D?  is  now  a non-negative  variable.  Then  from  Eq.  (5.31), 
follows  that 


Note 

it 
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where 


X = > X.,  X = y X.  and  $ = -r—  £ A.4>.. 
3=1  J 3=1  J o i=l  J 


probability  distribution  of  IK(t)  is  obtained  from  Eq.  (5.38). 

5.3.2  Case  (B):  Pr{U.(t)  = s.ib}n 

] 1 B 

Here  again,  we  obtain  the  probability  distributions  for 
1 2 

Uj(t)  and  IL(t),  the  two  components  of  IK(t)  in  case  3. 

By  definition,  as  given  in  Section  4.2 


Hence , 


U^(t)  = DC(t  , t-R.)  + D°(t_,  t-R.)  + DB(t-R. , t). 
1 3 2 3 3 2 3 33 


s .+b-d 


(5.40)  Pr{U^(t)  = s.+b-d}  = ^ 7 CP[s  .+b-d-k  I A.,  (t  .+R.-R. ) , <f>?] 

3 3 v-n  3 3 3 0 3 3 


3P[k ! X j R j * * .]• 


As  indicated  earlier,  in  the  unit  model  each  requisition  at 

the  depot  is  associated  with  a demand  for  replacement  for  some 

depot  repairable  units  and  some  condemned  units.  To  obtain  the 

2 

probability  distribution  of  U_.(t),  we  follow  the  steps  described 
in  Section  5. 2. 2. 2.  Let 


RM  = (HjUj.tp  = n.i  N,( VV  * V W 
RB2  = = dC0;  Dj(t1>t2)  = dj) 


W>- 


R3  = KB1  u RB2 


Then  similar  to  Eqs.  (5.9),  (5.10)  and  (5.11),  we  have 
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(5.41)  Pr{U?(t)  = d|RB;  EX  = 0;  N^t^tg)  = nr,  = n^} 


.0  (n.-n.),,v 
L 1 J (d). 


(5.42)  Pr{U?(t)  = d|KB;  EX=  e > 0;  ^(t^tg)  = n.;  N^t^tg)  no*  °} 


0 (n.-n.),,. 
h D J (d) 


and 


(5.43)  Pr{U j (t) 


= d|KB;  EX  = e > 0;  N^^.tg)  = n^;  ^(t^tg)  =nQ;  I=j> 


v .0 ,,  x .O(n.-n.-l),  , * \ 

= } d>.(k+e)  4> . j J (d-k). 

L I 1 

J J 


k>0 


Following  the  arguments  used  in  deriving  Eqs.  (5.25)  and  (5.2b),  we  obtain 


(5.44)  Pr{EX  = e;  lUt^tg)  = ; N/VV  = V I = k|RB),  k = j,o 

= Pr{(Y1+Y2+...+YnO  + (Y1+Y2+...+YnO  = z0(tl)+do"e’ 

» » 

Y . > e;  out  of  first  (n.+n  ) requistions  at  the 

V1  ’ 3 

depot,  are  from  source  j;  the  next  demand  is 

from  source  k | D^t^.tg)  = dQ;  D^t^tg)  - dQ,  RBl), 
k = j,  o 


nk'nk 


r » 

(n.+n  ) 


(n 


l v3  ’ Mr  *o(ti)t4;— di> 
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> I I ♦k(«*vd2>  V 


'2  2 


■N 


• ♦ ^3+n  "i  n°  1)(do"dl"e"k2+d2’dl+e~d2_Z0^tl^/'/ 


(n.+n  ),  ,C  ,D>. 

*o  3 ° (do’V’ 


where. 


i E 

nk  <_  n^(  >^1),  k = j,a;  and  e = 0,1, . . . ,zQ(t1)  + dc 


From  the  property  of  <pn(a,b),  we  see  that  the  ranges  for  the 
indicics  in  the  summations  must  satisfy 

0 < dx  < zQ(t1)  +dQ-e,  lik2<dQ-d1-e  and  0 1 d2  < ^ + e - zQ<t 
Thus  from  Eqs.  (5.41  - 5.44)  we  get 

s .+b 

(5.45)  Pr{U . (t ) = s.+b}  = t ^ Pr{U^(t)  = s +b-d} 

J J D r j j 


do>zo(ti)  d=0 


r 

r 

(y 

l 

-v 

t^3 

l 

l 

!k/ 

i 

t 

/ n . 

+ n \ 

1 n . 

n ' 

In . 

n 

: 

a\ 

3 

a 

a 

■ 

l n! 

+ n'  / 

l l 

jD  „ e 

do=° 


1 O' 


f , . l 

(n  -n  ) • US (a)  + (n.-n.)  • US(j)/  ( 
^ a o 33  Jj 


PCnj 1x3(t0-r0)3  • PCno|Xo(t0-R0)]j  , 
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where 


US(o)  = l 4>onj+no)(d1,  z0(ti)  +dQ-e-d1)  l l 4>°(e+k2-d2,  d2) 


k2  d2 


^ o (ni+nl),  . . . o(VV1),,c 

J l ^ 3 3(d-k3’k3)^0 


>3-°  3 


(d0“dre_k2+d2-d+k3’ 


dl+e-d2_Z0(tl)_k3) 


and 


OS(3)  - l ♦J-K'm,.  l d2> 


k2>0  d2 


^ o (n 

, z ♦ 3 

'<k3  3 


. -n.-l)  0(n  -n')  c 

J 1 (d-k2-k3,  k3)  • «“  0 “(^-d,-e+d,- 


0-1  —2  d+k3- 


Thc  ranges  for  the  indicies  in  the  summations  can  be  easily  obtained 

using  the  fact  that  in  an  expressinn  4>n(a,b),  a,  b>0  for  n>l, 
and  a = b = 0 for  n = 0 . 

Thus  Pr(B.(*)  = b}  can  now  be  obtained  from  Eqs.  (5.36),  (5.39) 


and  (5.45). 


CHAPTER  V 


A COMPARISON  OF  TWO-ECHELON  MODELS 


In  this  chapter  we  compare  the  two-echelon  inventory  models  studied 
in  the  previous  two  chapters  to  Sherbrooke's  METRIC  (Multi-Echelon 
Technique  for  Recoverable  Item  Control)  model  [17].  In  Section  6.1., 
we  briefly  describe  the  METRIC  model.  In  Section  6.2,  a general 
comparison  is  presented,  while  a computational  comparison  is  made 
in  Section  6.3. 


6.1  The  METRIC  Model 

Almost  a decade  ago  Sherbrooke  [17]  developed  the  well-known 
METRIC  model  for  a two-echelon  system  similar  to  that  described  in 
section  1.1.2.  A major  purpose  of  the  model  is  to  determine  optimal 
base  and  depot  stock  levels  that  minimize  the  expected  number  of 
total  system  backorders  subject  to  a constraint  on  system  investment 
or  system  performance. 

Demand  at  base  j is  assumed  to  be  represented  by  a compound 
Poisson  process  with  parameter  A_.  ( j=l,2, . . . ,J).  Upon  arrival  of  a 
requisition,  one  or  several  units  are  demanded  for  replacement  and 
a like  number  of  failed  units  are  turned  in  for  repair  at  the  base. 

A batch  of  failed  units  at  base  j is  repaired  at  the  base  with 
probability  r.j  and  is  shipped  to  the  d^pot  for  repair  with  probability 
(1-r^).  Thus  there  are  no  condemnations;  that  is,  the  system  is  by 
definition  conservative.  Bases  use  an  (s-1,  s)  policy  for 
procurement  of  serviceable  units  from  the  depot.  Repair  time  at  location 
j is  a random  variable  with  finite  and  known  mean  R^  ( j=0,l, . . . , J). 
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It  is  assumed  that  the  repair  time  is  the  same  for  each  unit  in  a 
failed  batch.  The  order-and-ship  time  at  base  j is  also  a random 
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variable  with  mean  t j=l, 2, . . . , J) . Furthermore,  it  is  assumed  that 

the  order  size  distribution  is  the  same  at  all  bases.  Let  4>  he 

the  mean  order  size  per  requisition  at  the  bases.  Thun  the  mean  depot 
J 

demand  rate  is  V A.(l-r.)<j>  = where  A.  = T (l-r.)A.. 

j=l  3 3 0-  ° j=l  : 3 

In  the  METRIC  model,  the  stationary  distribution  of  the  number 
of  backorders  at  a base  is  derived  from  Eq.  (12)  of  Feeney  and  Sherbrooke 
[5].  Let  S be  the  spare  stock  (inventory  on  hand  + on  order  + in 
repair  - backorders)  for  an  item  where  demands  are  compound  Poisson 
with  parameter  A,  and  resupply  time  is  a random  variable  with  mean 
T.  Also,  assume  that  the  resupply  time  is  the  same  for  ail  units 
demanded  by  a requisition.  Then  in  the  case  where  backlogging  is 
allowed,  the  number  of  units  in  resupply  has  a Poisson  distribution 
with  mean  AT. 

Let  Sq  be  the  spare  stock  at  the  depot  and  T_.(SQ)  be  the 
average  response  time  to  a demand  from  base  j . Sherbrooke  shows 
that  Tj(SQ)  can  be  expressed  as 

(6.1)  T.(SQ)  = rjRj  + (1-rj)(Tj  + &(-So)  ' V 


where  P[*|  m]  denotes  the  Fbisson  distribution  with  mean  m. 
tiie  METRIC  model,  the  stationary  distribution  of  the  number  of 
backorders  at  base  j is  given  by 


Thus  in 
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(6.2) 


Pr{B^(*)  = b}  = P[Sj+b| AjTj(SQ)] 


PCs.+b|XBR.+A?T.+A?  6(S  ) "R-] 
]]]]]]  0 0 

b=~s_.  0 j i ^ ^ 


where  AB  = r.A.  and  A?  = (l-r.)A.. 
3 3 3 3 3 3 


6.2  General  Comparison 

We  now  compare  some  of  the  features  of  the  METRIC  model  and  the 
two-echelon  models  studied  in  chapter  IV  and  V.  Our  models  are  more 
general  than  METRIC  in  that  they  permit  non-recoverability  and  positive 
condemnation  rates.  Our  analysis  includes  both  batch  and  unit  models 
for  inspection  of  failed  units,  whereas  METRIC  considers  only  the 
batch  model.  We  also  examined  the  case  where  order  size  distributions 
are  different  at  the  base.  The  METRIC  model  is  confined  to  the  case 
where  the  order  size  distribution  is  the  same  at  all  bases. 

The  METRIC  model,  on  the  other  hand,  is  more  general  in  that  it 
allows  random  repair  and  order-and-ship  times.  Furthermore,  METRIC 
provides  simple  but  approximate  expressions  for  the  stationary 
distribution  of  the  number  of  backorders  outstanding  at  a location  at 
any  point  in  time.  In  our  analysis,  the  expressions  are  relatively 
more  complex  and  a large  amount  of  enumeration  is  required: to. compute 
the  probability  distributions. 
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6.3  A Computat io nal  Comparison 

In  this  section,  we  assess  the  degree  to  which  the  METRIC 
model  could  serve  as  a useful  approximation  to  our  model  which  we 

shall  refer  to  as  the  EXACT  model.  Also,  we  compare  the  computational 

complexities  of  the  two  models  in  terms  of  their  execution  times  on 
a digital  computer. 

We  consider  the  case  of  unit  demands  at  the  bases.  For  the 
EXACT  model,  the  conservative  system  is  discussed  in  section  4.4.1  and 
the  stationary  distribution  of  the  number  of  backorders  at  the  bases 

is  given  by  Eq.  (4.25).  For  the  METRIC  model  it  is  computed  from 
Eqs.  (6.1)  and  (6.2)  upon  substituxing  ^ - 1 into  Eq.  (6.1).  Furthermore, 
the  term  5(S  )*R  is  simplified  as  a finite  sum  and  the  resulting 
expression  is  given  by 


(6.3)  5(SQ)-R0 


U 


x=0 


P [x 


'W} 


S 

SQ  {1-  l P[x|X0R0]}J 
x=0 


For  the  purpose  of  computational  comparison,  we  considered  several 
numerical  examples  covering  a wide  range  of  the  system  parameters  - 
demand  rates  at  bases,  depot  and  base  repair  times,  depot  and  base 
spare  stock  levels  and  probability  of  a unit  being  repairable  at  a 
base.  When  the  comparison  was  made,  the  difference  between  the  METRIC 
model  and  the  EXACT  model  was  found  to  have  the  same  pattern  in  all  cases. 

The  extent  of  discrepancy  between  the  two  models,  however,  did  vary  from 
one  example  to  the  other  depending  on  the  numerical  values  of  their 
parameters.  To  describe  the  comparison  between  the  two  models,  the  following 
example  was  selected  as  representative  of  all  situations  considered. 
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Example 


Number  of  bases  = 5 


Base 

(j) 

Arrival  rate  (per  day) 

<V 

Gi’der-snd -ship  time  (days; 

<V 

P 

Repair  time  (days) 

(R.) 

J 

1 

0.07098 

12 

6 

2 

0.14766 

12 

6 

3 

0.01497 

12 

6 

4 

0.04591 

12 

6 

5 

0.11019 

12 

6 

Repair  time  at  the  aepot  (RQ)  = ^6  days 


The  stationary  distribution  for  the  number  of  units  in  resupply  at 
the  bases  is  computed  for  the  two  models.  Tables  I,  II  and  III  show 
such  distributions  for  r^  = 0.10,  i-v  = 0.50  and  r^  = 0.90,  respectively. 
It  is  assumed  that  r^  is  the  same  for  all  bases.  The  corresponding 
arrival  rate  at  the  depot  (A  ^)  is  stated  at  the  beginning  of  each  table. 
The  depot  spare  stock  level  (SQ)  for  each  case  is  chosen  to  be 


approximately  equal  to  AqRq,  the  average  number  of  units  in  resupply 
at  the  depot.  Similarly  the  value  of  s ^ , the  spare  stock  level  at 
base  j , is  chosen  to  be  the  least  integer  greater  than  or  equal  to 
A^Cd-r^Jr.  + r^R,.].  The  letter  x refers  to  the  number  of  units  in 
resupply.  The  expected  number  of  backorders  (EB0)  is  computed  in  the 
last  row  of  each  table.  The  execution  time  for  the  two  models  is 
stated  below  each  table.  All  programs  were  written  in  FORTRAN  IV  and  run 
on  Cornell's  IBM  370/168  computer. 


The  data  displayed  in  the  tables  indicated  that  the  METRIC  model 


The  Probability  Distribution  ot  the  Humber  of  Units  in  Resupply  at  the  Bases 
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underestimates  the  stationary  probability  of  zero  unit  in  resupply  (s_. 

units  on  hand  at  base  j).  Also,  in  the  EXACT  model,  th<  stationary  distri- 
bution of  the  number  of  units  in  resupply  has  louder  tails.  Furthermore , 

we  wish  to  investigate  how  well  the  METRIC  model  approximates  the 
stationary  distribution  when  the  depot  snare  stock  level  (S^)  changes. 

A similar  investigation  will  be  made  when  r_. , the  probability  that  a 
failed  unit  is  repairable  at  base  j,  changes.  For  this  purpose,  we 
arbitrarily  examine  base  2. 

Figures  6.1,  6.2  and  6.3  show  the  probability  distribution  of  the 
number  of  units  in  resupply  at  base  2 when  S.  = 15,  19  and  23, 

J 

respectively.  The  value  of  r is  fixed  at  0.10.  It  is  clear  from 
these  figures  that  the  discrepancy  between  the  METRIC  model  and  the  EXACT 


model  decreases  as  the  value  of  SQ  increases.  This  can  be  explained 


as  follows.  As  S increases,  in  Eq.  (4.25),  the  contribution  of  the 

f /dO(t2,t3)“So\  d (t  t )-d-S 

^ ' ) • CA?/Ajd[l-A?/Aj  0 2’  3 0 


term 


l - 


d0^t2*  t3'*=‘:>0+d 


\ 


j oJ 


j 0‘ 


PCdo(t2,t  ,)  | AqRq],  d >_  0 decreases  and  Pr{B_.(*)  = b}  is  dominated 


B 0 

by  the  Poisson  term  P[s.+b|A.R.  + A.t.]  • 

3 3 3 13 


. Jl  ...  PW0(t2’t3)l10R0:li- 

-WV'0 
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For  the  METRIC  model,  as  Sq  increases  the  term  6(SQ)  * RQ  representing 

the  average  delay  per  depot  demand  decreases  and  consequently  T.. (SQ),  the 

average  response  time,  decree '•rs.  In  Eq.  (6.2),  Pr(B^(*)  = b)  is 

dominated  by  the  Poisson  term  P[Sj+b|A^Rj  + A^t^].  Thus  we  conclude 

that  for  the  values  of  SQ  sufficiently  large  (>  AQR0),  the  Poisson 

approximation  is  close  to  the  results  given  by  EXACT  model.  As  ■+  », 

B 0 

T.(s-)  = A.P..  + A.t..  Thus,  the  two  models  are  the  same  and 
3 0 3 3 3 3 

Pr { B . ( * ) - b}  = P[s.+b|A?R.  + A?t.}. 

3 3 3 3 3 3 


probability 


METRIC 


number  of  units  in  resupply 


Lgure  6.4:  The  probubility  distribution  of  the  number  of  units 

in  resupply  at  base  2 (r  = 0.5,  S - 11). 


137 


! 


figure  6. 


EXACT 

METRIC 


2 3 4 5 6 78  9 10 


number  of  units  in  resupply 

The  probability  distribution  of  the  number  of  uni  ‘.s 
in  resupply  at  base  2 ( r = 0.90,  SQ  = 2). 
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We  now  investigate  how  well  the  METRIC  model  approximates  the  EXACT 
model  when  r^,  changes.  Figures  6.2,  6.4  and  6.5  show  the  probability 
distribution  of  the  number  of  units  in  resupply  at  base  2 for  the  case 
when  r2  = 0.10,  r2  = 0.  50  andr2  = 0.90,  respectively.  The  depot'spare  stock 
level  for  these  cases  is  fixed  at  19,11  and  2,  respectively.  The  stock  level 
is  approximately  equal  to  the  corresponding  value  of  A^R^  each  case. 

It  is  clear  from  Figures  6.2,  6.4  and  6.5  that  the  discrepancy  between 
the  METRIC  and  the  EXACT  model  decreases  as  r^  increases.  This  can 
be  explained  in  a way  similar  to  the  previous  case  of  comparing 

g 

the  two  models  with  respect  to  S„.  When  r.  increases  A.  increases  and 

0 d 1 

0 

A..  decreases.  Consequently,  the  proportion  of  supply  received  at  the 
base  from  its  repair  process  increases.  Thus,  in  both  Eqs.  (4.25)  and  (6.2), 

g 

Pr{B_.  (*)  = b)  is  dominated  by  the  Poisson  term  P[s^+b|AjR..].  As 

r^  -*•  1,  both  Eqs.  (4.25)  and  (6.2)  reduce  to  P[s_.+b  | A_.R_.  ] and  the  two 

models  are  the  same. 

As  mentioned  by  Sherbrooke  Cl7],  the  expected  number  of  backorders 
at  the  bases  (EBO)  is  an  important  measure  of  system  performance.  We 
must  make  an  assessment  of  how  wdl  the  METRIC  model  approximates  the 
computation  of  EBO.  From  tables  I, II  and  III  it  is  clear  that  the  METRIC 
model  always  gives  the  values  of  EBO  which  are  less  than  or  equal  to  the 
corresponding  values  given  by  the  EXACT  model.  Furthermore,  it  is  clear 
from  Figures  (6.1-6. 5)  that  the  METRIC  curve  is  always  below  the  EXACT 
curve  near  the  tail  of  the  distribution.  In  other  words,  the  METRIC 
model  always  underestimates  the  probability  for  the  larger  numbers  of 
units  in  resupply  at  a base.  Consequently,  as  the  base  spare  stock  level  is 
increased  the  difference  between  the  values  of  EBO  given  by  the  two  rood  ALs 
increases.  Table  IV  displays  such  a discrepancy  for  base  2 of  the  example. 
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Table  IV 

The  Expected  Number  of  Backorders  at  Base  2 

r = 0.50.  S_  = 11 

NJ 

c 

0 

1.808638 

1.808638 

1 

0.972516 

0.999878 

2 

0.432787 

0.485504 

3 

0.161094 

0.213471 

4 

0.050994 

0.086882 

5 

0.013959 

0.033091 

6 

0.003354 

0.011839 

7 

0.000717 

0.003984 

8 

0.000138 

0.001262 

9 

0.000024 

0.000377 

Finally,  the  EXACT  model  requires  more  execution  time  than  the 
METRIC  model.  This  is  so  because  in  the  EXACT  model,  the  expression 


140 


for  Pr{IK(*)  = b}  is  relatively  complex.  It  involves  the  computations 
of  the  combinatorial  terms  in  addition  to  the  Poisson  terms,  and  a 
large  number  of  calculations  are  required  to  evaluate  Pr{B_.(*)  = b}. 

We  conclude  this  chapter  by  emphasizing  that  although  the  METRIC 
model  is  computationally  simple,  its  use  in  certain  situations  may 
not  be  appropriate.  When  the  depot  spare  stock  level  Sq  is  low  (less 
than  A R , the  mean  demand  during  the  repair  time)  and  the  proportion 
of  repairs  done  at  the  depot  is  high  (r^  is  low),  there  is  a 
pronounced  difference  between  the  two  models.  Thus  in  this  situation 
the  METRIC  results  may  be  misleading.  For  higher  values  of  and 

r..,  however,  the  approximate  results  given  by  the  METRIC  model  are 
fairly  close  to  the  results  given  by  the  EXACT  model.  Furthermore, 
when  an  optimization  algorithm  is  applied  to  determine  the  optimal 
stock  levels  of  several  items  at  different  locations  in  the  system 
under  a budget  constraint,  the  solution  may  be  wrong  if  the  METRIC 
results  are  used.  This  may  happen  because  of  the  above  mentioned 
varying  discrepancy  between  the  two  models  with  respect  to  S(.  and 
v. . Suppose  the  problem  is  to  determine  the  optimal  stock  levels  for 
two  items,  for  one  of  which  the  METRIC  and  the  EXACT  results  are  close 
while  for  the  other  there  is  a considerable  difference  between  the 
results  given  by  the  two  models.  In  other  words,  for  the  second  item 
the  METRIC  curve  is  steeper  than  the  EXACT  curve  towards  the  tail 
of  the  distribution  of  the  number  of  units  in  resupply.  Thus  the 
reduction  in  the  expecteo  number  of  backorders  per  additional  unit 


' JfcrV:  •*■££ 


of  the  second  item  is  larger  for  the  METRIC  model.  Consequently, 
it  will  erroneously  result  in  a larger  appropriation  of  spare  stocks 
for  the  second  item  if  the  METRIC  results  are  used.  A similar  conclusion 
can  be  drawn  for  the  case  when  the  problem  is  to  allocate  optimally 
the  spare  stocks  of  a single  item  at  different  bases  in  the  system. 

Here,  it  easily  follows  from  the  above  arguments  that  for  the  bases 
at  which  the  repair  capability  is  poor  (r^  is  low),  the  METRIC 
model  will  suggest  larger  spare  stock  levels  than  those  given  by  the 


EXACT  model. 


CHAPTER  VII 


CONCLUDING  COMMENTS 

In  this  study,  we  have  presented  an  analysis  for  continuous 
review  models  of  single  location  and  two-echelon  inventory  systems 
for  recoverable  items  with  random  demands. 

For  the  single  location  system,  the  (s,S)  and  (s,nQ)  policies 
have  been  examined.  Under  each  policy,  stationary  distributions 
for  inventory  position,  on-hand  inventory  and  backorders  have  been 
obtained  for  both  independent  and  dependent  demand  processes.  For 
the  diagnosis  of  failed  units  two  appealing  models  - batch  and  unit  - 
were  examined.  Our  analysis  is  more  general  than  any  previous  work 
since  only  independent  demand  processes  and  batch  inspection  were 
considered  before.  We  have  also  extended  our  results  to  the  case 
where  demand  processes  are  compound  Poisson  processes.  In  addition, 
we  have  demonstrated  how  the  special  cases  of  complete  recoverability  and 
complete  non-recoverability  are  obtained  from  our  results  with  a 
simple  change  of  parameters. 

As  mentioned  in  Section  3.1,  the  analysis  of  a single  location 

system,  in  addition  to  offering  the  solution  to  the  system  itself,  is 
of  great  importance  for  the  analysis  of  multi-echelon  systems.  The 
results  that  we  obtained  do  not  offer  a complete  solution  to  the  complex 
problems  facing  inventory  managers  in  such  a system.  Using  the  approach 
and  results  obtained  here  , further  analysis  may  be  pursued  to  evaluate 
cost/benefit  trade-offs  among  the  various  system  parameters  and  to  examine 
alternative  repair  disciplines. 
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For  the  two-echelon  system,  we  have  mainly  concentrated  on  obtaining 
an  exact  expression  for  stationary  distribution  of  the  number  of 
backorders  at  each  base.  The  approach  suggested  by  Kruse  and  Kaplan  [9] 
has  been  extended  to  the  case  of  an  arbitrary  order  size  distribution 
at  the  bases.  Here,  too,  the  batch  and  unit  inspection  models  have  been 
considered.  The  results  have  been  indicated  for  the  two  extreme  cases  of 
complete  recoverability  and  complete  non-recoverabilitv  of  the  item.  For 
the  case  of  complete  recoverability,  we  have  demonstrated  how  the  METRIC 
model  may  produce  a misleading  solution  to  the  problem  of  allocating 
items  in  the  system. 

Our  results  for  the  two-echelon  system  depend  heavily  on  the  assumption 
that  failures  are  generated  by  Poisson  processes,  repair  and  lead  times  are 
constant,  and  partial  backlogging  is  allowed.  Also,  we  have  considered 
and  (s-l,s)  policy  at  the  bases  and  (s,S)  policy  at  the  depot.  We 
have  assumed  that  a first-come,  first-served  policy  is  used  at  all 
locations  to  fill  the  backlogged  demands . We  make  no  claim  that  these 
assumptions  are  always  true.  Nevertheless,  our  approach  should  be 
useful  as  a basis  for  further  analysis. 

The  methodology  developed  in  chapter  IV  and  V may  be  used  to  obtain 
the  stationary  distribution  of  the  number  of  backorders  when  an  (r,Q) 
or  an  (s,S)  policy  is  used  at  the  bases.  It  will  be  interesting  to 
examine  the  situation  where  the  inter-arrival  time  at  each  base  has  an 
Erlang  distribution.  The  Erlang  distribution,  as  mentioned  by  Gross  and 
Harris  ([6],  p.  162-163),  provides  much  more  modeling  flexibiltty  than 
does  the  exponential.  It  can  be  used  to  provide  good  approximation  for 
many  different  inter-arrival  distributions. 
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An  attempt  may  be  made  to  analyze  the  situation  where  the  depot 
does  not  strictly  follow  the  first-come,  first-served  policy  to  fill 
the  base  demands.  On  the  other  hand,  a priority  scheme  is  followed 
at  the  depot  such  that  the  bases  with  low  mean  demand  would  have 
priority  over  the  bases  with  high  mean  demand.  Certainly,  this  will  not 
affect  the  total  expected  number  of  system  backorders,  but  will  result 
in  lower  spare  stock  levels  at  the  low  demand  bases. 

Finally,  our  approach  can  be  directly  used  to  analyze  arborescence 
systems  with  more  than  two  echelons  where  all  the  bases  and  intermediate 
echelons  use  an  (s-l,s)  policy  and  the  depot  uses  an  (s,S)  policy. 

In  addition  to  bases  and  the  depot,  repair  can  be  performed  at  intermediate 


echelons . 


APPENDIX  A 


ARRIVAL  SEQUENCES 

In  this  appendix,  we  derive  some  results  concerning  the  sequence 
of  arrivals  generated  by  p (_>2)  independent  Poisson  processes.  Vie 
first  consider  the  case  when  p = 2.  Let  {N^t),  t ■>  0}  and 
{^(t),  t >_  0}  be  the  two  independent  Poisson  processes  with  parameters 
X and  X2,  respectively.  Also,  let 

N(t)  = total  number  of  arrivals  during  the  interval  (0,t]; 

= the  time  of  the  kt‘'1  arrival  from  process  i,  i = 1,2  and 
k I 1 (Tq  = °>  1 = l»2)i 

and 

W = the  time  of  the  x1-^  arrival  in  the  total  arrival  process 
x 

{N(t ) , t > 0}. 

We  have  the  following  theorem. 

Theorem  1.  Let  n^  and  n^  be  non-negative  integers  such  that 
n^  + n2  > 0.  Also,  suppose  we  are  given  non-negative  integers 
n|  and  n^  such  that 

(a)  nj  <_  nj  + n^  < n^  + n? 

(b)  n|  + n^  > 0. 

Then 
(A.l) 
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th 


Proof : We  know  that  T-1", , the  time  of  n'  arrival  from  process  1, 

nl  1 

has  a gamma  distribution  with  parameters  n|  and  X^  (see 


Ro; 


[14],  p.  16-17);  that-  is-,  the  probability  density  function  of  T is 

n' 

?>iven  by  , n'  n'-l 

_ 1 (X, ) (u) 


(A. 2) 


0 < u < t 


(n^-1) ! 


We  can  vrrite 


(A. 3)  Pr{N2(Tb,)  = n^|N1(t)  = n^  NjCt)  = n2> 

Pr{N2(Ti.)  = n2 ; N^t)  = N2(t)  = n2} 

FrlN^t)  = n^,  N 2 ( t ) = n2> 


The  numerator  of  Eq.  (A. 3)  can  be  evaluated  by  conditioning  on  T 

nl 

Thus  we  have. 


(A. 4)  Pr{N„(T  ,)  = n';  N.(t)  = n. ; N„(t)  = n.} 

^ n^  2 1 12  2 

= Jg  Pr{N2(u)  = n2 ; N^t)  = ^ ; N2(t)  = n2|Tb,  = u}  • f(Tb,,  u) 
= Jq  Pr{N2(u)  = n2;  N^t-u)  = n1~n|;  N2(t-u)  = n2-np«  f(T^,, 


-Xu  n^  -X  (t-u) 
t e (X2u)  " e [ X? (t-u)] * ~ e 

^0  n^!  ’ (n^-n,'  Vi 


n -n  -X0(t-u)  n0-n' 

11  1 ~ - [>2(t-u)] 

(n2~n2)! 


-Xu  nf  n'-l 
e 1 (X1)  x(u)  1 

(nj-1)! 


du 
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"V  “l  'V  . n2 

e X,  e X, 


1 2 

(n'-l):  n’!  (nrnp!  (n2-n > 0 U 


t ni+n2_1  , . nl+n2-nl-n2  . 

(t-u)  du 


-(X  +X  )t  n 
e 1 “ (X1)  1 (X2) 


2 

r‘2 


(n^-1)!  n£!  (n-j-nJ  ) ! (n2~np! 


. t 


n,j+n2  (nj+n^-1)!  (n^n^n^-n^ ) ! 

(n1+n2) ! 


V.'e  know  for  the  denonimator  of  Eq.  (A. 3)  that 


(A. 5)  Pr{N1(t)  = n^-  N2(t)  = n2>  = 


-X,t  n -X  t n 
e * (Xxt)  e (X2t) 


n. 


n. 


Substituting  Eq.  (A.. 4)  and  (A.  5)  into  (A.  3)  and  upon  simplifying  we 
obtain  Eq.  (A.l). 

Q.E.D. 

We  now  have  the  following  tv?o  corollaries. 

Corollary  1.1  Under  tht  assumptions  of  Theorem  1,  let  x,  y be  non- 
negative integers  such  that  y <_  x < n1  + n2  and  x > 0.  Then 


(A. 6)  Pr{N1(W}{)  = y | N^(t ) = n^  N2(t)  = n2> 


Proof: 

We  note 

that  the  event 

N (W  ) = y is  equivalent  to  the  event 

X X 

W 

= x-y. 

The  waiting  time 

W will  coincide  either  with 

X 

T or 

T2  . 

Therefore, 

y 

x-y 

1 1 
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PrtN^W  ) = yjN^t)  = i^;  N2(t)  = n2> 

* PrtN^T*)  = x-ylN^t)  = nx;  N?(t)  = n2} 

+ Pr{N1(T^_  ) = xlN^t)  = ; N2(t)  = n2> 


rp) 

- ry  , x-y  \ .\y/vx-y/  ( 

(•■)(•■) 

' y ' ' x-y 


from  Theorem  1) 


y ' ' x-y 

(nl  + n2) 

\ x ' Q.E.D. 

The  above  corollary  has  been  also  proved  by  Simon  [20]  under 
more  general  conditions. 

Corollary  1.2.  Let  I(s)  denote  the  index  of  the  process  from  which 
the  next  arrival  occurs  after  time  s (0  <_  s < t).  Then  under  the 
assumptions  of  Theorem  1, 

(A. 7)  Pr{I(Wx)  = i|Ni(W}{)  = n!,,  N^t)  = N^t)  = n2> 


n.-n. 
i i 

(n1+n2>  - x 


i = 1,2. 


Proof:  We  prove  for  the  case  when  i = 1.  Given  that  N^(t)  = n^, 
j^(t)  = n2  and  N^CW^)  = n^,  the  event  ICW^)  = 1 means  that  the 
next  arrival  after  n^  arrivals  from  process  1 and  (x-n|)  arrivals 


t 


frot.'.  process  2 is  from  process  1. 


In  other  words,  by  the  time 


(n^+l)*"^  arrival  from  process  1 occurs,  there  are  (x  - n^)  arrivals 
frcm  process  2.  Thus  we  have. 


Pr{I (W  ) = ljJUW  ) = n!;  N.(t)  = n ; N (t)  = n } 

X 1 i.  X 11  1 * * 

Pr{I(W  > = 1;  N3(w J = n^N^t)  = ry,  N2(t)  = y} 
Pr{N1(Wx>  = (t)  = n^  N2(t)  = n2> 


The  numerator  of  the  above  expression  is  equivalent  to 


(A. 8)  Pr{N  (T1  )=  x-n!|N  (t)  = n ; N (t)  = n } 

2 (n;+l)  11  1 


/ni  V "2  A 
nl+1  \ni+1A  x_njV 


x+l 


' x+l  / 


(from  Theorem  1) 


Vnl 

Vn2"X 


From  Corollary  1.2,  we  have 


(A. 9)  Pr{N1(W}{)  = nJ_|N1(t) 

gc;) 


= ty  N2(t2)  = n2} 


15C 


From  Eqs.  (A. 8)  and  (A. 9)  we  obtain  Eq.  (4.7) 

Q.E.D. 

We  now  extend  the  above  results  to  the  case  when  there  are 
more  than  two  independent  sources  generating  the  arrivals.  Let 
{N.(t),  t>0},  i = l,2,...,p  be  the  independent  Poisson  processes 
with  parameter  A^»  respectively.  We  use  the  fold owing  notation 
similar  to  those  used  in  the  previous  case. 

N(t)  = f N . ( t ) , t>0. 

i=l  1 

T,1  = the  time  of  the  kth  arrival  from  process  i,  i=l,2,...,p 

k 

and  k > 1,  ( Tq  — 0 , i-l,2,..»,p). 

W = the  time  of  the  xth  arrival  in  the  total  arrival  process 
x 

(N(t),  t>0}. 

Similar  to  theorem  1,  we  have 

Theorem  2.  Let  (nj,  {n|}  i=l,2,...,p  be  non-negative  integers 
such  that 

(a)  f n.  > 0,  l n!  > 0 
i=.l  i=l 


and 


(b) 


n!  < ? n! 


i = 1,2, 


,p; 


Then 
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(A. 10)  PrtN.Or1 *  ) = n.\  j-1,2,. 

3 n.'  3 

1 


, . ,p;  j i i ! WR(t)  = k = 1,2, . . . ,p) 


nj\ 


n . . !1  \ n ! 

i 0 = 1 x 1 


> n.’  / y n. 
3-1  3 / M 3 

l E 

\ l ni 

\j=i  3‘ 


Proof : Similar  to  the  proof  of  theorem  1.2. 

Parallel  to  Corollaries  1.1  and  1.2  we  have  the  following. 
Corollary  2.1.  Under  the  assumptions  of  Theorem  2, 


Pr{  N 


.(VJ  ) = n'.,  i = 1,2, . . . ,p  |N . (t)  = n.,  i=l,2,. 

1 (n^w^ +. . .n^)  x 


• ,P> 


P /“i 

n ( , 

i=l  \ 1/ 


Proof:  Similar  to  the  proof  of  Corollary  1.1. 

Corollary  2.2.  Let  I(s)  denote  the  index  of  the  process  from  which 
the  next  arrival  occurs  after  time  s (0  <^s  < t),  then 

Pr(l(L'x)  = i|N^(Wx)  = n'^,  M^(t)  = n<f  k=l,2,...,p} 


n . - n . 

l l 


I n - I n ' 

i=l  3 3=1  3 


l = 1>2, • . . ,p. 
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APPENDIX  3 

SIMPLIFICATION  OF  COMBINATORIAL  E) PRI 3SIONS 
1.  Derivation  of  Eg.  (4.19) 

We  first  simplify  Eq.  (4.16).  To  guarantee  that  all  values  within 
the  combinatorial  expressions  are  non-negative,  the  range  of  dj^'t2,'C3^ 
satisfies:  0 < d^(t0,t3'  < z^t^)  - d^t^tg).  Expanding  the  combina- 

torial expressions  and  simplifying  the  summation  over  d^tjjt^),  we 
obtain 


v . . i -<vv  } 


[X 


r 0 . WV  " Z0(tl)  + d0(tl,t2) 

= U - Aj/AQJ 


Z0(tl)  " d0(tl’t2)  ^ 
\~ 

/ 


( z0(tl)  " d0(tl,t2) 


d-<VV  - o U V W 


..CxX]d°(t2,t8V-W°(tl>  ’ d°(tl’t2)  ’ dj(t2’t3>  \ 


J J 


= [I-AJ/Aq] 


0,  d0(t2,t3)  “ z0(tl)  + d0(tl,t2) 


Similarly,  to  simplify  Eq.  (4.17),  we  find  that  the  ranges  of 
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j 
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d.  (t2,t3>  and  dQ(t2,t3)  yielding  non-negative  values  within  th« 


combinatorial  expressions  are 


(i)  d < d°(t2,t3)  ^ " do(ti,1:2^  + d’ 


and  ( i i ) do(t2’t3)  - Z0(tl)  " <Vtl’t2)  + d' 


Expanding  the  combinatorial  expressions  and  simplifying  we  obtain 


Z0(tl)"d0(Vt2)+d  j / C \ 

j \(d0(t2,t3)  " Z0(tl}  + 10(tl,t2)j 

! V a / 


dj(t2’t3)=d 


/W  - d0(tl,t2) 


dj  Ct 2 * C3  > 


r o dj(t2’t3)+d 

• tX-j/x0J 


0 d0(t23t3)_dn(t2,t3)'d  ) 

a - 0 2 3 3 2 3 f 


d0(t2’V  " z0(tl)  ‘ d0(tl,t2) 


0 d 0 d0^2’^3^  ~ Z0^1^  d + do^l’^2^ 

[x:/x  ] [i-x:/x  ] 0 1 d 01 

jo  3 u 


zQ(ti)  - dQ(t1,t2) 


d-<VV  = 0 


,0(tl)  " d0(tl’t2) 


dj(w 


d/xjdi<t2’t3>ci-x?/x„]Zo<tl)  ■ do<tl’t2> ' dJ<t2’t3V> 


^ V'.-: 


Substituting  these  simiplif ications  and  taking  the  lirr.it  as 


t -»•  00 , we  obtain  Eq . (4.19). 

2.  Derivation  of  Eq.  4.24) 

0.  . 

From  Eqs . (4.21)  and  (4.23)  we  obtain  the  ranges  of  dd^.tj) 
d (t^tj)  such  that  the  values  within  the  combinatorial  expressions 

nonnegative. 

(i)  d < d”(t1,t2)  < z0(tx)  + d°(t1,t2)  + d 

and  (ii)  d0^tl,t2^  >_  z^t^)  + d- 

Upon  expanding  the  combinatorial  expressions  and  simplifying  m 
the  manner  similar  to  that  in  the  previous  case  we  obtain  Eq.  (4.24) 


and 

are 


BIBLIOGRAPHY 


[1]  S.G.  Allen  and  D.A.  D'Esopo,  "An  Ordering  Policy  for  Repairable 
Stock  Items,"  Opns . Res . 1G,  p.  669-774  (1968). 

[2]  P.  Billingsley.  Convergence  of  Probability  Measures,  Wiley,  Now 
York  (1968). 

[3]  A.J.  Clark,  "An  Informal  Survey  of  Multi-Echelon  Inventory  Theory", 
Naval  Res.  Log.  Quart.,  19,  p.  521-650  (1972). 

[4]  G.F.  Feeney  and  C.C.  Sherbrooke,  "A  System  Approach  to  Base 
Stockage  of  Recoverable  Items",  The  Rand  Corporation,  RM-4720-PR 
(DDC  No.  AD  626744).  (December  1965). 

[5]  G.F.  Feeney  and  C.C.  Sherbrooke,  "The  (s-1,  s)  Inventory  Policy 
Under  Compound  Poisson  Demand,"  Hgmt.  Sc . , 12,  p.  391-411  (1966). 

[6]  D.  Gross  and  C.ti.  Harris,  Fundamentals  of  Queueing  Theory, 

John  Wiley  and  Sons  (1974). 

[7]  G.  hadiey  and  T.M.  Whit in.  Analysis  of  Inventory  Systems, 

Prentice- Hall,  Inc.,  Englewood  Cliffs,  New  Jersey  (1963). 

[8]  D.L.  Isaacson  and  R.W.  Madsen,  Markov  Chains  Theory  and 
Applications  , John  Wiley  (1375). 

[9]  W.K.  Kruse  and  A.J.  Kaplan,  "On  a Paper  by  Simon",  Opns . Res . , 

21,  p.  1318-1322  (1973). 

[10]  Francois  bureau,  "A  Queueing  Theoretic  Analysis  of  Logistics 
Repair  Models  with  Spare  Units",  Technical  Report  No.  167, 

Department  of  Operations  Research  and  Department  of  Statistics, 
Stanford  University,  Stanford,  California  (December,  1974). 

[11]  J.A.  Muckstadt,  "A  Model  for  a Multi-Item,  Multi-Echelon, 

Multi- Indenture  Inventory  System,"  Mgmt.  Sc.,  20,  p.  472-481, 

(1973). 

[12]  J.A.  Muckstadt,  "Some  Approximations  in  Multi-Item,  Multi-Echelon 
Inventory  Systems  for  Recoverable  Items",  The  Rand  Corporation, 

F-5763  (December  1976). 

[13]  F.R.  Richards,  "A  Repairable  Item  Inventory  System",  Opns.  Res., 

24,  p.  118-130  (1976). 

[14]  S.M.  Ross,  Applied  Probability  Models  With  Optimization  Applications, 
Holden- Day  (1970). 

[15]  D.A.  Schrady,  "A  Deterministic  Inventory  Model  for  Repairable 
Items,"  Naval  Res.  Log.  Quart.  14,  p.  391-398  (1967). 


[16]  D.A.  Schrady,  "Mathematical  Models  of  the  Repairable  Item 
Inventory  System" , Report  No.  NFS  55S07081A,  U.S.  Naval 
Postgraduate  School,  Monterey,  California  (August  1967). 

[17]  C.C.  Sherbrooke,  "METRIC:  A Mult.i-Echelcn  Technique  for 

Recoverable  Item  Control",  Opns.  Res.,  16,  p.  122-141  (1968). 

[18]  C.C.  Sherbrooke,  "Discrete  Compound  Poisson  Processes  and 
Tables  of  Geometric  Poisson  Distribution",  Tr.c  Rand  Corporation, 
RM- 4831- PR  (July  1966). 

[19]  R.M.  Simon,  "The  Uniform  Distribution  of  Inventory  Fosition  tor 
Continuous  Review  (s,Q)  Policies",  P-3938,  The  Rand  Corp., 

Santa  Monica,  California  (1968). 

[20]  R.M.  Simon,  "Stationary  Properties  of  a Two-Echelon  Inventory 
Model  for  Low-Demand  Items",  Opns.  Res.  19,  p.  761-773  (1971). 

[21]  R.M.  Simon, and  D.A.  D’Esopo,  "Comments  on  a Paper  by  S.G.  Allen 

and  D.A.  D'Esopo:  An  Ordering  Policy  for  Repairable  Stock  Items", 

Opns.  Res.,  19,  p.  986-988  (1971). 

[22]  B.D.  Sivaslian,  "A  Continuous  Review  (s,S)  Inventory  System  with 
Arbitrary  Interarrival  Distribution  between  Unit  Demand'  , Opns. 
Res.  22,  p.  65-71  (1974). 

[23]  H.C.  Tijms,  Analysis  of  _(s, S)  Inventory  Models,  Mathematisch 

centrum,  Amsterdam  (1972). 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

RKAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1 REPORT  NUMBER 

341 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Snbtltle) 

AN  ANALYSIS  OP  A TWO-ECHELON  INVENTORY  SYSTEM  FOR 
RECOVERABLE  ITEMS p 

5.  TYPE  OF  REPORT  a PERIOD  COVERED 

d.  Technical  Repaijt 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfa) 

Kripa/shanker  v 

8.  CONTRACT  OR  GRANT  NUMBERf*) 

,C  ‘ 

-Ld'N0Q14-75-C-1172  "1 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

School  of  Operations  Research  and  Industrial 
Engineering,  College  of  Engineering,  Cornell 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

Task  NR  042-335 

. 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Sponsoring  Military  Activity 
Statistics  and  Probability  Program 

,2- 

J ‘ y — 

n.  NUMBER  OF  P«fi£S 

165  . IL  //?. 

14  MONlVoRlNG  Att^NC?  NAMfe  i Cc*ifrofj/ni 

IS.  SECURITY  ClXSS.  (ofay,  r.pofl^f 

Unclassified  \ — 

15  •.  DECLASSI  FI  CATION/ DOWNGRADING 
SCHEDULE 

y 


16.  DISTRIBUTION  STATEMENT  (ot  thlm  Heport) 


Approved  for  public  release;  distribution  unlimited. 


17  DISTRIBUTION  STATEMENT  (of  the  abstract  antorad  In  Block  20,  It  dltfarent  trom  Report) 


10.  SUH^- EMENT  ARY  NOTES 


19.  KEY  WOHOS  (Contlnua  on  teveraa  aids  If  necessary  and  Identify  by  block  number) 


Inventory  Theory 
Single  Location  System 
Two-Echelon  System  ’ 


Arrival  Sequences 
Unit  Model 
Batch  Model 


2(L  A TR  ACT  rCurrl fmjm  on  rarer* a aide  If  neceeeary  and  Identify  by  block  number) 

In  this  dissertation,'  we  present;  an.  analysis  of  continuous  review  models 
of  i two- o jhelon  inventory  system  for  recoverable  items.  The  system  consists 
of  a depot  anc  a set  of  bases.  Primary  demands  occur  at  the  bases  for  one  or 
;everal  units  at  a time.  It  is  assumed  that  demands  arrive  in  a Poisson  manner 
boon  arrival  of  a demand  for  certain  units,  a like  number  of  failed  units  are 
turned  in  it  the  base.  An  inspection  of  the  failed  units  is  carried  out  to 
decide  whether  the  units  will  be  repaired  at  the  base  or  at  the  depot  or  will 


DO  1473 


EDITION  OF  * MOV  Si  IS  OBSOLETE 


Unclassified  /, 


f 


SECURITY  CLASSIFICATION  or  THIS  PAGE  fMNar.  Data  Enl«r»rf> 


'4-  W 


SECURITY  CLASSIFICATION  OF  THIS  PAGE(»Ti«.  Dmlm  Enfnd) 


. 

• * j 

—sbe  removed  from  the  system  in  case  repair  is  not  economical.  The  bases  use 
an  (s-1,  s)  policy  tor  procurement  of  serviceable  units  from  the  depot,  arid 
the  depot  uses  an  (s,S)  policy  to  procure  from  the  external  supplier. 

Demands  in  an  out-of-stock  situation  are  backlogged.  It  is  assumed  that  all  th 
locations  have  infinite  repair  capacities  and  repair  and  procurement  lead 
times  are  constant. 

A common  problem  in  inventory  management  is  to  specify  the  policy  parameter 
that  will  minimize  expected  cost  per  unit  time  for  operating  the  systejn 
subject  to  constraints  of  certain  performance  measures.  To  formulate  such  a 
problem  we  must  find  the  stationary  distributions  for  inventory  position, 
on-’nand  inventory,  backorders  and  in-repair  inventory.  Our  main  objective 
is  to  find  exact  expressions  for  these  distributions^-^' 

The  investigation  begins  with  an  extensive  analysis  of  a single  location 
system.  The  procurement  policy  is  a continuous  review  \(s,S)  policy.  The 
inter-arrival  times  between  successive  requisitions  are  independent  and  iden- 
tically distributed  random  variables.  The  system  experiences  two  types  of 
demands  - recoverable  and  non-recoverable.  The  two  demand  processes  may  be 
independent  or  dependent.  For  the  inspection  of  failed  units,  two  models  - 
batch  and  unit  - are  considered.  In  the  batch  model,  the  entire  batch  is 
either  recoverable  or  non-recoverable,  whereas,  in  the  unit  model  each  unit  in 
a batch  is  inspected  independently,  the  special  cases  of  compound  Poisson 
demands,  (s,nQ)  procurement  policy,  complete  recoverability  and  complete 
non-recoverability  are  also  considered. 

For  the  two-echelon  system  we  first  consider  the  case  where  demands 
at  the  bases  occur  for  a single  unit  at  a time.  The  approach  is  then  applied 
to  a general  situation  where  demands  at  the  bases  are  random.  Both  the 
batch  and  unit  inspection  models  are  considered.  For  the  case  when  there  are 
no  condemnations  of  the  item,  results  are  compared  with  the  METRIC  model . The 
METRIC  model  provides  a simple  but  approximate  expression  for  the  probability 
distribuat ion  of  system  backorders.  The  comparison  indicates  that  there  is  a 
considerable  discrepancy  between  the  METRIC  results  and  our  results  when  the 
depot  spare  stock  is  low  or  when  a major  proportion  of  the  repair  is  done 
at  the  depot. 
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